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Effects of Isomerism on Infrared Absorption Spectra 
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The absorption spectra of four organic cyanides, two 
organic thiocyanates, and the corresponding isomers have 
been studied in the region from 1.0u to 6.5u. The intense 
band appearing near 4.4u in the normal cyanides appears 
at a longer wave-length and with doubled intensity in the 
isocyanides, the variation in wave-length and intensity 
being nearly constant for all compounds. In the case of 
the thiocyanates the variation in wave-length of the 4.6u 


band for the two configurations is less than in the cyanides, 
whereas the intensity differences are greater. No other 
bands in these organic compounds were appreciably 
affected by isomerism. In the absorption spectrum of 
aqueous solutions of hydrocyanic gas there appears an 
intense band at 4.75 and also a less intense band at 4.95. 
The band at 4.95 can be attributed to a small amount of 
the isomer HNC. 





OBLENTZ! in his studies of the infrared 

absorption of isomeric organic compounds 
found that differences in the structure of the 
molecules are accompanied by differences in 
spectra, the effects being more pronounced in the 
region of longer wave-lengths. In many cases the 
spectra of isomers are identical up to about 6 u, 
at which point marked shifts begin to appear. 
Other studies? made by Weniger on alcohols, 
acids, and esters show that the greater the 
difference in structure between isomers, the 
greater the difference in their spectra. 

Recently the authors have investigated the 
absorption spectra of thirteen cyanides and 
thiocyanates.* For each of these compounds a 
well-defined band was observed in the region 
between 4.4 uw and 4.9 yw, the exact position of the 
band being different for different compounds. The 
absorption in this region was attributed to funda- 
mental frequencies of the cyanide and _ thio- 

1'W. W. Coblentz, Investigations of Infrared Spectra 1-2, 


ag Institution of Washington, 35, 1905), pp. 
? W. Weniger, Phys. Rev. 31, 388 (1910). 
(1935) Gordy and D. Williams, J. Chem. Phys. 3, 664 
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cyanate molecules in which a CN vibration is 
active. It was the purpose of the present investi- 
gation to ascertain the effects of isomerism on 
this characteristic vibration. The isocyanides 
are different from isomers previously studied in 
that the carbon is divalent in the isocyanides and 
tetravalent in the normal cyanides. 

The experimental method was the same as that 
used in the previous investigation.* A Hilger 
spectrometer with a fluorite prism was used for 
the entire region from 1.0 u to 6.5 uw. Cell windows 
were of fluorite, and the thickness of the absorb- 
ing layer was 0.03 mm for all materials studied 
unless otherwise specified. 

Absorption curves for four cyanides and the 
corresponding isocyanides are shown in Fig. 1. 
The solid lines represent the percent transmission 
of the cyanides; the broken lines, that of the 
isocyanides. The band in the region of 3.25 u is 
characteristic of the group associated with the 
cyanide group and is found in other methyl, 
ethyl, butyl, and phenyl compounds.?: * Within 
the limits of experimental error the position and 


4E. K. Plyler and T. Burdine, Phys. Rev 35, 605. 


(1930). 
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Fic. 1. (a) Transmission of methyl cyanide and methyl isocyanide. (b) Transmission of ethyl 
cyanide and isocyanide. (c) Transmission of butyl cyanide and isocyanide. (d) Transmission of 


intensity of this band are the same for both iso- 
meric forms, the band arising from the methyl 
group being a possible exception. The character- 
istic cyanide band in the region of 4.6 u is shifted 
to the longer wave-lengths by about 0.2 u and 


phenyl cyanide, phenyl isocyanide, and phenyl thiocyanate 2.74 to 5.1u. 


TABLE I. 








BAND POSITIONS 


Cyanide 


VARIATION IN RATIO OF 
WaveE-LENGTH BAND INTENSITY 
Isocyanide to 





Methyl 3.25 


4.38 


Ethyl 3.22 
4.43 


Butyl 3.30 
4.45 


Phenyl 3.16 
4.48 


Thio- 

cyanate 

Methyl 3.32 
4.64 


Ethyl 3.28 
4.63 


Isocyanide AX Cyanide 
3.25u 0.00u 
4.58 0.20 
3.25 0.03 
4.63 0.20 
3.30 0.00 
4.65 0.20 
3.14 —0.02 
4.71 0.23 

I sothiocya- 
nate 
3.32 0.00 
4.70 0.06 
3.29 0.01 
4.75 0.12 
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its intensity is approximately doubled in the 
isocyanides. Furthermore, these variations in 
wave-length and intensity are approximately 
constant for all the compounds studied. The 
shift in wave-length is in agreement with the 
chemical data which indicate divalent carbon in 
the isocyanides and hence weaker binding forces. 
Table I gives a numerical comparison of band 
positions and intensities for corresponding com- 
pounds. 

The earlier work of Coblentz! on two thiocya- 
nates and the corresponding mustard oils was re- 
peated for purposes of intensity comparison and 
for more accurate determination of band position. 
The absorption curves for methyl and ethyl thio- 
cyanates and their isomers are given in Fig. 2. 
As may be seen from the figure, the shifts due to 
isomerism are different for the two and in both 
cases are less than the corresponding shift for the 
cyanides. On the other hand, the intensity varia- 
tion is much greater than in the case of the cya- 
nides. As before there is no apparent shift in the 
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_ Fic. 2. (a) Transmission of methyl thiocyanate and 
isothiocyanate. (b) Transmission of ethyl thiocyanate and 
isothiocyanate. Region: 2.74 to 5.14. 


band characteristic of the associated methyl and 
ethyl groups. The absorption of phenyl! thiocya- 
nate was measured and is given in Fig. 1 with the 
other phenyl compounds. 

The fact that the band in the 4.64 region 
varies while the 3.25 « band remains constant in- 
dicates that both these bands depend primarily 
on the structure of certain’ groups rather than 
upon the structure of the molecule as a whole. 
Although a survey was made of the whole region 
from 1.04 to 6.54, no bands were observed 
which could be classified either as harmonics of 
the 4.6 « band or as combinations involving this 
characteristic frequency. This is further indica- 
tion of a simple type of binding between the 
carbon and nitrogen atoms, which is only slightly 
modified by associated groups in the molecule. 

Nef> has proposed the two configurations, 
HCN and HNC, for hydrocyanic gas. If both 
configurations be present in a given sample, it 
should be possible to observe two bands in the 
region of 4.7 u. In the previous investigations, 


‘J. V. Nef, J. Am. Chem. Soc. 26, 1549 (1904). 
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lic. 3. The ratio of the transmission of hydrocyanic acid 
solution to that of an equal thickness of water for two 
cell thicknesses. Region: 4.1 to 5.5. 


the authors obtained an absorption curve show- 
ing a band at 4.75 uw characteristic of HCN in 
aqueous solution. In this absorption curve there 
is also a slight depression at 4.95 yw. For greater 
cell thickness this depression becomes deeper and 
can be interpreted as arising from the isomer 
HNC. In Fig. 3 are given the transmission curves 
of hydrocyanic acid for two cell thicknesses. On 
account of the different absorption coefficients 
for CN and NC compounds it is impossible to 
ascertain very accurately the relative abundance 
of the isomers, but it is evident from the figure 
that the amount of HNC present is very small, 
as is to be expected from chemical theory.® As- 
suming the absorption coefficient of HNC to be 
twice that of HCN, as is the case for the other 
isomers studied, we estimate the concentration of 
HNC to be from two to four percent that of HCN. 

Barker’ and Choi and Barker® have studied 
the absorption of gaseous HCN, but have re- 
ported no absorption in the region of 4.9 uw. This 
region may have been beyond the range of their 
investigation, as their published curves for this 
part of the spectrum extend only from 4.68 u 
to 4.86 uw. A second possibility is that the isomer 
HNC may be more abundant in aqueous solution 
than in the gas. 

The authors wish to express their thanks to 
Professor E. K. Plyler for his helpful discussion of 
the manuscript and to Mr. G. C. Kyker of the 
Chemistry Department, who prepared the chem- 


icals and offered valuable suggestions. 


6 J. Schmidt, Text Book of Organic Chemistry, (D. Van 
Nostrand Co., 1932), p. 323. 
7E. F. Barker, Phys. Rev. 23, 200 (1924). 


8K. N. Choi and E. F. Barker, Phys. Rev. 42, 777 
(1932). 











FEBRUARY, 1936 















JOURNAL OF CHEMICAL PHYSICS 


Raman Spectrum of Crystalline Ammonium Chloride* 


Crystalline NH,Cl was studied in its room temperature and high temperature forms. Cy- 
lindrical crystals were grown by sublimation in a special furnace. At room temperature vi- 
brational frequencies were observed at 3146, 3041, 2824, 2010, 1768, 1709, 1407 and subsidiary 
frequencies at 3300, 3232 cm—. In the high temperature modification vibrational frequencies 
were observed at 3102, 2801, 1931, 1655, 1400 cm=4. 
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F, T. Hotmes, Rouss Physical Laboratory, University of Virginia 
(Received September 9, 1935) 





RYSTALLINE ammonium salts have many 

curious physical properties,! some of which 
are ascribable to characteristics of the ammonium 
radical not yet fully understood. This work was 
undertaken in the hope that light would be 
thrown on some of these characteristics by more 
complete observations of the Raman effect in 
NH.Cl. 

The data obtained have been combined with 
the infrared observations of previous workers and 
interpreted with the help of two reasonable 
assumptions. Probably the most important result 
is the indication that the NH, group in crystalline 
NH.CI either has not the geometrical structure 
of a tetrahedron of cubic symmetry or else has 
this structure but is strongly perturbed in a 
nongeometrical sense. 























GROWING CRYSTALS BY SUBLIMATION 






A major difficulty was encountered in pre- 
paring a relatively pure and optically clear 
crystal of suitable size in the room temperature 
modification of this substance. Many attempts to 
grow crystals in the high temperature form and 
cool them through the 184°C transition! pro- 
duced merely mediocre crystals. However, it was 
found that an optically clear cylindrical crystal 
3 cm long could be grown by sublimation, at a 
temperature below? 184°C, in about a month. 

The essentials of the glass-walled furnace 
constructed for this purpose are shown in Fig. 1. 
The heating wire was wound directly on tube 1, 



















* Portion of a dissertation presented for the degree of 
Doctor of Philosophy at Yale University. 

1 Most of the peculiarities of the ammonium halides are 
mentioned by A. C. Menzies and H. R. Mills, Proc. Roy. 
Soc. A148, 407 (1935). 

UB 18 International Critical Tables, Vol. 4, p. 7. 
a ?R. Pohlman, Zeits. f. Physik 79, 394 (1932). 
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this assembly being separated from tube 2 by a 
coarse winding of asbestos string. Annular 
windings of asbestos, as at A, separated tubes 2 
and 3. Heating current was obtained from a 
commercial power line, the voltage fluctuations 
being sufficiently reduced by the use of “‘ballast 
tubes.” 

The tube 7 was half filled with NH,Cl (pre- 
viously sublimed several times in vacuum), 
evacuated, heated for some time, and sealed off. 


It was inserted in the furnace as shown and, 


moved toward the left. After the rounded end of 
T had passed M, the region of maximum 
temperature, a crystal C started to form. The 
rate of motion was then determined by the 
appearance of the surface of C, minute facets 
indicating dangerously fast growth. 

For high temperature observations a similar 
furnace was inserted in the arc chamber. Es- 
sential modifications consisted in replacing 2 by a 
tube of Corning Red Ultra glass* and sealing a 
flat window on the end of 7. Crystals were grown 
and then maintained at temperatures well above 
184°C during exposures. 


Room TEMPERATURE RESULTS 


Spectra of the room temperature modification 
were obtained by using Hilger constant deviation 
type instruments fitted with extra-dense flint 
prisms and long and short focus cameras. 
Commercial mercury and helium arcs were used 


3 L. J. Buttolph, Rev. Sci. Inst. 1, 650 (1930). 
















RAMAN SPECTRUM OF 


with standard filters in the arc chamber described 
by Buttolph.’ 

The vibrational frequencies in cm, with 
estimated relative visual intensities, are:3146, 15; 
3041, 20; 2824, 7; 2010, 2; 1768, 2; 1709, 5; 1407, 
2 (average deviation 4 cm or less). The first 
three are inherently of intermediate width, the 
next two are broad‘ and without cores, 1709 is 
sharp, and 1407 is sharp on its low wave-length 
side, as is the wider line 2824 (and possibly 3146). 

The infrared absorption frequencies observed 
by Wilberg® at approximately 2350 and 2250 
cm were not found, although 2350 (He 3888A), 
if present, might have been masked by an 
observed Raman line assigned to 2824 (He 
3819A). 

Considering the widths of most lines, these 
results agree satisfactorily with the less complete 
results of Schaefer, Matossi, and Aderhold,® 
Kastler,? and Menzies and Mills,‘ and, except 
for 2824 and 1407, with weighted means of the 
infrared results of Wilberg® and Pohlman.? One 
should note, however, that 1709 is inactive in 
infrared absorption at room temperature. 

Subsidiary frequencies were observed at 3300 
+6 and 3232+5 cm. These may correspond to 
fine-structure lines observed by Pohlman at 
—80°C. 


HiGH TEMPERATURE RESULTS 


The spectrum of the high temperature modifi- 
cation was obtained only at low dispersion with 
helium excitation. 

Observed frequencies in cm are: 3102, 30; 


2801, 4; 1931, 2; 1655, 2; 1400, 2 (average 
deviation 6 cm or less). The first, a very wide 
line, corresponds to one or both components of 
the high frequency pair which are well resolved in 
crystals at room temperature with this dis- 
persion. On the basis of their intensities and 


*Compare A. C. Menzies and H. R. Mills, Proc. Roy. 
Soc. A148, 407 (1935). 

a Wilberg, Zeits. f. Physik 64, 304 (1930). See also O. 
Reinkober, ibid. 5, 192 (1921). 

°C. Schaefer, F. Matossi, and H. Aderhold, Zeits. f. 
Physik 65, 289 (1930). 

7A, Kastler, Proc.-verb. 
Bordeaux (1931-1932). 
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w weak; 
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AMMONIUM CHLORIDE 89 
widths 1931 and 1655 correspond to 2010 and 
1768. Apparently the room temperature 1709 
cm! frequency is much weaker or has shifted in 
the direction of shorter wave-lengths to an 
unobservable position. No asymmetries were 
observed. 


DISCUSSION 


The 1407 cm~ frequency seems to be an 
unresolved doublet.?: > This view is supported by 
the fact that the corresponding NH,Br frequency 
is definitely double at low temperature.? Thus 
ten vibrational frequencies can be said to have 
been observed at room temperature in crystalline 
NH.Cl. 

Two assumptions will be applied throughout 
the rest of this paper: (1) That of most authors, 
namely, that all frequencies in this spectral 
region are produced by NH,'; (2) that no NH,* 
frequency in crystalline NH,Cl lies below*® 1407 
cm, These assumptions necessitate regarding at 
least one of the observed frequencies as due to a 
harmonic, a combination, or a perturbation. 
Consideration of all relevant data indicates that 
in NH,Cl 2824 may be the harmonic of 1407 
and 3146 may be a combination frequency, 
possibly due to 1407 plus 1709. 

The structure usually assumed for NH,* 
that of a tetrahedron of cubic symmetry. Using 
this model Hettich® and Menzies and Mills‘ have 
presented explanations of the complex macro- 
scopic phenomena in ammonium halide crystals. 
Also, Pohlman? has deduced from observed fine 
structure (without allowing for the possibility 
of vibration-rotation interaction’®) a value for 
the N—H distance which agrees with that 
calculated by Laschkarew and Usyskin"™ from 
electron diffraction patterns. However, this 
structure is spectroscopically unsatisfactory un- 
less seriously perturbed in a nongeometrical 
sense, since without such a perturbation only 
four fundamental vibrational frequencies could 
be observed.” Unfortunately the available spec- 


8 See, however, A. da Silveira, Comptes rendus 195, 521 
(1932) and Wilberg, reference 5. 

9 A. Hettich, Zeits. f. physik. Chemie A168, 353 (1934). 

10D. M. Dennison and M. Johnston, Phys. Rev. 47, 93 
raat E. Teller and L. Tisza, Zeits. f. Physik 73, 791 
(1932). 

iW. E. Laschkarew and I. D. Usyskin, Zeits. f. Physik 
85, 618 (1933). 

12D, M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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troscopic data do not strongly indicate the 
correctness of any particular geometrical con- 
figuration. 

Comparison of the Raman spectra of the room 
temperature and high temperature forms of 
NH,Cl suggests various possibilities. Of these, 
the simplest are that the frequencies occurring at 
room temperature at 2010 and 1768 cm™ origi- 
nate in similar types of vibration and that the 
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properties of NH,* undergo only minor changes 
at the high temperature transition point. 

The author is indebted to the physics de- 
partments at Lehigh University and the Uni- 
versity of Virginia for laboratory facilities, and 
wishes to thank Dr. J. E. Rosenthal of Columbia 
University for discussions, and Professor W. W. 
Watson of Yale University for his interest and 
advice. 
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The infrared absorption spectra of mixtures of water 
and acetic propionic and of propionic anhydrides have 
been measured in the region from 5u to 6.4u. The intensity 
of the bands was found to decrease with time. By using 
the relative intensity of the anhydride band the reaction 
rate was calculated for three concentrations. The reactions 
were unimolecular but the value of the reaction constant 





T# E reaction velocity of acetic anhydride has 
been measured! recently by means of its 
infrared absorption spectra. The results showed 
that the reaction constant increased with tem- 
perature and that the value of the constant de- 
creased in a given set of observations. The latter 
result may have been produced by the method of 
calculation. In order to study further the effect of 
concentration the work has been extended to 
propionic and acetic propionic anhydrides. 
Before the absorption spectra of different con- 
centrations were measured the region from 5 u 
to 6.5 wu was studied. It was found that the ab- 
sorption spectra of acetic and of propionic acid 
were similar in this region and that the bands, 
which were at about 5.7 wu, varied in position by 
about 0.03 u. The propionic and acetic propionic 
anhydrides were studied and it was found that 
there were two absorption bands at 5.48 w and 
5.66 uw, and 5.50 w and 5.68 yu, respectively. When 
: cells of about 0.01 mm in thickness were used, 
i the two bands could be separated by the fluorite 


1 Plyler and Barr, J. Chem. Phys. 3, 679 (1935). 
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The Reaction Rates of Propionic and Acetic Propionic Anhydrides 


E. S. BARR AND E. K. PLyLer, Department of Physics, University of North Carolina 
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was different for each concentration, its value being 
approximately proportional to the concentration of the 
water. The value of the reaction constant of acetic pro- 
pionic anhydride was about twice the value of the constant 
for propionic anhydride. A table is given of the values of 
the reaction constants for the different concentrations. 
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prism. The two anhydrides have similar spectra, 
as can be seen by the absorption curves shown in 
Fig. 1. The fact that the 5.504 band can be 
separated from the other anhydride band at 5.68 
pw and the acid bands at 5.7 uw made it suitable for 
the study of reaction rates. 

In the study of the reaction velocity, cells were 
made of three concentrations. The first solution 
had equivalent amounts of water and anyhdride, 
the second two equivalent parts of water to one 
of anhydride, and the third three of water and one 
of anhydride. The actual concentrations of water 
and anhydride are given in Table I. The same 
concentrations were used for each anhydride. The 
reactions were allowed to progress in the absorp- 
tion cell. After preliminary trials it was found 
that the reactions were slow and that the absorp- 
tion curves should be measured at time intervals 
of from one to two hours. In this way six to eight 
sets of observations were obtained for each con- 
centration. In Fig. 2 are shown the curves for 
the three concentrations of propionic anhydride. 
The bands in the section on the left are for 
equivalent amounts and the time of observation 
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Fic. 1. The absorption spectra of propionic and acetic 
propionic anhydrides in the region from 4.8u to 6.1n. 
The cell thickness was 0.01 mm. 


is given beside each curve. The fourth curve in 
the first group was placed inaccurately on the 
drawing. It should be shifted 0.06 » to the shorter 
wave-lengths. The middle set of observations are 
for two parts water and one part anhydride. The 
lowest curve in each set is for the shortest interval 
of time after mixing, and shows the 5.48 » band 
with the greatest intensity: These curves show 
that in any concentration the absorption in the 
region of 5.54 decreases with time. Also this 
decrease in intensity with time is greater when 
the concentration of the anhydride is smaller. 
Three sets of observations were also made for 
the acetic propionic anhydride. The absorption 
spectra were similar to those shown in Fig. 2. 
It was found by considering the reaction as uni- 
molecular, that approximately consistent values 
were obtained for the reaction constant. In order 
to find the concentration of the anhydride at any 
time, it was necessary to know the relative in- 
tensity of the band at the beginning of the reac- 
tion and also the relative intensity at the time of 
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Fig. 2. The absorption of propionic anhydride and water 
after different intervals of time. The different sets of 
curves represent mixtures of 1, 2 and 3 parts water to one 
part anhydride. 


TABLE I. The reaction constants of propionic and acetic pro- 
pionic anhydrides at 27°C, for different concentrations. 
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the observation. The relative intensity of the 
bands for different times of observation was 
found by determining the area under the absorp- 
tion curve. The relative intensity at zero time 
was found by extrapolation. 

In Table I the values of the reaction constant 
are given for the reaction rate of propionic and 
acetic propionic anhydrides. 

The constant for acetic propionic anhydride 
is roughly twice that of propionic anhydride for 
each concentration. For a given concentration 
the values of k tend to decrease with time. 
However, this effect is not so great as found in the 
case of acetic anhydride. The change in the value 


2 A more complete description of the method of calcula- 
tion can be found in the paper of reference 1. 
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PLYLER 


of the constant may not be a real property of the 
reaction. The method of measurement and calcu- 
lation may produce the change. All measurements 
were made at 27°C. 

When the concentration of water in the solu- 
tion was increased the reaction rate increased. 
It was found that the value of k for the different 
concentrations was approximately proportional 
to the equivalent amounts of water in the 
mixture. 

On account of the anhydride band being less 
overlapped by the acid band, it is believed that 
the accuracy is greater than in the previous 
study of acetic anhydride. The percent error in 
the results are probably between 5 and 10 per- 
cent. 
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The Near Infrared Absorption Spectra of Acid and Basic Solutions 


E. S. BARR AND E. K. PLYLER, Physics Department, University of North Carolina 


The infrared absorption spectra of H2O and aqueous solutions of HCl, HBr, NaOH, KOH, 
ZnBr2, ZnClz and NasCO; have been measured from 1.5u to 2.8u. Bands were observed at 
2.30u and 2.45y in hydroxides, at 2.30u and 2.55u in acids, and at 2.30u, 2.454 and 2.55 in 
hydrolyzing salts. The band at 2.30u is produced by hydration and 2.454 band by the OH ion. 


TUDIES! have already been made on the in- 
frared absorption of acid and basic solutions 
from ip to 7p. In the study on acid solutions a 
broad band was observed at 2.44 and in the 
hydroxide solutions at 2.34. When the hydrolyz- 
ing salts were studied, the two bands could not 
be separated. With greater resolution it would be 
possible to separate the bands in the 2.4y region. 
A large type Hilger infrared spectrometer, with 

a quartz prism, was used for the resolving instru- 
ment. The deflections were increased by using a 
Moll thermorelay. The effective slit width at 2.5u 
was 0.015. The percent absorption of water was 
first determined in order that a comparison could 
be made with the absorption curves of the solu- 


1E. K. Plyler and E. S. Barr, J. Chem. Phys. 2, 306 
(1934); E. K. Plyler and W. Gordy, J. Chem. Phys. 2, 
470 (1934). 


tions, the thickness of all cells being 0.1 mm. The 
absorption spectra for some of the solutions which 
were studied are shown in Fig. 1. It can be seen 
that the absorption of all the solutions in the 
region of 2.4u is greater than that of water. 
Aqueous solutions of KOH, NaOH, ZnBre, and 
ZnCl, of 5 and 10 N concentrations, 5 N Na2CQOs, 
18 and 36 percent HCl, and 17 and 34 percent 
HBr were studied. The ratio of the percent ab- 
sorption of the solutions to that of water was 
plotted and the results are shown in Fig. 2. 
The water content of the solution cells was less 
than in the water cell and the transmission is 
greater than 100 percent at some wave-lengths. 
Also, because of the shift of the 1.954 water band 
to longer wave-lengths in some of the solutions, a 
minimum appears at about 2y, but this is not an 
actual band. 
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The HCI and HBr solutions have bands at 
2.30u and 2.55u and the NaOH and KOH solu- 
tions at 2.30u and 2.45u. The hydrolyzing salts 
have bands at 2.30, 2.45y and 2.55. The absorp- 
tion spectra for the more concentrated solutions 
are shown in the second part of Fig. 2. The same 
bands are present with greater intensity. The 
absorption band at 2.30u is present in all the 
solutions and is probably due to hydration. A 
comparison of the intensity of the 2.454 band in 
the different hydrolyzing salts shows that the 
absorption is due to the OH ion. The nature of 
the absorption at 2.55u in the acids cannot be 
determined from the hydrolyzing salts. 

All of the solutions studied in this work have 
intense bands in the region of 1.85u. This band is 
sharp in the hydroxide and acid solutions, but 
much broader in the case of the hydrolyzing salts. 
These bands are due to combinations of the 
fundamental frequencies at longer wave-lengths 
and have been classified in a previous study." 
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Sodium Resonance Radiation and the Polymerization of Ethylene 


Joseru C. JUNGERS AND HuGu S. Taytor, Frick Chemical Laboratory, Princeton University 
(Received November 17, 1935) 


Ethylene has been found to quench the resonance radiation of sodium between 130 and 
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250°C. No measurable polymerization of ethylene results as a consequence of quenching, in 
contrast to previous results with excited mercury and cadmium. A catalytic polymerization 
of ethylene at sodium surfaces occurs above 250°C with an activation energy of 20 kcal. The 
inefficiency of the energy of sodium resonance radiation in producing polymerization is ascribed 


to unfavorable localization of the energy received by ethylene in the quenching process. 





T is known that the polymerization of ethylene 
can be accelerated by both excited mercury! 
and cadmium? by using the resonance radiation 
of these elements as exciting radiation. With 
mercury this occurs at room temperature; with 
cadmium around 200°C. Thermally, polymeriza- 
tion may be secured at temperatures in the 
neighborhood of 350—-400°C with an observed 
activation energy of 35-42 kcal.* 4 Quite recently 
the subject of ethylene polymerization has ac- 
quired a new interest by reason of the demon- 
stration that it is a chain reaction that can be 
brought about by the production in the ethylene 
system, at temperatures of 200—300°C, of free 
radicals; for example by thermal decomposition 
of metal alkyls* and azomethane,*® or by photo 
production of free radicals.” It was, therefore, of 
interest to test the effect of the resonance radi- 
ation of sodium as exciting agent, in an ethylene 
system containing sodium vapor. The choice of 
this source of energy derives its interest from the 
magnitude of the energy quantum involved 
(A=5897A; E=48 kcal.) since this energy is 
greater than the necessary thermal activation 
energy but too small of itself or even with the 
accompanying formation of sodium hydride 
(Onan =12.1 kcal.) to effect the breaking of a 
carbon-hydrogen linkage thus producing a free 
radical. 

We have, accordingly, studied the effect of 
ethylene on the quenching of sodium resonance 

1QOlson and Myers, J. Am. Chem. Soc. 48, 389 (1926); 
Bates and Taylor, ibid. 49, 2438 (1927). 

2 Bates and Taylor, J. Am. Chem. Soc. 50, 771 (1928). 

3 Pease, J. Am. Chem. Soc. 53, 613 (1931). 

4 Storch, J. Am. Chem. Soc. 56, 374 (1934). 

5 Taylor and Jones, J. Am. Chem. Soc. 52, 1111 (1930). 
1938) K. Rice and Sickman, J. Am. Chem. Soc. 57, 1374 


7Taylor and Emeleus, J. Am. Chem. Soc. 53, 562 
(1931). 
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radiation, the rate of polymerization of ethylene 
in presence and absence of such radiation and 
also a polymerization which occurs in somewhat 
higher ranges of temperature due to the presence 
of sodium surfaces or vapor in such reaction 
systems. 

EXPERIMENTAL 


Light source 

As source of resonance radiation we have used 
the commercial Zeiss sodium lamp of the stand- 
ard type used in polarimetric and other optical 
work. The lamp operates on a 110-volt a.c. 
circuit with suitable external resistance, at 1.3 
amperes. We have used such lamps for illumina- 
tion periods of many hours with good satis- 
faction. While we have made no quantitative 
measurements of light intensity, visual observa- 
tion showed a reasonable constancy of this light 
source. 


Reaction system 

We have used a Pyrex reaction system, the 
actual reaction vessel a cylinder 10 cm long by 2 
cm diameter, fitted with suitable manometric 
device and pump system for evacuation. Into 
the reaction vessel, after a standard high tem- 
perature exhaustion, the sodium metal was 
introduced by distillation and formed on the 
walls of the vessel a bright metallic mirror. 
Large amounts of the metal were always present. 


The temperature of the reaction system was 


controlled by means of suitable vapor baths in 
Pyrex tubes surrounding the reaction space. 
The ethylene used was obtained from cylinders, 
passed through a carbon dioxide-ether trap and 
purified before use by several fractional con- 
densations in liquid air with removal of per- 
manent gases by evacuation. In all cases the gas 
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was condensed onto bright metallic sodium 
before introduction into the reaction system, in 
order to minimize still further the possible oxygen 
concentration. 


Quenching of resonance radiation by ethylene 


Upon illumination of the heated reaction tube 
containing sodium and its vapor alone with the 
light from the resonance lamp the whole vessel 
became suffused with the characteristic yellow 
glow of the resonance radiation. This glow was 
especially striking in the temperature region of 
150—200°C. At the temperature of boiling di- 
phenyl (255°C) the concentration of sodium 
vapor is already sufficiently high that complete 
absorption occurs in the layers immediately 
adjacent to the lamp and the fluorescence is 
localized. Under conditions in which the glow is 
especially pronounced, the admission of even 5 
mm of ethylene sufficed to extinguish the flu- 
orescence as observed visually. Such experiments 
sufficed to indicate (a) the emission of resonance 
radiation from the lamp source and (b) the 
transfer of the energy to ethylene by collision. 


Photo-polymerization experiments 

The following Table I records the conditions 
under which the polymerization of ethylene 
under the influence of resonance radiation was 
attempted. In no case, between the reaction tem- 
perature limits of 132 and 255°C, was any 
acceleration of the polymerization by photo- 
sensitization observed. The results were uni- 
formly negative nor was any “‘permanent”’ gas, 
showing a pressure at liquid-air temperatures, 
produced during the several runs. 

Even at a pressure of 133 mm recorded in 
run 5 the only pressure diminution recorded can 
be accounted for by the thermal reaction dis- 
cussed in the succeeding section. An attempt to 
study the sensitized polymerization in presence 


TABLE I. Photosensitized polymerization of ethylene by 
sodium resonance radiation. 











TIME OF PRESSURE 

RUN Temp. PRESSURE ILLUMINATION CHANGE 
No. (°C) (mm at 25°C) (Hours) (mm) 

1 132 3.3 1.66 <0.1 

2 184 4.8 10 <0.1 

3 218 Le 5.5 <0.1 

4 255 5.25 5 <0.05 

5 255 133 20 ~!1 
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TABLE II. Thermal polymerization of ethylene in Pyrex 
vessels with and without a coating of sodium. 








PREs- 


Run Temp. SURE AP Al 





No. (°C) (mm) (mm) (hours) AP/At REMARKS 

1 302 182 20.6 44 (0.46 Sodium Present 
2 302 154.8 26.4 63 0.42 Sodium Present 
3 302 #167 2.4 38 0.06 Sodium Absent 
4 354 166 1.0 23 0.04 Sodium Absent 
5 354 180 22.0 17. 1.3 Sodium Present 
6 354 183 26.5 m «i Sodium Present 
7 302 184 23.8 70 0.34 Sodium Present 
8 354 182 30.3 21 1.4 Sodium Present 
9 255 181 23.0 232 0.1 Sodium Present 








of hydrogen or the sensitized hydrogenation of 
ethylene is prevented by the ready reaction of 
sodium and hydrogen, in which respect this agent 
differs from the previously studied cases of 
mercury and cadmium. 


Catalytic polymerization of ethylene 


It was thought that sodium vapor might be an 
efficient agent for the removal of traces of oxygen 
from ethylene gas and that, therefore, the 
thermal polymerization of ethylene might be 
reduced thereby, since, as is well known, traces 
of oxygen accelerate the polymerization process.*® 
A series of experiments on the rate of thermal 
polymerization carried out in reaction vessels 
coated with a bright metallic sodium deposit are 
recorded in Table II. 

Comparison of Experiments 3 and 4, in the 
absence of sodium, with the others in presence of 
sodium shows a marked accelerating influence on 
polymerization due to the presence of sodium. 
Experiment 9 indicates that the influence is 
measurable already at 255°C and accounts for 
the observed pressure change in the photo- 
sensitization experiment No. 5 of Table I under 
similar temperature and pressure conditions. 
That the influence of the sodium is catalytic is 
evident from the activation energy of the 
process. From Experiments 6, 7 and 8 we deduce 
E=20 kcal. The less accurate combination of 
Experiments 7 and 9 gives E=17 kcal. These 
values are much lower than the accepted value 
for thermal polymerization of 35-42 kcal.* 4 
This catalytic reaction makes difficult the study 
of photosensitization in the higher temperature 
ranges, but illumination with resonance radiation 
in this temperature interval produced no notice- 


able acceleration. 


8 Lenher, J. Am. Chem. Soc. 53, 3752 (1931). 
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DISCUSSION 


The absence of any polymerization of ethylene 
photosensitized by sodium diminishes markedly 
any expectation that examples of photosensitiza- 
tion using sodium will be found. Acetylene and 
hydrogen both react with sodium vapor and 
many other common reactants are similarly 
reactive towards sodium. It would seem, there- 
fore, that there is little prospect of extending the 
list of atom-photosensitized reactions to this 
element. 

The quenching of the sodium resonance radi- 
ation by ethylene indicates definitely that in the 
collisions with ethylene, at least 48 kcal. of 
energy per mole are transferred to the ethylene 
molecules. Although this is a larger energy than 
that required as activation energy in the thermal 
polymerization process no observable polymer- 
ization results. This is another example of the 
relative efficiency of thermal activation as com- 
pared with photoactivation. The negative result 
in this photosensitization confirms the view re- 
sulting from the studies with mercury and cad- 
mium as activating agents that the efficiency of 
these latter is to be associated with their ability 
to produce, in the reaction system, free radicals 
which are efficient agents in the initiation of a 
polymerization chain. With respect to the optical 
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S. CALDER 
and thermal methods of securing a particular 
reaction it should be observed that the method 
of calculation of activation energy in thermal 
processes takes account only of those collisions 
which, having the necessary activation energy, 
are successful, whereas, in the photosensitization 
studies, all the collisions are counted whether the 
energy is favorably located for reaction or not. 
Hence the lower apparent efficiency of light 
energy; hence, also, the observation that in 
thermal processes many collisions may occur 
with the necessary activation ernergy but un- 
favorable in the location of the activation. It is 
significant in this regard that Pease* computed 
from his data that the A factor in the velocity 
equation Ae-¥/2T was of the order of 0.001 of 
that to be expected of a straight bimolecular 
process and that the existence of a chain of any 
marked length would cause a still greater dis- 
crepancy between the actual and theoretical 
magnitudes of A. These considerations do not, 
however, exclude the possibility that the activa- 
tion energy deduced from thermal data may have 
a low value due to a reaction of ethylene with 
another substance of unknown nature, which 
reaction provides for the initiation of the chain. 
To this extent the thermal polymerization 
process is still obscure. 
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The Homogeneous Unimolecular Decomposition of Gaseous Alkyl Nitrites 


V. The Decomposition of Methyl Nitrite at Low Pressures 


E. W. R. STEACIE AND D. S. CALDER, Physical Chemistry Laboratory, McGill University, Montreal 
(Received November 19, 1935) 


The decomposition of methyl nitrite has been investigated from 0.005 to 5.0 cm in the 
temperature range 210—240°C. In this region the rate falls off steadily with diminishing pressure. 
The reaction is appreciably heterogeneous at the lower pressures. The results are in satisfactory 
agreement with the previous observations of Steacie and Shaw at higher pressures. The falling- 
off in rate with decreasing pressure is satisfactorily accounted for by the Rice-Ramsperger 
theory, assuming 13 degrees of freedom and a molecular diameter of 5.0 x 10-8 cm. 


INTRODUCTION 


N previous papers it has been shown that the 
decompositions of the gaseous alkyl nitrites 
constitute a series of simple first order reactions.! 


1 Steacie and Shaw, (I) Proc. Roy. Soc. A146, 388 (1934) ; 
(II) J. Chem. Phys. 2, 345 (1934); (III) ibid. 3, 344 
(1935); (IV) Proc. Roy. Soc. A151, 685 (1935). 





Recently Rice and Herzfeld? have suggested that 
the majority of first-order reactions may really 
be chain processes involving free radicals, and in 
a number of cases their suggestion has received 
striking experimental confirmation. A_ large 


2 Rice and Herzfeld, J. Am. Chem. Soc. 56, 284 (1934). 
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number of first-order reactions are therefore 
under suspicion. In the case of the nitrite de- 
compositions, however, the work of Rice and 
Radowskas? on the free radical decomposition of 
ethyl nitrite makes it almost certain that no 
complicating chain process occurs. These reac- 
tions remain, therefore, among the few that are 
unequivocally unimolecular. It thus seems ad- 
visable to obtain more comprehensive data on 


these reactions, and the present communication ° 


deals with the decomposition of methyl nitrite at 
low pressures. 


EXPERIMENTAL 


The progress of the reaction was followed by 
observing the rate of change of pressure in a 
constant volume system. A mercury manometer 
and a set of McLeod gauges were used to cover 
the various parts of the pressure range investi- 
gated. Two Pyrex reaction bulbs were used. One 
of these had a volume of 336 cc and a surface- 
volume ratio of about 0.30 cm, the other was 
packed with short lengths of 3-mm Pyrex tubing 
and had a surface-volume ratio of 4.5 cm~!. The 
reaction bulb was immersed in a well-stirred oil 
bath, the temperature of which could be hand 
regulated to 0.1°C. The reaction bulb was con- 
nected by means of the shortest possible length 
of 1/4-inch tubing to the McLeod gauge or to 
the manometer. The connections to the pumping 
system and methyl nitrite supply were made 
through mercury seals, no stopcock grease being 
used in or near the reaction system. 

Methyl nitrite was prepared and purified as 
previously described. It was stored as a liquid in 
a trap cooled with solid carbon dioxide. The de- 
sired amount for an experiment was obtained by 
filling a short length of capillary tubing with the 
gas to a predetermined pressure. 

As previously mentioned, mercury seals were 
used in or near the reaction system. In other 
parts of the apparatus stopcocks were used, and 
these were lubricated with Apiezon L grease on 
account of its negligible vapor pressure. All con- 
nections were made of wide tubing to facilitate 
evacuation. Between runs the system was evacu- 
ated to at least 10-° mm. If left standing after 
evacuation, the pressure at the end of 24 hours 
was always below 10-‘ mm with the empty reac- 





* Rice and Radowskas, J. Am. Chem. Soc. 57, 350 
(1935). 


tion vessel, and hence no trouble was experienced 
because of slow distillation from the walls of the 
vessel as reported by Rice and Sickman*‘ in their 
investigation of the decomposition of ethyl ether 
at low pressures. In the packed vessel the effect 
was noticed to a very slight extent. 

In carrying out an experiment, methyl] nitrite 
was introduced by quickly lowering and raising 
the mercury seal separating the supply from the 
reaction system. After 15 seconds the mercury in 
the McLeod gauge was raised and the initial 
pressure read. Blank experiments made with the 
reaction vessel cold showed that this allowed 
ample time for the establishment of pressure 
equilibrium. After reading the initial pressure, 
the mercury in the gauge was lowered to a fixed 
point so located that the gas in the gauge was still 
cut off from the reaction vessel. After the lapse 
of a definite time the mercury was lowered so as 
to connect the gauge to the reaction vessel, and 
then raised again and the new pressure deter- 
mined. By the use of this method only one point 
on the amount reacted-time curve was obtained 
from each experiment. It had the advantage, 
however, of permitting the use of gauges of large 
capacity, and at the same time keeping the 
dead space down to a few percent of the volume 
of the reaction system. 


CALCULATION OF RATES 


From the initial and final pressure we may 
calculate the percent pressure increase which has 
occurred during the reaction time. This must be 
corrected for the fact that not all the gas in the 
system was in the reaction vessel. The correction 
for the volume of gas contained in the gauge is 
simple and certain, since this gas has been isolated 
during the entire experiment. The correction for 
the gas in the connecting tubing between the 
vessel and the gauge is less definite, however, 
since it is uncertain what assumption should be 
made regarding diffusion. Since the tubing was 
short and of 1/4-inch diameter, we have as- 
sumed that complete mixing occurred. In any 
case the error from this cause is of the same order 
of magnitude as the experimental error. 

From the corrected pressure increase the “‘per- 
cent reaction” is calculated, assuming that 87 
percent pressure increase corresponds to com- 


4 Rice and Sickman, J. Am. Chem. Soc. 56, 1444 (1934). 
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plete reaction. The evidence for this will be 
given later. We thus have the percent reaction 
corresponding to a given time. From this the 
value of T4 (the time for 40 percent reaction) is 
calculated assuming that the reaction is first 
order. This is not quite justified, since we are 
dealing with a pressure region in which the reac- 
tion is tending to a higher order. No appreciable 
error is introduced, however, for the following 
reasons: (a) The times used were chosen so that 
the percent reaction was almost always in the 
range from 25 to 60. Theoretically the ratio 
T/T; for a first-order reaction is 1.77 and for a 
second-order reaction 2.00. Hence even if the 
order were as high as 1.5, the error involved in the 
method of calculation would not be greater than 
5 percent. (b) By trial it was found that in any 
single run the order did not deviate greatly from 
the first. (c) The results show no systematic 
variation with the value of the “percent re- 
action.”’ The rate of reaction has therefore been 
reported throughout as 104/T sec.—!. 


EXPERIMENTAL RESULTS 


In all about 150 runs were made at four differ- 
ent temperatures. Table I gives a summary of 
the data for representative runs in the empty 
reaction vessel. The data for all runs are given in 
Fig. 1 in the form of a log P—log 104/T graph. 

The results for the packed reaction vessel at 
494.3°K are given in Table II. Similar results 
were obtained at other temperatures, but there 
seems to be no object in giving them in detail. 

Some typical values of the total pressure in- 
crease at completion are given in Table III. The 
mean value for all runs was 87.0 percent, and 
this value has been used throughout in calculat- 
ing the amount of reaction from the observed 
pressure increase. As in the investigation of 
Steacie and Shaw, it was found that after the 
pressure in any run had reached a maximum, 
there was a very slow decrease. The data in 
Table III represent the maximum pressure 
reached. 


DISCUSSION 


The mean value of the pressure increase at 
completion is about 4.5 percent higher than that 
found by Steacie and Shaw. They found also 
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TABLE I. Data for representative experiments. 








PER- PER- 

INITIAL CENT INITIAL CENT 
PRESSURE TIME REAC- Ts PRESSURE TIME ReEac- T 4 
(cm) (sec.) TION (sec.) . (cm) (sec.) TION (sec.) 





Temperature—484.3°K 


0.00533 7200 20.5 16000 0.131 7200 47.8 5640 
0.0125 7200 20.0 16500 0.169 2700 26.6 6110 
0.0257 7200 23.1 14000 0.468 2700 37.5 4010 
0.0403 7200 28.2 11100 0.834 2700 43.6 3300 


0.0537 7200 32.6 9310 1.43 2700 49.2 2790 
0.0727 5400 29.3 7950 2.46 2700 53.3 2480 
0.0943 7200 41.8 6790 5.79 manometer 1540 





Temperature—494.3°K 
0.00479 5400 15.3 16600 0.129 5400 62.1 2840 
0.0107 5400 19.8 12500 0.268 1800 38.6 1880 
0.0174 5400 25.0 9550 0.750 1800 54.1 1180 


1500 55.6 943 
manometer 1080 
manometer 590 


0.0381 5400 38.5 5670 1.67 
0.0604 5400 48.5 4150 1.87 
0.0748 5400 51.5 3800 6.00 





Temperature—504.4°K 


0.00687 3600 20.3 8100 0.288 1200 51.4 850 
0.0120 3600 26.9 5860 0.515 1200 61.8 636 
0.0246 3600 37.4 3920 0.766 1200 67.4 547 
0.0437 3600 50.7 2600 1.59 1200 72.9 469 
0.0615 3600 56.5 2210 1.16 manometer 540 
0.0806 1200 32.9 1540 3.99 manometer 290 
0.104 3600 69.9 1530 6.06 manometer 280 





Temperature—514.4°K 


0.00455 1500 11.3 6390 0.250 600 48.1 467 
0.00999 1500 18.6 3720 0.493 600 61.1 324 
0.0164 1500 23.9 2810 0.786 600 66.5 280 


600 71.0 247 
manometer 160 
manometer 135 


0.0310 1500 34.7 1820 1.14 
0.0550 1500 39.3 1530 3.70 
0.0840 1500 56.1 931 6.18 
0.134 1500 64.0 751 








that after rising to its maximum value the pres- 
sure decreased again slowly due to some sort of 
condensation. It seems probable that the higher 
value found here is due to diminished condensa- 
tion at low pressures. 

In Fig. 1 the reaction velocity data are plotted 
together with the high pressure rate calculated 
from Steacie and Shaw’s equation 


K=1.84 X 10%? e~ 86400 RT, 


It will be seen that their rates are definitely 
higher than those found here. The discrepancy is 
equivalent to a difference of from 1 to 1.5° in the 
absolute temperature scale, and is therefore no 
greater than that which is usually found in com- 
paring the results of different investigations. 

A comparison of the results of Tables II and III 
shows that the reaction is appreciably heterogene- 
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Fic. 1. Effect of pressure on rate. Curve (1) 514.4°K; 
(2) 504.4°K; (3) 494.3°K; (4) 484.3°K. Full circles—this 
investigation; open circles—Steacie and Shaw’s high 
pressure rates. 


ous at low pressures. Since the surface-volume 
ratio in the packed vessel was 15 times that in 
the empty vessel, we may calculate roughly that 
the reaction is about 5 percent heterogeneous at 
an initial pressure of 1.5 cm, about 15 percent at 
0.05 cm, and about 40 to 50 percent heterogene- 
ous at the lowest pressures employed. It follows 
that below about 0.05 cm the results have no 
great significance. Above this pressure an error 
will be introduced into the form of the log P—log 
104/T4 curve, in the sense that the falling off in 
rate with diminishing pressure should be some- 
what greater than that found. The general form 
of the curve, however, will not be altered suffi- 
ciently to affect seriously the subsequent calcu- 
lations. 

The average value of the energy of activation 
at pressures from 0.05 to 1.5 cm is found to be 
36,100 calories per gram molecule in excellent 
agreement with Steacie and Shaw’s value of 
36,400. At lower pressures the heat of activation 
decreases rapidly with diminishing pressures. 
Such a decrease is to be expected on theoretical 
grounds, but it is doubtful in this case to what 
extent it is due merely to the increasing hetero- 
geneity of the reaction. . 

It is to be noted that the curves in Fig. 1 show 
no sign of upward concavity as is the case with 
ethyl ether at low pressures.‘ There is therefore 


TABLE II. Experiments in packed reaction vessel at 494.3°K. 








INITIAL PRESSURE TIME 
(cm) (sec.) 


0.0224 1800 
0.0449 2700 
0.0745 3600 
0.0998 1800 
0.600 1080 
1.53 1080 
1.96 manometer 
5.26 manometer 


PERCENT Tw 
REACTION (sec.) 


71.5 1100 
64.6 1330 
81.8 1080 
61.1 973 
62.2 566 
65.3 521 

360 

370 











TABLE III. Pressure increase at completion. 








PERCENT 
PRESSURE 
INCREASE AT 
COMPLETION 


INITIAL 
PREsS- 
SURE 
(cm) 


TEm- 
PERA- 

TURE 

(°K) 


PERCENT 
PRESSURE 
INCREASE AT 
COMPLETION 


INITIAL 
PRES- 
SURE 
(cm) 





Empty Reaction Vessel 


82.2 494.3 3.12 
87.5 504.4 0.0599 
89.2 504.4 4.49 
83.2 514.4 0.00670 
88.7 514.4 0.0213 
83.4 514.4 6.41 


0.0175 
0.0603 
0.0782 
1.66 
2.70 
1.86 


87.7 
88.2 
85.6 
91.7 
87.6 


Packed Reaction Vessel 


85.0 504.4 0.0162 
79.8 


0.0532 
0.0758 








no indication of a second first-order region coming 
into existence at low pressures. This, however, 
would hardly be expected with as simple a mole- 
cule as methy] nitrite. 

The falling off in rate as the pressure diminishes 
can be satisfactorily accounted for on the basis 
of Theory II of Rice and Ramsperger® if we as- 
sume a molecular diameter of 5.0 10-* cm, and 
that 13 degrees of freedom are involved. In 
making the calculation the rates have been 
divided by 2, since on the basis of the Rice 
mechanism’ two molecules are decomposed for 
each elementary act. As with the azo compounds 
the number of degrees of freedom involved is 
approximately half the total number avial- 
able (30). 


5 Rice and Ramsperger, J. Am. Chem. Soc. 49, 1617 
(1927). 
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The free energies of the diatomic alkali vapors Ke, Naz 
and Liz are calculated from spectroscopic data, and the 
equilibrium constants for the dissociation into the mon- 
atomic form are computed. With the aid of these constants, 
equations are deduced from the observed vapor pressures 
obtained by various investigators, which give the partial 
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pressures of atoms and molecules in the saturated vapors 
for temperatures up to the normal boiling point. In the 
case of sodium and potassium, the entropy of the solid 
metal is computed from the vapor pressure equation, and 
is compared with that predicted by the third law of 
thermodynamics. 








N spite of the numerous investigations that 
have been carried out to determine the vapor 
pressures of the alkalis, there has still been some 
difficulty in reconciling the results that have been 
obtained by various methods, and a correspond- 
ing uncertainty as to the correct form of the 
vapor pressure equation. One reason for this is 
that molecule formation may be appreciable in 
the vapor phase, and although attempts have 
been made both experimentally and by means of 
the Gibson-Heitler equation to determine the 
degree of such association, there has still been a 
chance of appreciable error in allowing for this 
effect. Thanks primarily to the work of Loomis 
and his associates,! however, there are available 
data sufficient to permit a precise calculation of 
the equilibrium between atoms and molecules in 
the vapors, and the results of such calculations, 
together with their application to the vapor 
pressure data are recorded briefly here. 

The spectroscopic constants used in the calcu- 
lation of (F°—E,°)/T for Ke, Naz and Lie are 
listed in cm= in Table I; they are weighted 
averages computed for a nonisotopic gas from 
the table of Loomis! and the isotope abundance 
ratios of Brewer and Kueck,? by means of the 
familiar relations connecting the reduced mass of 











TABLE I. 

Ke Naz Lie 
We 92.49 159.2 352.0 
WeXe 0.352 0.73 2.5 
Bo 0.05593 0.1543 0.6776 
— DoX108 0.0827 0.586 10.22 

0.000218 0.00082 0.00728 

AE)® 11800 17515 26275 








1 Loomis, Phys. Rev. 38, 2153 (1931). See this paper for 
a bibliography for the alkali molecular spectra. 
2 Brewer and Kueck, Phys. Rev. 46, 894 (1934). 
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a molecule with the spectroscopic constants.’ 
Even for the highest temperatures considered 
here, only the normal !2 state of the molecule 
need be taken into account, since the excited !Z 
and ‘II states make only negligible contributions 
to the state sum. The table also gives the 
values'!:* of AE,° in calories per mole for the 
reaction M,=2M. 

The resulting (F°—E,°)/T are listed in Table 
II,5 and were computed by methods already 
described ;* because of the small values of the 
fundamental frequencies such calculations would 
be exceedingly laborious if integral approxima- 


: TABLE IT. 








—(F° —Eo°)/T 
T°K Ke Naz Lie 
298.1 51.093 46.685 39.320 
300 = =51.148 46.738 39.369 


400 53.645 49.163 41.639 
500 55.606 51.077 43.447 


log K =log (Pyg)?/(Py.) 
Ke az Lie 





—5.247 —9.165 —15.540 
—5.190 —9.081 —15.414 
— 2.962 —5.797 —10.500 
—1.617 —3.816 — 7.539 


600 57.224 52.659 44.953 —0.715 —2.490 — 5.558 
700 58.603 54.009 46.246 —0.067 —1.538 — 4.138 
800 59.807 55.188 47.379 +0.420 —0.823 — 3.070 
900 60.877 56.235 48.388 0.800 —0.264 — 2.237 
1000 61.841 57.178 49.298 1.105 +0.184 — 1.570 
1200 63.53 58.83 50.89 1.562 0.857 — 0.564 
1400 64.98 60.24 52.26 1.887 1.339 + 0.155 
1600 66.25 61.47 53.45 2.129 1.702 0.697 
1800 67.39 62.58 54.52 2.315 1.981 1.118 
2000 68.41 63.58 55.48 +2.464 +2.204 + 1.456 
S%o9g.1 59.706 55.024 47.081 















3 For a discussion of the validity of averaging the spectro- 
scopic constants for the various isotopic varieties of mole- 
cule present, and then treating the particular element as 
nonisotopic, see Giauque and Overstreet, J. Am. Chem. 
Soc. 54, 1731 (1932); see also Gordon and Barnes, J. Chem. 
Phys. 1, 692 (1933). 

4 Loomis and Nusbaum, Phys. Rev. 40, 380 (1932). 

5 As is now customary, none of the entries includes the 
nuclear spin contribution. The constants used were 
hc/k=1.4324, R=1.9869, additive constant for the trans- 
lational free energy = —7.267. (F®°—E o°)/T can be inter- 
polated using not more than second differences; in the 
case of the equilibrium constants, it is better to interpolate 
T-log K. 

6 Gordon and Barnes, J. Chem. Phys. 1, 297 (1933). 
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tions’ were not used to avoid summation over all 
values of the vibrational quantum number. It 
should be noted in passing that with these mole- 
cules, the use of the rigid rotator-harmonic oscil- 
lator approximation to the energy levels (such as 
is employed in the simple Gibson-Heitler equa- 
tion) would lead to very considerable errors. For 
example, in the case of potassium, it would pre- 
dict for 2000° and 1000° 67.69 and 61.56, re- 
spectively, as the values of — (F°— Ey°)/T instead 
of 68.41 and 61.84 as given in Table IT; similarly 
for sodium, the entries for 2000° and 1000° would 
be 63.07 and 56.98 with resulting errors of 0.51 
and 0.20 cal./deg., respectively. 

Values of (F°—E,°)/T for the atoms can be 
readily computed since for 7<2000° only the 
normal ?S state (of quantum weight 2) need be 
considered. From these and the values of AE° 
given in Table I, log K =log (P)?/(Pm,) for the 
three metals can be obtained, and is entered in the 
last three columns of Table II. As an example of 
the errors introduced by the use of the simple 
Gibson-Heitler equation for these vapors, it 
should be noted that for 1000° and 2000° in the 
case of sodium, the equation would give K = 1.69 
and K = 207, respectively, as compared with the 
exact values from Table II, viz. 1.53 and 160. 
Experimental determinations of the degree of 
association of the alkali vapors have been made 
by Lewis* using a magnetic deflection method, 
but his values of the equilibrium constants are all 
considerably larger than those recorded in 
Table II. 

Of the three alkalis, sodium has by far the most 
complete and concordant vapor pressure data; 
the work of Edmondson and Egerton,’ of Thiele’® 
and of Rodebush and his associates! *: 3 covers 
a range of temperature of nearly 700° and a range 


of pressure from 1/10,000 of a millimeter to 


* To ensure numerical accuracy in the integral approxi- 
mations, the calculation was carried out for 2000° and 
1000° for each gas with two different values of M, see 
reference 6, Eq. (12), p. 304. Note also the misprint in the 
expression for J2 of Eq. (14), p. 304; the coefficient of the 
last term should be (1/27), not (47) as printed. 

8 Lewis, Zeits. f. Physik 69, 804 (1931). 

— and Egerton, Proc. Roy. Soc. A113, 520 
( ; 

10 Thiele, Ann. d. Physik (5) 14, 937 (1932). 
“—— and de Vries, J. Am. Chem. Soc. 47, 2488 
a Rodebush and Walters, J. Am. Chem. Soc. 52, 2654 
( s 
i930 and Henry, J. Am. Chem. Soc. 52, 3159 
( > 


nearly 500 millimeters. After allowing for the 
partial pressure of the molecules by means of 
Table II, the “best’’ equation" to fit the vapor 
pressures of the atoms (in mm of mercury) is 
given by 


log pPwra= 11.3245 —26,087/4.575T 
—1.174logT. (1) 


When used in conjunction with Table II, this 
leads to vapor pressures which in general agree 
with the results of experiment over the whole 
range to better than +0.5 percent and with no 
drift in the deviations. 

From Eq. (1), AEo® for the reaction Na(/) 
= Na(g) is 26,087 cal.; the latent heat of fusion of 
sodium" is 600 cal., and hence from the heat 
capacity data for crystalline sodium of Simon and 
Zeidler, AE,® for the reaction Na(s)=Na(g) is 
26,040 cal. Similarly, from Eq. (1), the entropy 
of Na(g) and the data of Simon and Zeidler for 
the range 298°—371°, S%93.; for crystalline sodium 
is 12.15 cal./deg.—in close but possibly fortui- 
tous agreement with the ‘third law’”’ value 12.2 
cal./deg. 

Fiock and Rodebush'® have measured the 
vapor pressure of potassium for the range 679°— 
1033°K (pressures from 4.6 to 783 mm); mole- 
cule formation is less prominent than with so- 
dium, but is still appreciable. After allowing for 
the partial pressure of Ke by means of Table II, 
px (in mm) can be represented by the equation!’ 


log px =11.410—22,210/4.575T 
—1.275 log T. (2) 


Eq. (2), in conjunction with Table II, predicts 
vapor pressures for this range which agree with 
Fiock and Rodebush’s results within their estim- 
ated experimental error of +1 percent. The 
agreement between these data for relatively high 
temperatures and those for the lower range 
(340°-470°K), however, is not nearly so satis- 


factory as in the case of sodium. For example, 


14 The coefficient of the last term assumes with Thiele 
(reference 10) Rengade’s value for heat capacity of liquid 
sodium—7.30 cal./deg./gram-atom (Rengade, Comptes 
rendus 156, 1897 (1913)). 

( po aee and Zeidler, Zeits. f. physik. Chemie 123, 383 
1926). 

16 . and Rodebush, J. Am. Chem. Soc. 48, 2522 
(1926). 

17 The coefficient of the last term corresponds to a heat 
capacity for liquid potassium of 7.50 cal./deg./gram- 
atom; see reference 15. 
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for 373.2°K, Eq. (2) predicts Px(=P)=1.32 
<10-° mm, while the equations fitted to their 
experimental results by Edmondson and Eger- 
ton,® Mayer !* and Neumann and Vélker"™ give 
1.82, 3.01 and 1.89X10-, respectively. 

From Eq. (2), from the heat of fusion of po- 
tassium (575 cal.) and from the thermal data of 
Simon and Zeidler,“ AE,° for the reaction K(s) 
=K(g) is 22,230 cal.; while from Eq. (2), the 
entropy of gaseous monatomic potassium and the 
heat capacity data of Simon and Zeidler, S%2x95.1 
for crystalline potassium is 14.60 cal./deg.— 
definitely Jess than the third law value, 15.2 
cal./deg. The reason for this discrepancy is not 
immediately apparent, but it should be noted 
that the Neumann-Volker vapor pressure equa- 
tion leads to a value of S°29., in agreement with 
the third law value, but is in definite disagreement 
with Fiock and Rodebush’s results even when al- 


18 Mayer, Zeits. f. Physik 67, 240 (1931). 
19 Neumann and Volker, Zeits. f. physik. Chemie A161, 
33 (1932). 
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lowance has been made for the molecule forma- 
tion. 

In the case of lithium, Hartmann and 
Schneider®® have measured the vapor pressure in 
the range 1204°-1353°K (pressures of from 17 to 
94 mm); after correcting their observed pressures 
for pri, by means of Table II, the pressure of the 
lithium atoms can be represented (as a rule to 
better than one percent by the equation) 


log pri = 11.889 — 38,900/4.575T 
~1.174 log T. (3) 


The normal boiling points for the three alkalis 
and the partial pressures of atoms and molecules 
at the boiling point in the saturated vapors, 
predicted by Eqs. (1), (2) and (3) and Table II 


aaa K Na Li 
Tp, 1030° 1162° 1599° 
bu 716 mm 658 mm 648 mm 
bu, 44 mm 102 mm 112 mm 


20 Hartmann and Schneider, Zeits. f. anorg. Chemie 180, 


275 (1929). 
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The electron distribution of (NH«4)2C20,-H.0 projected on (001) was obtained by a Fourier 


analysis making use of quantitative intensity measurements of x-ray scattering from (hk0). 
This analysis combined with a determination of the atomic parameters perpendicular to (001) 
gave the complete structure. Separations of covalently linked atoms in the C,O«" ~ group are 
C-—C, 1.581+0.01A; C—O, 1.25+0.02A; and C—O, 1.23+0.02A with the angle O—C—O 
=129°+2°. In this compound the O—C—O parts of an oxalate group are in planes at an 


INTRODUCTION 
HE crystal structures of a number of 
oxalates have been determined by the usual 
methods of analysis of x-ray diffraction data. 
In H.C,0,-2H,0,! a H:2C.0, and B H.C.02 the 
oxalate group apparently has a center of sym- 
metry and thus the four oxygen atoms of the 
group must be coplanar but cannot be crystal- 
lographically equivalent. Atomic separations as 


1W. H. Zachariasen, Zeits. f. Krist. 89, 442 (1934). 
2S. B. Hendricks, Zeits. f. Krist. 91, 48 (1935). 


angle of 28° while in H:C,0,-2H2O, a and 8 H2C.0, the parts are coplanar. 





found for these compounds and for KHC,O,, 
RbHC.O,, K.C,0,-H.O, and RbeC.0,: H2O are 
C—C, 1.57 to 1.60A; C—O, 1.20 to 1.30A; and 
angle O—C—O, 125°-134°.8 A more accurate 
determination of these values and a knowledge 
of the possible equivalence of the oxygen atoms 
can be obtained from a Fourier analysis of the 
electron distribution on a given plane. 


3 Note X-rays in Theory and Experiment, by A. H. 
Compton and S. K. Allison (New York, 1935), p. 444 
et seq., for a logical development of the subject. 
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ELECTRON DISTRIBUTION IN (NH,4)2C20,4-H20 103 


W. L. Bragg‘ has shown that the electron 
density per unit area at any point (x, y) in a 
projection with a center of symmetry is given by 


le « 
a(x, y) =— Lond mF im cos 27(nx+mz), 
eer 


where u and m are the indices of planes in a zone 
whose axis is normal to the plane of projection, 
Fim is the structure amplitude of the scattered 
wave, A is the area of the plane of projection, 
and Z is the total number of electrons in the 
projected volume. To evaluate this series it is 
necessary to have values of Fy» for all planes 
of a zone and these can be derived from the 
measured intensities which are proportional to 
Fin?.® The signs of Fim, however, are only to be 
found from an approximate analysis of the 
structure by the usual methods. 

It is perhaps apparent that the plane selected 
for the projection should be one that gives the 
minimum overlapping of atoms and groups in 
the projection. Further, in the case of an oxalate 
group with a center of symmetry, the most 
desirable plane would be one that is parallel to 
the plane of the oxygen atoms. These conditions 
are best fulfilled by a projection of the structure 
of (NH4)eC:04-H2O on (001). This compound, 
as it will develop, is of particular interest in that 
the oxalate group has a twofold axis rather than 
a center of symmetry and there is the possibility 
of formation of hydrogen bonds with the am- 
monium ion and the water molecule. 


EXPERIMENTAL PROCEDURE AND APPROXIMATE 
ANALYSIS OF THE STRUCTURE PROJECTED 


Previous work® on ammonium oxalate mono- 
hydrate led to the space group P2,2; with 
2 (NH4)2C204-H20 ] in a unit of structure having 
a=8.06A, b=10.34A, c=3.82A. We have verified 
the space group assignment and have found 
a=8.04A, b=10.27A, and c=3.82A; the first 
two values being determined from high order 
reflections of the pinacoids and the last from a 
layer line photograph about the c axis. 

Intensities of reflections for (hkO) were meas- 


ured for Mo Ka radiation from a constant source 


‘W. L. Bragg, Proc. Roy. Soc. A123, 537 (1929). 
°W. L. Bragg, Zeits. f. Krist. 69, 118 (1928). 
1938 *¢ Wood, Proc. Univ. of Durham Phil. Soc. 7, 111 





with an argon-filled ionization chamber; the 
ionization current being measured by an FP54 
vacuum tube direct current amplifier. The 
sample of (NH4)2eC:0,- HO was a cylindrical one, 
0.1 cm in radius (7), with its axis parallel to 
(khO) and it was illuminated through a narrow 
slit. The linear absorption coefficient (u) meas- 
ured with monochromatic Mo Ka radiation is 
1.42, and thus yr=0.14. With this arrangement 
the maximum error that could be introduced 
into the intensity measurement due to varying 
absorption in the sample’ would be less than 
one percent. These measurements were reduced 
to absolute values by comparison with the 
integrated reflection from (400) for Mo Ka radi- 
ation on an extended face of NaCl, this being 
taken as 0.98 X10~‘ from the work of James and 
Firth. The square of the structure amplitude, 
Fin?, was found according to the following 
paces for a mosaic crystal: 


nm = integrated reflection = 2Qre-?+" 


non 2 ,t +cos? 20 
Q= (-<F wn) —————=2.1910-°6F,,,,2, 
mc? 2 sin 20 


1+cos? 20 
gin 20 


where 


N is the number of scattering units per unit 
volume, and the other terms have their usual 
significance. » was taken in the first analysis as 
the linear absorption coefficient since it was 
thought that a correction for extinction could 
best be based on a knowledge of the structure 
and the atomic scattering factors. In the end it 
appeared, however, that extinction must be 
quite small. 

The general positions of the space group P2,2; 
- xyz; 1/2+x, 1/2-y, Z, 

Zyz; 1/2—x, 1/2+y, Z, 

which lead to special positions on the twofold 
axes, 00Z, 1/21/2Z; and 1/20Z; 01/22Z. 
Because of interatomic distance requirements 
the oxalate groups must be at one of these sets 
and the H,O molecules at the other, say 1/2 0 Z; 
01/2 Z. The NH, C, O; and O: atoms are in 
sets of general positions, thus giving thirteen 
parameters defining the structure. 


7F. C. Blake, Rev. Mod. Phys. 5, 169 (1933). 
8 James and Firth, Proc. Roy. Soc. A117, 62 (1927). 
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TABLE I. Observed and calculated values of the structure amplitudes (F) for (hkO) of (NH4)2C204-H,0-Mo Ka. 

















h\k 0 1 2 3 4 5 6 7 8 9 10 
0 obs 25.2 27.0 11.2 11.8 11.8 
calc 18.5 30.6 8.0 6.7 10.0 
1 obs 23.5 15.5 13.4 28.9 21.6 17.4 $13 0.0 13.1 12.2 
calc 18.4 10.0 10.1 33.2 13.8 11.9 7.9 0.4 11.7 7.5 
2 obs 3.2 12.6 0.0 29.9 19.4 36.1 21.9 6.4 18.4 14.8 9.9 
calc 6.0 10.1 3.1 25.5 12.0 44.6 13.3 6.7 13.2 7.2 7.2 
3 obs 32.5 29.7 15.2 6.7 0.0 13.3 Rae 9.4 6.6 10.6 
calc 38.6 33.8 9.9 3.6 1.4 10.3 o1 6.9 7.0 6.0 
4 obs 29.4 8.3 25.2 25.1 10.9 14.5 0.0 12.3 10.6 (3.3) (9.4) 
calc 21.6 2.6 25.8 23.9 8.1 9.5 0.5 9.3 6.5 2.0 5.7 
5 obs 10.2 8.0 0.0 (3.8) (5.7) (3.0) 23.5 0.0 
calc 5.8 9.1 1.1 3.8 3.3 2.0 19.8 1.4 
6 obs aod 8.2 6.8 22.0 0.0 (2.8) 0.0 8.6 0.0 (7.3) 
calc 4.5 4.8 5.3 19.0 0.4 3.7 2.6 7.7 0.2 3.5 
7 obs 13.8 15.5 16.6 0.0 (6.6) (7.8) (7.2) (3.8) 
calc 9.8 10.0 12.0 1.5 4.2 2.6 3.4 2.5 
8 obs 8.2 0.0 (7.5) 7.8 
calc 7.8 1.6 21 5.8 
10 obs 16.5 0.0 (7.1) 
calc 8.2 1.0 4.4 
Parenthetical values are visually estimated from an equatorial zone (hkO) Weissenberg photograph. Bars indicate negative values. 


In deriving the approximate projection on 
(001) it was assumed that the oxalate group has 
a coplanar arrangement of atoms of the configu- 
ration previously found. It would then follow 
that this plane is parallel to (001) and this is 
supported by the cleavage parallel to that 
direction and c as the direction of vibration of 
the slow ray. The high polarizability of the 
NH,* ion makes impossible further use of the 
optical data in the structure analysis.® To find 
the projection on (001) it is necessary to find 
the parameters of the ammonium group and the 
orientation of the oxalate group, which requires 
two parameters in both the a and 0d directions. 
The approximately normal intensities for reflec- 
tions from (10.0.0) and (0.10.0) (note Table I) 
suggest that the majority of the atoms scatter in 
phase from these planes. This fact, the intensities 
from other orders of the pinacoids, and the ab- 
sence of (220) require the following approximate 
parameters: 


Xc, 0.095 ; Xo,5 0.21; X 0,5 0.11; XNA, 0.40; 
Yo, 0.03; Yo,; —0.05; Yo,; 0.145; Ynu,; 0.24. 


The above parameter values were used in 
calculating the signs of Fx), use being made of 
assumed scattering curves for the various atoms 
and groups as follows: NH,*t as in NH,Cl,'0 

®S. B. Hendricks, Zeits. f. Krist. (A) 91, 290 (1935). 


10 Wyckoff and Armstrong, Zeits. f. Krist. 72, 324 
(1929-30). 


O and C as in urea! and water below sin 6/A 
=0.20 similar to NH,* and above that value 
similar to O. A Fourier analysis was made using 
these signs and the resulting improved parame- 
ters were used in a final calculation of the signs 
and values of ,F(10) listed in Table I, as well as 
in the final Fourier analysis. It is to be empha- 
sized that the final analysis does not depend 
upon arbitrary assumptions used in deriving the 
approximate structure. 

The method suggested by Beevers and 
Lipson” was used for summing the double 
Fourier series and this greatly shortened the 
labor involved. Projected electron densities were 
calculated for values of x and y by steps of 0.02 
save in the region near atomic centers where the 
interval was reduced to 0.01. Final values on 
the open scale are listed in Table II and the 
resulting contour map is shown as Fig. 1. The 
final parameters were taken as the average 
values for the centers of the contour areas and 
they are listed in Table III. The variation in 
these as determined from various contours did 
not exceed 0.001. An arbitrary splitting of the 
projection into areas per atom led to the results 
shown in Table III. The projected number of 
electrons exceeds that required by a neutral 
water molecule and is less than that required 


11 Wyckoff, Zeits. f. Krist. 81, 108 (1932). The carbon 


data are combined with Bernal’s results from graphite. 
12 Beevers and Lipson, Phil. Mag. 17, 855 (1934). 
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TABLE II. Calculated values of the projected electron density (¢) on (001) for (NH4)2C204- 1,0. 


ELECTRON DISTRIBUTION 









IN H:O 








(NH4)2C20,- 






























x\y 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 
0.00 2.00 1.69 1.11 0.92 1.27 1.70 1.69 1.38 1.20 1.33 1.49 1.41 1.29 
0.02 2.11 2.00 1.36 0.94 1.06 1.44 1.61 1.33 0.96 0.97 1.18 1.16 1.00 
0.04 2.60 2.85 2.27 1.43 1.06 1.32 1.70 1.79 1.48 1.14 0.86 0.91 0.73 
0.06 3.51 4.37 3.89 2.52 1.59 1.79 2.79 3.35 2.86 2.07 1.13 0.79 0.67 
0.08 4.39 5.79 5.38 3.53 2.51 2.79 4.59 5.73 4.96 2.96 1.30 0.75 0.84 
0.10 4.56 6.07 5.62 3.64 2.42 3.63 6.27 7.86 6.70 3.80 1.38 0.73 1.17 
0.12 3.46 4.80 4.38 2.65 1.89 3.63 6.79 8.60 7.20 3.82 1.13 0.51 1.28 
0.14 2.22 2.84 2.55 1.45 1.20 3.04 6.03 7.63 6.30 3.11 0.52 0.07 0.98 
0.16 1.66 1.38 1.27 0.90 1.03 5 4.48 5.49 4.46 2.11 0.15 0.00 0.61 
0.18 1.97 0.93 0.96 1.05 1.40 2.06 2.87 3.24 2.62 1.27 0.22 0.00 0.43 
0.20 2.34 0.99 1.03 1.58 1.92 1.76 1.70 1.59 1.30 0.86 0.60 0.59 0.71 
0.22 2.05 0.86 1.00 1.56 1.80 1.48 1.10 0.88 0.86 0.96 1.07 1.17 1.11 
0.24 1.36 0.61 0.79 1.30 1.47 1.30 1.10 1.06 1.13 1.24 1.43 1.49 1.39 
0.26 0.84 0.68 0.93 1.22 1.36 1.41 1.53 1.71 1.69 1.55 1.45 1.40 1.33 
0.28 0.86 1.06 1.40 1.54 1.58 1.64 1.96 2.19 2.13 1.75 1.44 1.28 1.23 
0.30 1.01 1.35 1.75 1.78 1.62 1.71 1.95 2.16 2.11 1.74 1.60 1.62 1.52 
0.32 1.02 1.29 1.59 1.64 1.47 1.47 1.55 1.54 1.47 1.63 2.15 2.67 2.57 
0.34 0.94 1.03 1.21 1.12 1.07 1.16 1.19 1.05 1.10 1.74 3.10 4.26 4.28 
0.36 1.06 1.08 1.18 1.05 1.06 1.19 1.19 0.93 0.97 2.08 4.22 5.83 5.94 - 
0.38 1.40 1.39 1.32 1.19 1.12 1.42 1.45 1.19 1.32 2.64 4.94 6.81 6.69 
0.40 1.74 1.73 1.45 1.07 1.03 1.31 1.50 1.33 1.48 2.73 4.98 6.49 6.27 
0.42 2.20 2.12 1.43 0.68 0.51 0.81 1.13 1.13 1.21 2.30 4.07 5.15 4.83 
0.44 3.37 3.01 1.64 0.33 0.00 0.25 0.68 0.75 0.80 1.52 2.77 3.57 3.27 
0.46 5.49 4.74 2.60 0.57 —0.24 0.09 0.54 0.77 0.43 0.88 1.77 2.31 2.03 
0.48 7.76 6.64 3.89 1.29 0.24 0.36 0.68 0.58 0.33 0.67 1.19 1.56 1.57 
x\y 0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50 
0.00 1.25 1.30 1.11 0.75 0.54 0.64 0.75 0.59 0.64 1.83 4.49 7.45 8.75 
0.02 0.97 1.09 Liz 0.94 0.62 0.47 0.48 0.52 0.88 1.79 3.93 6.51 7.76 
0.04 0.70 0.84 1.13 1.07 0.72 0.34 0.12 0.26 0.47 1.17 2.55 4.50 5.49 
0.06 0.67 0.72 0.93 1.34 0.88 0.53 0.36 0.34 0.51 0.85 1.64 2.67 3.37 
0.08 0.97 0.88 0.82 0.96 1.10 1.02 0.86 0.81 0.86 1.31 1.27 1.77 2.20 
0.10 1.44 1.05 0.67 0.84 1.27 1.48 1.39 1.17 1.19 1.25 1.37 1.53 1.74 
0.12 1.65 1.28 0.59 0.58 1.26 1.56 1.40 1.09 1.07 Be 1.28 1.36 1.40 
0.14 1.48 1.13 0.63 0.71 1.21 1.37 1.01 0.71 0.83 1.21 1.33 1.16 1.06 
0.16 1.02 0.94 0.76 0.89 1.23 1.20 0.71 0.51 0.81 1.26 1.27 1.10 0.94 
0.18 0.77 0.76 0.82 1.20 1.36 1.13 0.69 0.60 1.00 1.42 1.40 1.13 1.02 
0.20 0.62 0.71 0.93 1.24 1.36 1.10 0.78 0.87 1.23 1.50 1.38 1.05 1.01 
0.22 0.97 0.83 0.88 1.06 1.02 1.03 0.90 0.96 1.53 1.77 1.53 1.02 0.86 
0.24 1.24 0.96 0.86 0.86 0.70 0.54 0.64 1.44 2.30 2.82 2.50 1.64 0.84 
0.26 1.22 1.03 0.95 0.85 0.59 0.37 0.79 2.06 3.81 4.95 4.54 2.98 1.36 
0.28 1.15 1.11 1.17 1.08 0.71 0.33 0.83 2.86 5.47 7.18 6.73 4.48 2.05 
0.30 1.26 1.15 1.29 1.38 0.99 0.49 0.96 3.13 6.23 8.22 7.67 5.10 2.34 
0.32 1.93 1.32 1.32 1.45 1.35 0.66 0.90 2.74 5.65 7.25 6.63 4.17 1.97 
0.34 3.07 1.75 1.25 1.40 1.35 0.87 1.26 1.44 3.84 4.82 4.23 2.72 1.66 
0.36 4.11 2.17 1.23 1.28 1.44 1.15 0.96 1.47 2.33 3.43 2.09 1.80 2.22 
0.38 4.71 2.39 1.27 1.36 1.51 1.52 1.18 1.29 1.58 1.40 1.17 1.84 3.46 
0.40 4.37 2.36 1.44 1.47 1.60 1.16 1.14 1.29 1.37 1.16 1.13 2.27 4.56 
0.42 3.42 2.03 1.58 1.66 1.54 1.07 0.93 1.22 1.39 1.28 1.39 2.45 4.39 
0.44 2.41 1.73 1.68 1.79 1.52 0.88 0.79 1.02 1.33 1.31 1.34 2.15 3.51 
0.46 1.75 1.57 1.90 2.00 1.33 1.20 1.05 1.35 1.33 1.21 1.22 1.79 2.60 
0.48 1.44 1.54 1.83 1.70 1.25 1.38 1.50 1.63 1.38 0.94 1.10 1.64 2.11 












for carbon, while the other values correspond to 
O,-, Oz. and NHy,*t. These particular values 
depend upon choice of area but it is to be noted 
that the value of o is greater for O, than for O2 
and that this would be more pronounced if the 





area of Os» were increased or that of O, decreased. 
The oxygen atoms thus apparently have neither 
equivalent surroundings nor similar electron 
distributions. 

Values of Fox) calculated using the above 
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Fic. 1. The approximate electron distribution of 


(NH4)2C204-H:O projected on (001), with a contour 
interval of 1 electron/A?, the lowest contour being 2. 
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parameters and the indicated atomic scattering 
curves are listed in Table I. The agreement is 
not perfect, a major part of the discrepancy 
undoubtedly being due to the use of the particu- 
lar scattering curves. All absent reflections, 
however, are in close agreement and these are 
particularly sensitive for testing the possible 
correctness of a structure. Agreement between 
calculated and observed values of Fiaxo) for 
strong reflections indicates the absence of ex- 
tinction. 

In this projection the separation of the atoms 
in a C,0;,-~ group are C—C, 1.581A; O,—C, 
1.216A; and O.—C, 1.200A with the angle 
O—C—O=124°. Since the oxalate group has a 
twofold axis the projected separation of the 
carbon atoms is their true separation while the 
O-—C distances and the angle O—C—O are 
minimum values. 


COMPLETE STRUCTURE ANALYSIS 
AND DISCUSSION 


In order to complete the structure analysis it 
is necessary to determine the parameters along 
the c axis, the normal to the plane of the Fourier 
analysis. Trial and error methods can safely be 
used at least in determining the geometry of the 
oxalate group. The observed intensities (004)a., 
(014)m.w., (034)m. cannot be explained if the 
atoms of the C204 group are coplanar. Systematic 


TABLE III. Parameter values and electron counts for 
(NH 4)2C204- H20. 











ATOM PARAMETERS ELECTRON COUNT 
Area, Av.a@ No.of 
x y Z A? (001) Electrons 
Cc 0.092 0.027 0.07 2.25 2.218 5.0 
Or 0.200 —0.056 0.14 5.00 1.870 9.3 
Or 0.118 0.142 0.00 4.58 1.763 8.1 
NH, 0.386 0.228 0.43 5.70 1.763 10.1 
HO —0.20 6.32 1.740 11.0 








consideration of the intensities of (00/) and (Ok/) 
leads to the parameters listed in Table III. 
Agreement between intensities estimated visually 
from a Weissenberg photograph and calculated 
values of Fox1)/fo are listed in Table IV; the 
assumption being made that fun, =fo =fn,0 = 2fc, 
where f; is the scattering power. These parameter 
values are probably accurate to 0.01 for the 
atoms of the C,0,-~ group and 0.03 for NH,* 
and H.O. Distances within the oxalate group 
are thus found to be: C—C=1.581+0.01A; 
O,—C=1.25+40.02A; and O.—C=1.23+0.02A 
with the angle, O-C—O=129°+2°. Measure- 
ments on other oxalates mentioned in the 
introduction agree with these values within 
limits of experimental error; this being particu- 
larly true for 'H2C,0,-2H,O as analyzed by 
Zachariasen who found C—C=1.59+0.07A; 
O—C=1.25+0.05A and the angle O—C—O 
=126°+3°. The angle between the O—C—O 
plane and (001) is 14° and thus the angle between 
the two O—C—O planes of an oxalate group is 
28°, this being diagrammatically shown in Fig. 2. 

Arrangements of atoms surrounding . NH,* 
ions and water molecules are diagrammatically 
shown in Fig. 2. Water molecules are tetra- 


TABLE IV. Estimated intensities of reflection for (Okl), 
Cu Ke radiation, of (NH4)2C204- H20 and calculated 
values of F/fo. 











k/1 1 2 3 4 

0 1.4, s. La Ge D. 0.4, a. 0.0, a. 

1 0.5, m. 1.0, m.-m.s. 1.3, m.-m.s. 0.7, w. 
2 = 1.8, v.s. 1.0, m 1.0, m.w. 0.5, v.w. 
3. 0.0, a. 0.9, m. 0.5, w. 1.7, m. 
4 10,ms. 0.7,w 1.5,m.-m.s. 1.0; m.w.-w. 
5 02, a. 0.7, w. Li, m. 0.2, a. 

6. i3,ms. 6.1, vw. 0.6, a. 

7 13,m 0.3, a. 

8 0.2, a 0.7, v.w. 

9 0.8, v.w. 1.8, m.s. 
10 1.1, Ww 1.4, m.w. 
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Fic. 2. The configuration of the oxalate group and the atomic surroundings of the water 
molecule and the ammonium ion in (NH,4)sC,0,-H,0. 


hedrally surrounded by two NH,* ions and two 
oxygen atoms, (O:), having the higher value of co. 
The minimum distances of approach NH,—H,O 
and H.O—OQ, considerably exceed the value of 
2.55A required for hydrogen bond formation, 
but the tetrahedral arrangement is in agreement 
with Bernal’s and Fowler’s'® model of the water 
molecule. Ammonium ions too are tetrahedrally 
surrounded by two sets of oxygen atoms as 
shown in the figure. Interatomic distances for one 
set are near 2.80A and for the other near 3.20A. 
The observation that the planes of the 
O-—C-—O parts of the C,0O,-~ group are at an 
angle of 28° to one another is particularly 
interesting since the group is a plane one in 
H:C,0,-2H,O and a and 8 HeC.0,. Indeed it is 
the first instance yet encountered in crystal 
structure analysis where a group has different 
configurations in various salts. However, such 
a result is neither improbable nor unexpected 
due to the possibility of rotation about the C—C 
single bond. The effect is similar to that which 
leads to the optical activity of diphenyl deriva- 
tives'* as brought about by the steric effects of 
appropriate groups. 
~ 48 Bernal and Fowler, J. Chem. Phys. 1, 515 (1933). 


14 See for instance R. Adams and H. C. Yuan, Chem. 
Rev. 12, 261 (1933). 





The separation of the carbon atoms in the 
C,O,-~ group, namely 1.58+0.01A, considerably 
exceeds the value 1.54A found for the carbon to 
carbon distance in other aliphatic compounds 
that have been studied. This is possibly a result 
of repulsion between the negatively charged 
(—CO,-) parts of the group. 

Several comments should be made relative to 
previous work on oxalates. Structure analyses 
of K.C.0,;-H.0, RbeC20,-H.0, KHC,O, and 
RbHC;0; were based on the assumption that 
the atoms of the C,0,-~ group were coplanar. 
Reconsideration of the data from the first two 
compounds shows that the CO, parts of a C20, 
group might be displaced by as much as 10° 
from a coplanar arrangement, the angle being 
0°+10°. This does not seem to be probable for 
the other two salts but it is perhaps the factor 
that prevented a successful analysis of the 
structure of NaHC,0,-H,0. Such a rotation of 
the CO, parts of the C,0,;-~ group in ammonium 
oxalate monohydrate partially accounts for the 
departure of its optical properties from values 
predicted on the basis of a coplanar C,0,-~ 
group. 

We are indebted to Lola S. Deming for aid in 
summing the various Fourier series. 
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The rates of evaporation of n-dibutyl phthalate droplets 
ranging from 3u to 0.54 in radius were measured in a 
Hoag type of Millikan’s oil-drop apparatus under different 
conditions of vapor pressures. The data obtained are 
interpreted in two ways. One interpretation is based on an 
equation derived from Fourier’s equation of diffusion; the 
results thus obtained show that the surface tensions of 
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the droplets decrease with increasing curvature. The other 
is based on an equation derived on the basis of an assump- 
tion of a vapor film surrounding the evaporating droplet ; 
the results obtained here lead to the validation of Lord 
Kelvin’s equation. An explanation is also given for the 
abnormally high values of surface tension obtained by 
Woodland and Mack. 








HE capillary condensation hypothesis sug- 

gested by Zsigmondy’ in explanation of 
adsorption in porous bodies prompted a number 
of scientific workers* to study the mechanism of 
the adsorption process by the application of the 
Kelvin‘ equation 


In (P/Po) =+420V/RTr, (1) 


where P/P, is the relative vapor pressure of a 
liquid surface of radius of curvature r with 
respect to a plane surface, V is the molar volume 
of the liquid, o is the surface tension which is 
assumed to be independent of curvature,. and 
R and T have the usual significance. The sign 
of the right-hand term depends on whether the 
surface is concave or convex, the minus referring 
to liquids in capillaries and the plus to droplets. 

The results of this application were rather 
perplexing, since in many instances® the calcu- 
lated radii of the pores were less than the known 
diameters of the molecules of the adsorptive. 
On the other hand, Patrick® and co-workers 
extended the Zsigmondy hypothesis by assuming 
that the surface tension in small closed capillaries 
is increased because of a negative hydrostatic 


1 Part of thesis submitted in partial fulfillment for the 
degree of Master of Science in Chemistry at Howard 
-University. 

2 Zsigmondy, Zeits. f. anorg. Chemie 71, 356 (1911). 

3 Anderson, Zeits. f. physik. Chemie 88, 191 (1914); 
Coolidge, J. Am. Chem. Soc. 46, 596 (1924). 

4 Thomson, Phil. Mag. [4] 42, 448 (1871). 

5J. W. McBain, The Sorption of Gases and Vapors by 
Solids (George Routledge & Sons, London, 1932), p. 437. 

6 McGavack, Jr. and Patrick, J. Am. Chem. Soc. 42, 946 
(1920); Patrick and Greider, J. Phys. Chem. 29, 1031 
(1925); Patrick, Preston and Owens, J. Phys. Chem. 
29, 421 (1925); Patrick and Eberman, J. Phys. Chem. 29, 
220 (1925). 
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pressure in the surface, and were able to show 
the independence of the volume of adsorption 
from temperature and from the chemical nature 
of the adsorptive, thus showing the applicability 
of a modified Kelvin relationship. 

In this connection, Shereshefsky,’ studying 
the vapor pressure of water in small capillaries 
several w in diameter, found the lowerings con- 
siderably greater than the corresponding values 
obtained from the Kelvin equation. These results 
were interpreted to indicate an increase in the 
surface tension,, and therefore the need for 
extending this relationship. However, Goodris 
and Kulikova,® in studying the evaporation of 
water droplets of radii approximately equal to 
those of the capillaries, found the Kelvin equa- 
tion to hold. 

Recently Woodland and Mack,® applying the 
same method in a study of the rate of evaporation 
of droplets of dibutyl phthalate and dibutyl 
tartrate, report exceedingly high values for the 
surface tensions of these substances when calcu- 
lated with the aid of the Kelvin equation. These 
high values they seek to explain by assuming 
that the evaporating droplets are surrounded by 
a layer of vapor from 0.5 to 0.6u in thickness. 

In view of these discordant results together 
with the high values of surface tension and the 
unlikely thickness of the vapor layer reported by 
the latter workers, we undertook this study of 
the Kelvin relationship. 


7 Shereshefsky, J. Am. Chem. Soc. 50, 2966 (1928). 


8 Goodris and Kulikova, J. Russ. Phys.-Chem. Soc., 
Phys. Part 56, 167 (1924). 
— and Mack, Jr., J. Am. Chem. Soc. 55, 3149 










mes elhUr,rlC(iC ll -.—~ -. me - ae 


_—- ~~ |S ~~ A et 


oes © = OQ ft aaeéda da 





EVAPORATION OF SMALL DROPS 


APPARATUS AND METHOD 


Similar to Goodris and Kulikova* and Wood- 
land and Mack® we studied the rate of evapora- 
tion of tiny liquid droplets by utilizing the 
method developed by Millikan’® in the oil-drop 
experiment. It was also considered desirable first 
to obtain data that would lend themselves to 
direct comparison with the work of Woodland 
and Mack. For this reason the Hoag!! modifica- 
tion of the oil-drop apparatus and droplets of 
n-butyl phthalate” were used in this preliminary 
investigation. 

The set-up included an a.c. power supply 
unit!® to furnish the required plate potentials, 
which could be varied in a continuous manner 
over a range of 0-850 volts. The output of the 
unit was stabilized with a voltage regulator“ 
connected in series with the input. This arrange- 
ment reduced the fluctuations in the line to 


+ one percent. 

Four series of experiments were carried out in 
which droplets ranging from 3.0 to 0.54 in radius 
were studied. In one series droplets were allowed 
to evaporate in a condenser chamber saturated 
with dibutyl phthalate vapor. This condition 


was obtained by placing micro-boats containing 
the liquid in convenient positions in the con- 
denser. In the second series the evaporation of a 
number of droplets was allowed to proceed under 
conditions in which the condenser chamber con- 
tained activated charcoal. In a third series of 
experiments the charcoal was replaced by small 
sheets of blotting paper which were placed 
against the upper and lower plates. In this 
connection it may be pointed out that the blot- 
ting paper had no dielectric properties and 
merely affected the distance between the plates. 
But as the ratio e/d, as evident from Eqs. (2) 
and (3) given in the following section, is deter- 
mined by the measured values of V and m, an 
error in the value of d is of no consequence to 
the purposes of this investigation. In the latter 
two series the condenser was taken apart, 
cleaned, and fresh charcoal or blotting paper 


10 Millikan, Phys. Rev. 32, 389 (1911); Millikan, The 
Electron (University of Chicago Press, 1917), p. 168. 

1 W. M. Welch Co., Chicago. 

12 Eastman Kodak Co., Rochester. 

13 G—M Laboratories, Chicago. 

4 Delta Manufacturing Co., Cambridge, Mass. 
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introduced for each droplet studied. In the fourth 
series similar measurements were carried out 
under reduced pressures. This necessitated a 
more complicated apparatus which will be de- 
scribed in a forthcoming paper. 


RESULTS 


As it is hardly possible to include in this 
paper all the data and calculations obtained for 
the large number of droplets investigated, we 
are presenting complete data for only two 
phthalate droplets; these were studied in a 
chamber containing blotting paper. They were 
chosen because of the relatively large change in 
their radii during the period of observation, as 
well as because they are approximately of the 
same radius as the phthalate particle reported 
by Woodland and Mack. A further reason for 
selecting them was the fact that they show 
considerable differences in the rate of evaporation 
at comparable radii. 

The data for these droplets are contained in 
Tables I and II. The observed values are given 


TABLE I. Data of evaporation of droplet No. 20. Temperature 
21°C; No. of charges, 28; density, 1.051. 








1 2 3 4 5 6 
PLATE 
POTEN- TANGENT 
TIME ‘TIAL (volt/| Raprus 
(min.) (volts) sec.) i" 


2 221.5 0.127 1.673 
4 206.0 125 = 1.634 
176.4 A119 1.550 
148.0 114 1.463 
122.0 167 1372 
97.0 101 1.270 

73.5 093 1.156 

52.6 .084 1.036 

34.0 074 0.895 

18.3 062 728 


RATE OF 

EVAPOR. 
(g/sec. cm?) 
3.358 x 1078 
3.470 
3.677 
3.939 
4.218 
4.616 
5.149 
5.809 
6.828 
8.610 


RATE OF 
EVAPOR. 
(g/sec. cm) 


7.060 x 10™™ 
7.125 
7.161 
7.240 
7.262 
7.368 
7.474 
7.565 
7.684 
7.880 











TABLE II. Data of evaporation of droplet No. 22. Tempera- 
ture 21°C; No. of charges, 5; density, 1.051. 








1 2 3 4 5 6 
PLATE 
PoTEN- TANGENT 
TIME ‘TIAL (volt/| Raprus 
(min.) (volts) sec.) BM 


323 1.068 
283 .660 1.022 
245 .637 0.974 
172 568 .866 
107 .508 .739 
77 .462 662 
51 405 Ry a | 
28 334 473 


RATE OF 
EVAPor. 
(g/sec. cm?) 


7.63 X 10-8 
8.35 
8.83 

10.01 

12.32 

13.88 

16.07 

19.70 


RATE OF 
EVAPOR. 
(g/sec. cm) 


10.25 X10-" 
10.72 
10.83 
10.89 
11.43 
11.55 
11.65 
11.72 





0.660 
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Fig. 1. The evaporation of droplet No. 20. Curve I— 
The time variation of the surface. Curve I1I—The variation 
of the plate potential with time. 
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Fic. 2. The evaporation of droplet No. 22. Curve I— 
The time variation of the surface. Curve II—The variation 
of the plate potential with time. 


in the first two columns; the values of the third 
column were obtained with a tangent meter from 
the voltage-time curves given in Figs. 1 and 2 
drawn on a larger scale. The remaining columns 
were calculated from the first three by the 
application of the equations given below: 

1. Stokes’ law 
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v=2gr'p/9n, (2) 
where v is the rate of free fall, 7 is the radius of 
the droplet, p is the density of the liquid, and 
n is the viscosity of air. 

2. The relationship between the electrical field, 

the charge and the mass of the droplet 

Ve/d=mg = (4/3)rr* pg, (3) 
where V is the plate potential, d is the distance 
between the plates, e is the charge of the droplet, 
g is the acceleration constant, and m is the mass 
of the particle. In the fifth column are given 
the rates of evaporation expressed in grams per 
second per cm radius of droplet. In the sixth 
column are given the rates of evaporation in 
grams per second per sq. cm of surface of the 
droplet. 

It is of importance to point out that differ- 
ences in the rates of evaporation as appreciable 
as those shown in the above tables were obtained 
with the largest number of droplets studied in 
these experiments. Furthermore, these differ- 
ences were observed both when the evaporation 
was taking place in a condenser saturated with 
phthalate vapor and when the condenser con- 
tained an adsorbing material. Thus droplets 1 
and 2 in Table III evaporated presumably under 
equal environmental conditions, since the cham- 
ber in each instance was saturated with phthalate 
vapor; but the first droplet showed a rate of 
evaporation almost twice as large as the second, 
in spite of the fact that the former was a larger 
droplet. The same is also true of droplets 7 and 
11 which were evaporating in a chamber con- 
taining charcoal, and of droplets 18 and 21 
which were evaporating in a chamber containing 
sheets of blotting paper. 

It was also observed that the rate of evapora- 
ration was affected by the relative number of 
droplets entering the condenser. This number, 
not easily controllable, varied widely from experi- 
ment to experiment, as it depended on the hand 
pressure applied to the rubber bulb of the 
atomizer and the relative direction of the nozzle 
with respect to the perforations in the upper 
plate of the condenser. When the number of the 
droplets thus introduced into the chamber was 
relatively large, the rate of evaporation of the 
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TABLE III. Evaporation of droplets of n-butyl phthalate under various conditions. 











1 2 3 4 5 6 7 8 
LENGTH OF INITIAL FINAL RATE OF EVAPOR. CONDITION 
No. oF Oss. RADIUS Raprus (g/radius) OF o(calc.) 
DROPLET (min.) m “ INITIAL FINAL CONDENSER (ergs/sq. cm) 
1 26 1.33 1.28 0.643 X 107" 0.670 10-8 Saturated 6.2 x 108 
2 46 1.27 1.22 0.371 .388 Saturated 5.8 
3 32 1.10 1.00 .860 .946 Saturated 5.1 
4 30 0.92 0.85 549 .595 Saturated 4.2 
7 58 2.20 2.01 1.347 1.457 Charcoal 10.2 
11 32 1.25 1.18 0.719 0.761 Charcoal 5.7 
12 21 0.63 0.45 951 1.331 Charcoal 2.7 
18 30 2.66 1.80 17.120 18.230 Blotting paper 6.9-1.2 
21 10.5 i Be 0.89 8.09 10.43 Blotting paper 6.0-3.5 
Woodland 
& Mack 29.5 1.25 0.75 4.08 $.27 Charcoal 1.4-2.3 








droplet under observation was generally slower, 
and vice versa. 

It is therefore very likely that the apparent 
inconsistencies pointed out above are due to a 
variation in the degree of saturation of the 
chamber and of the adsorbing material. All the 
droplets, except the ones held for observation, 
whether charged or uncharged, fell to one or the 
other plate. While many of them in reaching the 
surface of the plates may have spread, a large 
fraction undoubtedly adhered in some modified 
droplet form, and caused a_ supersaturated 
atmosphere. It also explains why in the case of 
the blotting paper where all the droplets reaching 
the plates were absorbed the rate of evaporation 
was so much higher. 

However, all droplets, each under its own 
environmental conditions, behaved consistently 
alike as evidenced by the linear time variation 
of the surface shown in curve 1, Figs. 1 and 2, 
where the square of the radius is plotted against 
the time. 


TREATMENT OF RESULTS 


Maxwell’ pointed out that the rate of diffusion 
and the rate of heat conduction are equivalent 
processes and that mathematically analogous 
equations can be applied to both phenomena. 
Therefore, in a hollow sphere which is large in 
comparison with the droplet and which is lined 
with an absorbing medium, a droplet held in the 
center will evaporate as fast as the vapors are 
carried away by diffusion to the absorbing walls. 


® Maxwell, Scientific Papers (1890), Vol. II, p. 639. 





As the conditions of the present experiments 
approximate this ideal arrangement, the radial 
flow of vapor from the droplet to the walls of 
the condenser may be described by the Fourier 
equation 

dp/Ddt=d?(rp)/rdr’, (4) 


where D is the diffusion coefficient, p is the 
density of the vapor, and 7 is the radial distance 
from the center of the droplet. 

Assuming that for short time intervals the 
densities of the vapor in the different parts of 
the system are constant, the above equation 
becomes 

d?(rp)/rdr’*=0, (5) 


which on integration gives equation 
p =A +B/r, (6) 


where A and B are constants. 

From this general equation relating the density 
of the vapor with the radial position, after 
introducing appropriate boundary conditions, 
there is obtained!* 


—dm/dt= —4A4xr’Ddp/dr 
=4nr;RoD(pi—po)/(Ro—n), (7) 


where dm/dt is the rate of evaporation of the 
droplet, p: is the density of the vapor at the 
distance 7; equal to the radius of the droplet, 
and po is the density of the vapor at a distance 
Ry equal to the radius of the hollow sphere. 
When 7 is small in comparison with Ro, 


Eq. (7) can be simplified to 


16 Ingersoll and Zobel, Mathematical Theory of Heat 
Conduction (Ginn and Company, 1913), p. 27. 
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—dm/dt=4rDr\(pi— po), (8) 


which on assuming the vapor to obey the gas 
laws can be put into the form 


—dm/di=4nrr,MD(pi— po) /RT. (9) 


This equation, in which ; is the vapor pressure 
of the droplet, becomes identical with that 
derived by Langmuir,!’? when the vapor 
pressure of the absorbing medium is negligible 
with respect to ~:. Langmuir applied this equa- 
tion to explain the results of Morse'*® who found 
the rates of evaporation of small iodine spheres 
to be proportional to the radius. 

Assuming the droplet to be a perfect sphere, 
its mass m is given by 


m= (4/3) mrp. (10) 


Neglecting the variation of the density with the 
radius and differentiating, we obtain 


dm/dt=4rr pdr/dt = 3rpdS/dt, (11) 


where S is the area of the droplet. Combining 
with Eq. (9), we obtain 


—dS/dt=8xrMD(pi— po) /pRT. (12) 
This equation, first derived by Whytlaw-Gray 
and Whitaker'® to explain results obtained by 
them in a study of the vapor pressure of small 
water droplets, is in agreement with the constant 
rate of change of surface also found in our 
experiments, as shown in Figs. 1 and 2. 

In view of the appreciable effect of the vapor 
pressure of the absorbing medium upon the rate 
of evaporation of a droplet discussed in the 
preceding section, it is evident that pp in Eq. (9) 
cannot be assumed to be negligible, and that the 
rate of evaporation per unit radius, (1/7r;)dm/dt, 
is proportional to the difference (p:— 0) of the 
vapor pressures. This fact most likely explains 
why Woodland and Mack® in assuming the rate 
of evaporation per unit radius to be directly 
proportional to the vapor pressure of the droplet 
have obtained from the Kelvin equation ab- 
normally high values of surface tension; also 
the equally high values shown in column 8 of 


17 Langmuir, Phys. Rev. 12, 368 (1918). 

18 Morse, Proc. Am. Acad. 45, 362 (1910). 

19 Whytlaw-Gray and Whitaker, Proc. Leeds Phil. Soc., 
Sci. Part I (1926). 


SHERESHEFSKY AND S. 


STECKLER 


Table III and in column 6 of Tables IV and V. 
There is, however, a difference in that our 
values unmistakably decrease with decrease in 
droplet size, while the values of Woodland and 
Mack show a tendency to vary in the opposite 
direction. 

One method’ of testing the validity of the 
Kelvin equation is to calculate the absolute 
surface tension values corresponding to different 
droplet radii. As the application of this method 
to the available data is prevented by the lack of 
knowledge of fo in Eq. (9), it is possible to 
obtain evidence as to the constancy of o from 
calculated relative values. 

Putting Eq. (9) into the form 


Am = Ki(pi1— po), (13) 


where Am,=(1/r:)dm/dt and K,=4rMD/RT, 
we obtain for the relative rate of evaporation 
per unit radius for two droplets of radii 7; and 7 
the expression 


Am:/Am,= (p2— po) /(pi-— po). (14) 


Rearranging we get 


p2/ pi = Ame/Am,+ po/pi(1—Ame/Am,) (15) 


or in generalized form 


Pn, ‘pi = Amy ‘Am, + po, ‘pili <—_ Am,,/Am,), (16) 


where the relative vapor pressure p,/; of the 
droplets of radii 7, and 7 are expressed in 
terms of the measured quantities Am and the 
constant but unknown quantity po/ 1. 

Combining the Kelvin equation in the expo- 
nential form with Eq. (16) we obtain 


e2eV(L/r_—l/ 71) /RT 


=Am,/Am,+ po/pi(1—Am,/Am). (17) 

It is evident that by taking an arbitrary 
value of o, preferably the measured value, in 
Eq. (17) we should, if the surface tension is 
independent of curvature, obtain a constant 
value of o/; for all values of rp. 

The application of this equation to the avail- 
able data is shown in Tables IV and V. In 
column 3 of these tables are given the values of 
the left-hand member of the equation ; in column 
4 the values of Am,/Am;, and in column 7 the 
values of the ratio of the vapor pressure of the 
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TABLE IV. Surface tension and curvature. Droplet No. 20. 








4 5 6 
RELATIVE RELATIVE 
RATE OF RATE OF 
EVAPOR. 

(unit surface) 


1.0000 
1.0334 
1.0950 
1.1730 
1.2561 
1.3746 
1.5334 
1.7299 
2.0334 
2.5640 


1 2 3 
Pn/P1 
RELATIVE 
v.p. 
(Kelvin) 
1.00000 
1.00010 
1.00034 
1.00061 
1.00093 
1.00134 
1.00189 
1.00260 
1.00368 
1.00550 


RApDIus 
(Obs.) 
B 


¢ calc. 


EVAPOR. 
(ergs/sq.’cm) 


(unit radius) 


1.0000 
1.0092 
1.0143 
1.0255 
1.0286 
1.0436 
1.0584 
1.0715 
1.0084 
1.1162 


TIME @ cale. 
(min.) Po/ pr (ergs/sq. cm) Po/ pi 





1.673 
1.634 
1.550 
1.463 
1.372 
1.270 
1.156 
1.036 
0.895 
0.728 


0.9970 
.9964 
.9965 
-9964 
.9964 
.9965 
.9964 
.9965 
.9965 


0.9890 
.9766 
.9758 
.9675 
9693 
.9676 
.9636 
9584 
9527 


5 x 108 


. 
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TABLE V. Surface tension and curvature. Droplet No. 22. 








4 5 6 8 
RELATIVE RELATIVE 
RATE OF RATE OF 

EVAPOR. EVAPOR. 
(radius) (surface) 


1.000 1.000 
1.047 1.095 
1.057 1.147 
1.063 1.312 
1.113 1.614 1.18 
1.127 1.819 0.96 
1.137 2.107 74 
1.143 2.583 57 


1 2 3 
Py /Pi 
RApDIUS RELATIVE 
(obs.) V.p. o calc. 
Ud (ergs/sq. cm) 


@ calc. 


TIME 
(ergs/sq. cm) 


(min.) (Kelvin) Po/ tri Po/ tr 





| 


1.00000 
1.00064 
1.00154 
1.00295 
1.00406 
1.00565 
1.00766 
1.01621 


— 


1.068 
1.022 
0.974 
.866 
.739 
.662 
577 
473 


0.9969 
9957 
9951 
9951 
9951 
.9949 
9952 


0.995 
.988 
977 
974 
.969 
.960 


5.04 X 108 
2.80 
1.29 


Se ON wWh 
ADU 
SOK COAWNA 








absorbing medium and that of the droplet at 
its largest radius. 

Since the vapor pressure of the absorbing 
medium is determined by the initial conditions 
of each experiment and is undoubtedly constant 
during the whole period of observation, the 
decreasing ratio po/pi is a definite indication of 
a decreasing value of o with decreasing droplet 
size. The actual values of surface tension based on 
the assumption that the surface tension of the 
droplet at its largest radius was 32.6 ergs/sq. cm 
are given in column 8. They were obtained from 
Eq. (15) by using the first value of fo/p, and 
calculating for ¢. A similar tendency is also to 
be observed in the values of o in column 6 
obtained from the Kelvin equation on the basis 
of the erroneous assumption that p)=0, that is, 
that Am,/Am,= pn/ fr. 

To explain the abnormally high values of « 
they have obtained, Woodland and Mack postu- 
lated the presence of a layer of vapor surrounding 
the evaporating droplet. The thickness of this 
layer they calculate by applying Eq. (9) in 


which they neglect the term fo. Following this 
method of calculation we obtain values for the 
thicknesses of the vapor layer considerably 
larger than theirs. Thus, as shown in Table VI 


TABLE VI. The thickness of vapor layer. 








Drop.et No. 22 
Raprus THICKNESS OF 
(calc.) Vapor LAYER 
B ue 


4.814 3.746 
5.035 4.013 
5.087 4.113 
5.115 4.249 
5.367 4.628 
5.425 4.763 
5.472 4.895 
5.505 5.032 


Drop.et No. 20 
Rapius Raprus THICKNESS OF 
(obs.) (calc.) Vapor LAYER 
& B u 


1.673 3.316 1.643 
1.634 3.347 1.713 
1.550 3.363 1.813 
1.463 3.401 1.938 
1.372 3.411 2.039 
1.270 3.461 2.191 
1.156 3.510 2.354 
1.036 3.553 2.517 
0.895 3.609 2.714 
0.728 3.701 2.973 


RapDIus 
(obs.) 
7" 





1.068 
1.022 
974 
.866 
739 
662 
577 
A473 











the values for droplet 20 increase from one to 
four times the radius of the droplet, and for 
droplet 22 the thicknesses are even greater. It 
seems hardly plausible, in view of the mean free 
path under the present experimental conditions 
being about 0.01y, that the postulated saturated 
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layer should exceed so considerably this value. 
This is also borne out by the agreement between 
our observations and Eq. (12). This equation as 
shown above was derived from Eq. (9) in which 
r;, on the assumption of a vapor layer, stands for 
the radius of the evaporating vapor sphere, and 
from Eq. (11) in which 7 stands for the observed 
radius of the liquid droplet. From this it is seen 
that the agreement of the observed rate of 
change of surface with theory is possible either 
on the assumption that the ratio 7;/r is constant 
but larger than one, or that 7;=r. The first 
assumption must be ruled out, since as pointed 
out above the ratio 7,/r constantly increases 
with decreasing r; it therefore follows that 
within the experimental error the radii of the 
liquid droplet and of the vapor sphere coincide, 
and that the possible difference between 7 and 7 
does not exceed 0.02u. However, it is possible 
that the value of 3.110-* mm for # used in 
these calculations is too small, as indicated by 
the extrapolated value of 7.8 X10-* mm obtained 
from the vapor pressure measurements of Hick- 
man.?° On the other hand, it is obvious that for 
droplet 7 in Table III we would obtain negative 
values for the thickness of the vapor layer. 

Though the values for the thickness of the 
layer are too large to be accepted as evidence, 
there are other indications pointing to the 
possible correctness of the vapor-layer idea. 
One such indication is the constant rate of 
change of surface exhibited by the droplets in 
these experiments. As evident from Eq. (12) this 
behavior requires that (p:— po) be constant ; and, 
therefore, that ~,, the vapor pressure of the 
evaporating droplet, be constant during the 
whole period of observation. This condition 
would seem possible in the light of the Kelvin 
relationship if the droplet at a certain distance 
from it is surrounded by a vapor film of constant 
pressure. 

Let us assume that this film is within one 
mean free path from the surface of the droplet. 
In the region between the film and the surface 
of the droplet the number of molecules of vapor 
per cc is given by N=Np,/RT, where N is 
Avogadro’s number and 4; is the vapor pressure 
of the droplet. Of these molecules moving with 


2° Hickman, J. Phys. Chem. 34, 627 (1930). 
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the velocity u, 1/6 will be directed towards the 
film. The number of molecules that will pass in 
time dt an element of area dxdy in the plane of 
the film is 


N'=(1/6)dxdyudt Np,/RT. (18a) 


Similarly, the number of vapor molecules passing 
dxdy from a region of thickness udt adjacent to 
the other side of the film is 


N" =(1/6)dxdyudtN po/RT, (18b) 


where #p is the partial pressure of the vapor in 
this region. The net transfer of molecules through 
the element dxdy in the film is 


Ni =(1/6)dxdyudtN(p, — po) /RT, 
and the mass transferred is 
wi =(1/6)dxdyudtM(pi—po)/RT, (19) 


where M is the molecular weight of the vapor. 
This net transfer represents the mass which 
escapes the region between the droplet and the 
film and which diffuses to the walls of the 
vessel. The rate of evaporation per sq. cm is 


therefore 
> 


Am,= (2M/9xrRT) ‘(pi — po) ’ (20) 


since u=(8RT/2M)}. For simplicity and similar 
to Eq. (13) it can be written in the form 


Am, = Ko(pi-— po). (21) 


There is one important distinction between this 
equation and Eq. (13). In the latter Am, is the 
rate of evaporation per unit radius, while here 
it represents the rate per unit area. 

Relating the rate of evaporation of a droplet 
of radius 7, and of vapor pressure p, to that of a 
droplet of radius 7, and of vapor pressure fi, 
we obtained an expression equal in form to 
Eq. (16). Then combining it with the Kelvin 
equation we obtain an expression which is 
analogous to Eq. (17), and which equally allows 
the evaluation of the variation of surface tension 
with curvature. 

The results obtained in following the same 
method of calculations as in the case of the 
diffusion equation are presented in Tables I, II, 
IV and V. In the sixth column of Tables I and 
II observed rates of evaporation expressed in 
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grams per second per sq. cm are given. In the 
fifth column of Tables IV and V are given the 
relative rates of evaporation per sq. cm, and in 
the ninth column of the same tables the values 
of po/p:. It is to be observed that the values of 
po/p: are constant to a remarkable degree, the 
maximum deviation being in Table IV about 
0.06 percent and in Table V about 0.2 percent, 
and the deviation from the average being 0.04 
and 0.15 percent, respectively. If we were to 
exclude in each table the first value of this ratio 
corresponding to time intervals in which the 
steady state was probably not yet established, 
the deviations from the average would be reduced 
to about 0.01 percent in the case of both droplets. 
This fact seems therefore to indicate the validity 
of the Kelvin relationship for droplets ranging 
in size from 1.674 to 0.474; and since there is 
no reason to suspect any different behavior for 
droplets larger than 1.67y, the Kelvin equation 
seems to hold for all droplet sizes down to 0.5y, 
and most likely also for radii reasonably smaller. 

The contradictory results obtained from these 
two interpretations of the process of evaporation 
of liquid spheres are not at all surprising. If 
there is no discontinuity at the surface of the 
evaporating droplet, in the sense that the density 
gradient varies continuously from the surface of 
the liquid down to the absorbent, then the 
assumption of a vapor film or vapor layer around 
the droplet together with the conclusion based 
upon this assumption are erroneous. On the 
other hand, if the idea of a film of vapor or an 
equivalent form of discontinuity is correct, then 
Eq. (17) and the conclusions based upon it are 
not valid, since the Kelvin equation in the 
present form is not applicable to vapor spheres. 
In this connection it is worth pointing out that 
for a similar reason there is a serious objection 
to the use of the Kelvin equation for the calcu- 
lation of surface tensions of solids from evapora- 
tion measurements. 
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On theoretical grounds there is little reason 
to suspect that in microscopic droplets surface 
forces would extend beyond a distance of a few 
molecular layers. However, there are numerous 
investigations recorded in the literature in which 
evidence of surface forces acting through rela- 
tively great distances is presented. For instance, 
Trillat®! and Perrin” from x-ray measurements 
and Taylor®* from measurements of the index of 
refraction on liquid films find evidence pointing 
to a possible orientation extending through a 
relatively large number of molecular diameters. 
Derjaguin** reaches similar conclusions from 
measurements of the resistance to the movement 
of an oscillating lens presented by water surfaces. 

The distances through which the surface forces 
are reported to act are of the same order of 
magnitude as the diameters of the droplets in 
the present experiments. Unfortunately, it was 
not possible from our data to determine whether 
the range of action of the surface forces of the 
droplets was equally great, since we were only 
able to obtain relative values of the surface 
tensions. However, the decreasing surface tension 
with increasing curvature obtained on the basis 
of one of the mechanisms of evaporation dis- 
cussed above can only be explained?’ on the 
assumption that the surface forces of the droplets 
extend through distances equally great as those 
reported by other investigators. But to establish 
this it would be necessary to determine first 
which of the mechanisms represents the process 
of evaporation of tiny droplets. 


*1 Trillat, Comptes rendus 180, 1839 (1935). 

2 Perrin, Kolloid Zeits. 51, 2 (1930). 

*3 Taylor, J. Opt. Soc. Am. 23, 309 (1933). 

*4 Derjaguin, Zeits. f. Physik 84, 657 (1933). 

2 The effect of electrification on surface tension was 
shown by J. J. Thomson in Application of Dynamics to 
Physics and Chemistry (Macmillan and Company, 1888), 
p. 165, to be negligible. Nor is it likely that contamination 
of the surface is responsible for the effect observed, since 
the surface tensions of substances considered to be con- 
taminants and of the liquid are nearly alike, and high 
concentrations would have to be assumed to approach 
quantitatively anywhere near the lowerings observed. 
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The thermal decomposition of divinyl ether has been studied at temperatures between 











VOLUME 4 





460° and 500°C, and at pressures from 30 to 1000 mm. That the reaction has a chain mechanism, 
chains being broken on the walls at low pressure and in the gas phase at high pressure, is 
concluded from the effects of (a) increased surface, (b) added nitrogen (c) varying temperature 
coefficient with pressure and (d) induced decomposition with azomethane. The over-all order 
at high pressure is three-halves judged both from the individual velocity constants during a 


run and from the variation of velocity constant with pressure. 





HE thermal decomposition of numerous 
relatively complex organic compounds has 
been shown to proceed with an over-all first 
order rate. This has been satisfactorily accounted 
for by the free radical chain theory in many 
instances for which experimental data are avail- 
able; a treatment confirmed by a specific rate 
enhancement upon introduction of free radicals 
either thermally or photochemically. The theo- 
retical treatment of certain reactions however 
from this point of view is not without difficulty 
and rival theories must not be overlooked. The 
ethers constitute a series, the decompositions of 
which seem undoubtedly to follow a chain 
mechanism, and analysis of the chain steps has 
furnished indirect evaluations of the activation 
energies of various radical-molecule reactions. 
Divinyl ether being unsaturated might be ex- 
pected to furnish additional evidence for or 
against the free radical mechanism as a general 
scheme for ether decompositions. This choice 
would seem particularly interesting since from 
the point of view of a direct decomposition into 
product molecules the reaction would appear to 
be very simple. The prevalent view of ether 
decomposition prior to the chain mechanism 
suggestion, considered the possibility of a hydro- 
gen transfer across the oxygen in the ether 
molecule followed by rupture to yield a hydro- 
carbon and an unstable residue which subse- 
quently decomposed. Thus for methyl ether the 
hydrogen transfer would yield methane and 
formaldehyde. The parallel for divinyl ether 
would suggest the formation of ethylene and 
ketene, the latter in turn yielding ethylene and 
carbon monoxide. If the reaction were to proceed 
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at temperatures below which ethylene poly- 
merizes and decomposes, the products expected 
would be very simple and the reaction kinetics 
therefore easier of interpretation even than those 
of the saturated ether. Such a procedure, how- 
ever, is not found. On the contrary the com- 
plexity appears to be even greater. Evidence is 
presented that the reaction has a chain mechan- 
ism and from a theoretical point of view at least 
will be equally as complex as diethyl ether in the 
contingency of multiple interacting chains with 
possibly doubtful establishment of a steady state. 


> EXPERIMENTAL 


The divinyl ether employed was a Merck 
product obtained through the courtesy of Dr. 
Major to whom thanks are due. It was furnished 
containing one-hundredth of a percent of an 
oxidation inhibitor and for the first series of 
experiments was distilled from this. The boiling 
point was 28.3°C in agreement with that found 
by Miles and Menzies.'! Subsequent experiments 
were performed with the sample direct without 
further purification with no change in the 
decomposition rates. This is to be expected 
in view of the nonvolatility of the inhibitor 
present. It may be well to remark at this stage 
that polymerization of the liquid ether was 
observed on a few occasions but in all such 
cases oxygen may have been present. The liquid 
under its own vapor pressure alone showed no 
tendency to polymerize and the vapor main- 
tained for long periods of time at temperatures 
where decomposition was just too slow to 
measure showed no volume decrease. 


1 Miles and Menzies, J. Phys. Chem. 37, 425 (1933). 
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DECOMPOSITION OF DIVINYL 


The experimental procedure was similar to 
that previously reported for static studies from 
this laboratory. The reaction vessel of Pyrex 
glass was about 250 cc capacity. During an 
attempt to find a limiting rate at high pressure 
condensation of some product of the reaction in 
the capillaries outside the furnace necessitated 
an auxiliary heating of these to about 80°C. 
In the lower pressure range this was not found 
necessary. 

To investigate the products of decomposition 
several dynamic experiments were made. The 
ether vapor was forced through a small capillary 
under constant pressure, thence through a 
furnace, by varying the temperature of which 
the products corresponding .to more or less 
reaction as desired could be found. 


RESULTS AND DISCUSSION 


The results of a typical static experiment in 
the early stages are given in Table I. 


TABLE I. Initial Pressure 339 mm; temperature 500°C. 








PRESSURE 
CHANGE Runi (sec.-1) 


48 2.54 X 1078 
87 2.47 

119 2.40 

147 


TIME 


(min.) pj (sec.~!) 


8.11 X10-* 
8.49 

8.86 

9.41 

169 9.78 

187 wx 10.1 

201 ; 10.2 

216 10.8 

237 11.4 

252 11.9 
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The endpoint of the reaction appeared to be 
somewhat erratic varying between 1.05 and 1.24 
for the ratio of pressure change to initial pressure. 
No systematic variation with temperature was 
found though the higher the initial pressure the 
greater in general was the endpoint ratio. In 
Table I, are included the values of the uni- 
molecular and bimolecular velocity constants 
calculated in the usual way but assuming an 
endpoint ratio of unity; the time unit is the 
second, the concentration unit in k,; being the 
pressure directly observed in mm _ mercury. 
Since only the relative variation is required this 


is satisfactory. It will be observed that Runi_ 


decreases whilst ky; increases, the latter some- 
what more rapidly than the former. The over-all 
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order during a run is thus between one and two. 
For comparison of results at different pressures 
and temperatures use has therefore been made 
of the extrapolated first order constant obtained 
graphically from the pressure-time ratios in the 
early stages. This at least may be expected to 
give a more nearly true initial decomposition 
rate of the pure substance even when the 
products of reaction do affect the rate, than any 
other simple method. 

Table II contains these data for the various 
pressures and temperatures studied. Numerous 
other experiments were made at lower pressures 
down to about one mm. In all these cases, 
however, a marked induction period was ob- 


TABLE IT. 








INITIAL 
PRESSURE 
(mm) 


INITIAL 
PRESSURE 
(mm) 


kX 104 


(sec.~!) 





480°C 


Besseee 
QU NK Of ~10 = AO ur=~1 00 


SCAN OL HK aS 


HM NrIADSMUADAUNORE 
KH ABANBOUANADTO 


— 


470°C 
added surface 


& Co 


39 
89 
133 
288 
474 


1. 
1. 
ay 
2. 
2. 
2. 
2: 
3. 
3. 
4. 
4. 
4. 
5. 
5. 
5. 
6. 
2 

2. 
3 

3 

4. 
5. 
6 

7. 
7. 
8 

9. 
0. 
7. 
8. 
1. 
5. 
6. 
0. 
6 

8 

8 

1. 
6 

9 


QR W W DO DO DN tO 
CHM nnwe Pe Bin 











H. AUSTIN TAYLOR 








1 
Log P 20 25 30 


Fic. 1. Curves showing the variation of initial reaction rate 
with pressure at 460°, 470°, 480°, 490° and 500°C. 


servable, greater at lower temperatures. The 
data are not included even though the slope of 
the pressure-time curve after the induction 
period would seem to fall well on the log k—log p 
plot. Experience with several reactions with 
induction periods has shown that such agreement 
is frequently fortuitous and always misleading. 

The above data are plotted in the figure, the 
resulting curves being reminiscent of the type 
found in diethyl ether and propylamine® but 
not at all like those to be expected for a simple 
unimolecular reaction which would be concave 
downwards. There does seem to be evidence 
that the curves are becoming concave down- 
wards at the highest pressures measured though 
the rate of approach to a limit is too slow to 
justify any extrapolation. At the lowest pressures 
studied too, there does not appear to be any 
indication of a limiting value. The average 
slope of the curves is about 1/3 corresponding 
to a 4/3 order and at the highest pressures the 
slope is 1/2, that is the reaction is 3/2 order. 

The effect of added surface is shown in Fig. 1 
by the dotted curve. The reaction vessel was 
filled with short lengths of 5-mm tubing, the 
ends being fire-polished, whereby the calculated 
surface to volume ratio was increased 3.5 times. 
At the highest pressures the effect is very small 
but in the lowest pressure range the rate is 
decreased to an extent approximately equivalent 
to a decrease of ten degrees in temperature. 
This fact alone is suggestive of a chain mechan- 
ism for the reaction. 


20. K. Rice and D. V. Sickman, J. Am. Chem. Soc. 56, 
1444 (1934); 57, 22 (1935). 


Table III contains the data obtained in a 
study of the effect of various gases added to the 
system. The gas labeled products was simply 
the residue after some divinyl ether had com- 
pletely decomposed. An examination of these 
data shows that the effect of added gases is in 
general to reduce the rate of reaction. With the 
exception of hydrogen there is no reason to 
suspect from the results that specific chemical 
action was occurring in any case. With hydrogen, 
although the initial rate was the same as that for 
the ether alone the pressure-time curve reached 
its endpoint sooner yielding a value 0.9 as 
compared with 1.1 to 1.2 as previously men- 
tioned. It seems likely that this reduction of the 
endpoint ratio, is due to hydrogenation of un- 
saturated bodies shown later to be present. This 
does not absolutely preclude of course the possi- 
bility that chemical action might have occurred 
for example with the products even though the 
endpoint is apparently unchanged. Thus the 
ratio for 100 mm ether alone at 500°C was 1.12 
and for 112 mm ether in presence of 112 mm of 
product was also 1.12. The point is mentioned 
because the effect of products in reducing the 
initial rate of decomposition seems considerably 
greater in the lower ether pressure range than is 
found for example with the more probably inert 
gas nitrogen. Since ethylene and carbon mon- 
oxide are among the products their effects were 
tested separately but showed no such marked 
reduction of the rate as did the products. When 
a sufficiently large excess of ethylene was used a 
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relatively small rate reduction was observed. 
The results for nitrogen are particularly illumi- 
nating. For pressures of the ether around 100 to 
200 mm even twice that amount of nitrogen has 
a relatively small effect. At higher ether pressures 
a considerably smaller amount of nitrogen causes 
a marked reduction in the initial decomposition 
rate. Now it is precisely in this higher range of 
ether pressures that the effect of surface is 
least, hence that the reaction is a chain reaction, 
chains being broken on the walls at lower 
pressures and in the gas phase at higher pressures 
would seem a possible explanation. It would 
thus appear possible that the curves shown in 
the figure are complex, the one showing the 
variation with pressure of the over-all rate of 
reaction heterogeneously terminated, the other 
homogeneously. The apparent inflection should 
then correspond to a pressure at which the 
homogeneous and heterogeneous effects become 
approximately equal. It would be interesting to 
study the reaction in presence of very extensive 
surfaces at considerably higher pressure and also 
in large volume with a correspondingly small 
surface at lower pressure, to give an idea of the 
variation of the heterogeneous and homogeneous 
effects over a larger pressure range. The indica- 
tion herein is that the log k—log p plot is concave 
upwards for the heterogeneous effect and the 
reverse for the homogeneous. The former is 
explicable on the basis of increasing saturation 
of the surface for deactivation or chain breaking ; 
the latter since the over-all order is less than 
two whereas deactivating or chain-breaking 
collisions are presumably bimolecular. 

An examination of the figure shows that the 
log k—log p plots for the different temperatures 
are not parallel but approach each other in the 
low pressure range and diverge from each other 
at higher pressures. This is tantamount to saying 
that the energy of activation increases with 
increasing pressure. Thus it is found that around 
30 mm pressure the mean energy of activation, 
calculated from the curves for the different 
temperatures at constant pressure, is 47,000 
calories; at 100 mm it is 50,000 whilst at 300 
mm, 51,000 calories. Such a variation is outside 
the limits of experimental error, which should 
not be more than 1000 calories, and would be 
inexplicable for a simple reaction. In a chain of 
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reactions with the possibility of more than one 
chain carrier such a variation is at least expli- 
cable if not in fact to be expected. 

In an effort to determine the products of the 
reaction a small amount of the liquid divinyl 
ether was sealed in a small bomb and completely 
decomposed at 450°C. The resulting gases were 
then analyzed. The bomb after the decomposi- 
tion was dark brown in color but showed no 
signs of free carbon. A few drops of a colorless 
liquid were present which possessed an odor 
suggesting a paraffin hydrocarbon which volatil- 
ized around 70°C. A C—H balance from the gas 
analyses gave an end product with the formula 
CH;,.;. The data for the analyses are given in 
Table IV. Included in the table are the analyses 


TABLE IV. 
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made in the dynamic experiments by the method 
outlined above. The effluent gases from the 
furnace were bubbled through water to remove 
soluble products including unchanged ether, the 
gases in the above analyses being collected over 
water. The presence of ketene in the effluent gas 
was demonstrated by treatment with alkali and 
subsequent identification of acetic acid. The 
multiplicity of products renders a complete 
systematization of mechanism extremely diffi- 
cult. However even a cursory review of the 
possibilities of chain reactions in this decompo- 
sition will indicate the complexities to be ex- 
pected. 

There is very little evidence to indicate 
definitely what types of disruptions may be 
expected in a molecule as unsaturated as divinyl 
ether. The only apparent rule relative to un- 
saturated compounds’ is that with a doubly 
bound carbon atom, the alpha-bonds are stronger 
than normal whereas the beta-bonds are weaker. 


Actually in divinyl ether there are no beta-bonds 


3F. O. Rice and K. K. Rice, Aliphatic Free Radicals 
(Johns Hopkins Press, 1935). 




















































unless the C—O bonds may be considered as 
beta to the farther double bond. It would seem 
probable though that the ether linkages would 
be slightly weaker in this case than for the 
saturated ethers and hence the primary break 











CH, : CH—O—CH 





: CH,—-CH, : CHO+CH, : CH 
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would certainly occur there. Beyond that there 
is little precedent to indicate the relative 
strengths of the remaining bonds. The following 
equations summarize the majority of the more 
obvious possibilities. 





CH, : CH+CH2 : CH—O—CH : CH2—-CH, : CH2+CH2 : COCH : CH: 
CH, : CH+CH2 : CH—O—CH : CH2-CH2 : CH2+CH : CH—O—CH : CH 
CHz : COCH : CHe-CH2 : CO+CH: : CH 
CH : CH—O—CH : CH2—-CH : CH+CH:2 : CHO 
CH;+CH2 : CH—O—CH : CH:-CH,+CH: : COCH : CH2 
CH;+CH2 : CH—O—CH : CH2-CH,+CH : CH—O—CH : CH: 
2CH.2 ° CH-C,H, 
2CH;—C2Hs 


There is therefore a high probability of at 
least two distinct chain carriers with the possi- 
bility of two alternative modes of reaction of 
each, one of which regenerates itself whilst the 
other regenerates the second carrier. The two 
schemes are therefore considerably interlocked 
and judicious assumptions could easily account 
at least qualitatively for a variation of reaction 
order and a variation of energy of activation 
with pressure. Furthermore, as numerous efforts, 
which find no place here, have shown, judicious 
choice from among the above can give a three- 
halves order reaction rate which in most cases 
corresponds to an over-all activation energy 
around 50,000 calories as observed. It seems 
useless to be more specific at the present time 
until further evidence is forthcoming. 

It should be remarked in passing that the 
presence in quantity of acetylene and butadiene 
in the final products should not be expected 
since decomposition of the former sets in 
markedly around 500°C, of the latter below 
400°C, whilst polymerization of acetylene could 
account for the presence of condensable products 


CH3;+CH: : CH->CH,+CH : CH 




















observed in the high pressure runs. Moreover it 
should not be overlooked that addition reactions 
as well as polymerizations are possible in a 
system as unsaturated as that of vinyl ether and 
its products, though it is doubtful if they are 
likely to any appreciable extent under the 
temperature conditions used herein. 

The most hopeful indication of further evi- 
dence seems to be forthcoming from work now 
in progress in this laboratory by F. P. Jahn on 
the decomposition of divinyl ether induced 
thermally by azomethane. The reaction proceeds 
quite rapidly at 330°C with less than ten percent 
of azomethane. The over-all order is about 
three-halves though a resolution of this into the 
separate orders relative to the ether and azo- 
methane alone has not yet been made. It is not 
certain that the induced reaction is identical 
with the decomposition here studied; a lower 
endpoint signifies a difference. Such a difference 
may however actually afford a means of dis- 
tinguishing between the many possible chains 
suggested. 
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The potential energies for linear symmetrical configura- 
tions of the triatomic hydrogen molecule have been calcu- 
lated by the variational method with hydrogen-like 1s 
atomic orbitals by varying the “effective charge’’ and the 
amount of all polar and homopolar states in the eigen- 
function. The lowest energies (in kcal./mole = 1/23.06 ev) 
and corresponding separations between neighboring atoms 
(in Bohr radii) are given in the table. The activation 
energies are for the reaction H+H.—~H2+H, and are 


calculated with the same type of approximation for Hy, 
and Hs. All of the difficult three-center integrals for the 
repulsions between pairs of electrons were calculated 
explicitly in ellipsoidal coordinates without the use of 
infinite series. 

Experi- 


ment 
— 103. 


Wang plus 
Wang Polar 
—56.16 —67.09 
1.89 1.84 
30.68 25.15 ve 


H.L. plus 
Polar 
—60.39 


2.00 
13.63 


Heitler- 
London 
—53.11 

2.00 


19.07 


En; —3EH 


rH3 ss 
Act. Energy 





INTRODUCTION 


HE simplest triatomic systems are those 
formed from three hydrogen nuclei. The 
neutral triatomic molecule is of particular interest 
since it is the intermediate in reactions of the 
type 
D+H:—-DHH—DH+H (1) 
and also in the ortho-para hydrogen conversion. 
Various approximations for the energy have 
been considered for Hs. London! proposed to use 
hydrogen atomic orbitals and to neglect all but 
the diatomic Coulombic and exchange integrals. 
This led to values for H; very much too high as is 
known by comparison with experiment and cer- 
tain of the terms neglected are of the same order 
of magnitude as those which are retained. 
Coolidge and James? have calculated the energy 
for a single configuration of the linear symmetric 
H; molecule using hydrogen atomic orbitals to 
form the lowest homopolar state. Their particular 
configuration gives an energy more than a half 
volt higher than they would have obtained for a 
somewhat larger interatomic distance. A discus- 
sion of the criticism of James and Coolidge of 
the semi-empirical method of calculating activa- 
tion energies will be given in a forthcoming paper 
in the light of this work and the results of the 
method itself. 
In the present paper the variational method 
for the energy is carried through for all the linear 


1928) London, Zeits. f. Physik 46, 455 (1928); 50, 24 
24. S. Coolidge and H. M. James, J. Chem. Phys. 2, 
811 (1934). 


symmetric configurations using hydrogen-like 
atomic orbitals (of the 1s type) to form the four 
following eigenfunctions: 

(a) the best homopolar function. 

(b) the best linear combination of polar and homopolar 
functions. 

(c) the eigenfunction of (a) with the “effective charge”’ 
varied to give the lowest energy. ; 

(d) the eigenfunction of the type in (b) obtained by 
simultaneous variation of the screening constant and the 
polar-homopolar composition. 


II. THE ENERGY OF THE TRIATOMIC HYDROGEN 
MOLECULE 


All of the integrals used in the calculation of 
the energy of the triatomic hydrogen molecule 
are defined and discussed in the appendix. The 
Hamiltonian operator for the electronic energy is: 


H=Hat+Hs+Hce+Hpo, (2) 
where 
Ha = —(h?/82?m)[(VP+V2?+V3"] 
=kinetic energy of the electrons, 
Hp=e*[ras + Pac +r! ] 


=energy of mutual repulsions of nuclei, 


3 
He=-e? [rai +r +8] 


i=1 


=energy of attraction of electrons to the 
nuclei, 


Hp= erie! + rig +re3- | 


=energy of mutual repulsions of electrons. 
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Here a, b and ¢ are atoms arranged in the sym- 
metrical linear configuration with 6 the central 
atom. The distance between neighboring atoms 
is Rap (ay9=0.528 Angstrom, i.e., the radius of 
the first Bohr orbit for atomic hydrogen). The 
electrons are denoted by 1, 2, and 3. r;; is the 
distance between the 7th and jth particle. The 
calculated energy for the molecule is given as: 


E=fffV*t¥HWVdridtedr3/ fff V*Vdridredr3. (3) 


By the variational principle £ is the actual energy 
of the molecule when VW is really an eigenfunction 
of H, and E is greater than the lowest energy of 
the molecule when W is only an approximate 
eigenfunction. In this treatment we obtain the 
energy using various types of approximate eigen- 
functions built from the 1s hydrogenic atomic 
orbitals: 


a(1) =(2?/m)te-#ra, (1) = (23 / ar) be #701, 


(4) 
c(1) = (28/1) te-*ra1, 


where z is the “effective nuclear charge.’’ These 
eigenfunctions must be antisymmetric with re- 
spect to an interchange of any two electrons 
because of the Pauli exclusion principle. Without 
loss of generality they are taken as eigenfunc- 
tions for S? and S,, the square and the z com- 
ponent of the total spin, respectively; the lowest 
energy usually corresponds to the lowest allowed 
value of S*. The fact that H is symmetric with 
respect to an interchange of atoms a and c¢, 
allows the eigenfunctions to be separated into 
symmetric and antisymmetric functions with 
respect to this interchange, and the two types of 
functions are noncombining. In the London ap- 
proximation, i.e., where the multiple exchange 
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and overlap integrals are neglected, it is easy to 
show that the antisymmetric eigenfunctions cor- 
respond to the states of lowest energy. This is also 
found to be the case in the present calculation. 


1. The Heitler-London approximation 

The simplest form of eigenfunction which we 
can take for the triatomic hydrogen molecule 
corresponds to a pure homopolar state with z= 1. 
This eigenfunction is the three-atom analog of 
that used by Heitler and London?’ (H.L.) for the 
normal diatomic hydrogen molecule, and we 
should expect the energy obtained in this manner 
for the triatomic molecule to be of the same order 
of accuracy. There are two eigenfunctions of this 
type which can be formed: y; is symmetric with 
respect to an interchange of atoms a and c¢ and 
corresponds to the bond a—c; y is antisym- 
metric in this respect and corresponds to a bond 
between a and 0 alternating with a bond between 
b and c. With the abbreviated notation: 


i a(1)a(1) B(1)B(1) e(t)a(1) 
( : )-6 1)-Ha(2)a(2) 6(2)8(2) c(2)a(2)|, (5) 
a(3)a(3) 6(3)6(3) ¢(3)a(3) 


> 
where a(1) and @(1) are spin functions of electron 
1 which are eigenfunctions, respectively, for 
Su= 3h/2x, —fh/2r. 


abe abc abc 
sc (00 Jes O00 bs OO ES 
apa Baa aap 
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abe abc 
w=" )-G..): ) 
These correspond to the electronic energies: 

_ SSS Vy rdridradrs _ a°b2c?+ a®bcbe — b2acac — abbcac (8) 

SSS Vididridredrs = 1+ o0?(R) — I00(2R) — Lo0(2R) L007(R) 
o SSS Ul ysdridradrs _ a®b?c?+ b*acac — a*bcbe — abbcac | (0) 

SSS Vsbsdridradts 1+ I0?(2R) —Io0?(R) — Lo0(2R)I00?(R) 

where we have used the notation : Ig9(R) = (1/7) fe~a~edr, 

a*bcbc = (a(1)b(2)c(3) | H| a(1)c(2)b(3)), ete. (10) 


These energy components are plotted in Figs. 1 and 2. 
The energy of the symmetrical state, y;, is shown in Fig. 3, curve V. It corresponds to a 


3 Heitler and London, Zeits. f. Physik 44, 455 (1927). 
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strong repulsion between the three hydrogen 
atoms. 

The energy of the state y; (antisymmetric with 
respect to a and c) is shown in Fig. 3 as curve I. 
The energy is less than for the three separated 
atoms and greater than for a diatomic hydrogen 
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molecule and a hydrogen atom. The lowest 
energy is —53.11 kcal./mole for the separation 
of 2.0 a between neighboring atoms. This is to 
be compared with the energy of the diatomic 
hydrogen molecule calculated by Heitler and 
London and Sugiura, —72.18 kcal./mole. If a 
linear symmetrical molecule really corresponds 
to the activated state, then the activation energy 
for the reaction H+H2—He2+H, will be 19.07 
kcal./mole by this calculation. However this 
approximation for the diatomic molecule is 36 - 
kcal. or 1.6 ev higher than the experimental value 
of —108.6 kcal. Comparing the experimental 
activation energy (5.5 to 7.3 kcal.),® one sees 
that the energy obtained here for the triatomic 
molecule must be about 0.6 ev less accurate than 
that for the diatomic molecule. 

James and Coolidge calculated the energy of 
the triatomic hydrogen molecule for the single 


4Y, Sugiura, Zeits. f. Physik 45, 484 (1927). 

5 It appears that the zero-point energy in the activated 
complex is within one or two kcal./mole the same as that 
for the normal He molecule (6.2 kcal./mole). See A. and 
L. Farkas, Proc. Roy. Soc. A152, 124 (1935); K. Geib and 
P. Harteck, Zeits. f. physik. Chemie Bodenstein Festband, 
849 (1931). 










































internuclear separation R=1.7 ay. For this point 
they obtained the energy of —39.7 kcal./mole 
which is in good agreement with our value for this 
separation. 


2. Heitler-London plus polar states 


A slightly better approximation to the energy 
of the triatomic hydrogen molecule may be ob- 
tained by including the polar states in the eigen- 
function. There are three polar states which are 
antisymmetric with respect to an interchange of 
a and ¢ and which must be considered. y2 corre- 
sponds to two electrons on atom a and one on ¢ 
alternating with two on c and one on a. ys cor- 
responds to two electrons on b and one electron 
alternating between a and c. 4 corresponds to 
one electron on b and two electrons alternating 
between a and c: 


aac cca 
COC 
aBa aBa 
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abb bbe 
V3 _ ( ) + ( ). (12) 
apa aBa 
bce aab 
w=(“*)+/( ) as 
aBa aBa 
The energy components are shown in Figs. 1 and 
2. By letting S;; be the matrix components of 


unity: 
Sis=SIS Vib dridredrs 
and //;; the matrix components of /7: 
Hi;; Sf f VF jd7 dredr3. 


The energies are given as the roots of the secular 
equation: 

We are interested only in the lowest root of this 
equation, the other roots corresponding to higher 
states of the molecule. The composition of the 
eigenfunction which leads to the lowest energy 
value may be written as 








W = ,(Si1)~ Wi tc2( See) Wotc3(S33)~!Wstea( Sas) Wa. (15) 


The c; are the solutions to the set of linear equations : 





The absolute values of the c; are chosen so that V 
is normalized to unity. Fig. 4 shows the c¢; 
plotted as functions of the internuclear separa- 
tion R. 

The energy of the pure polar states 2, 3 and 4 
are shown by curves II, III and IV, respectively, 
of Fig. 3. They correspond to much higher ener- 
gies than the homopolar state. However, these 
states interact and the best linear combination 


Values of coefficients of y's. 






in Bohr radi neighbors in ineor symmetric Hy 


Fic. 4. Coefficients of Y = cyitcaet+ cs t+cas. V and y's 
are normalized to unity. ¥;’s are not orthogonal. 





(S11) $e: = C1 (S11)~? + CoH ie(S22)? + 3H i3(S33)? + 08H i4( S44) 3, i=1, 2, 3,4. (16) 








of the four gives an energy level which is con- 
siderably lower than the homopolar. It is custom- 
ary to designate such improvement by the term 
“resonance energy’’; its physical significance is 
the improvement of the faulty charge distribution 
due to the approximate nature of the original 
eigenfunctions. The best linear combination is 
shown as curve I’ in Fig. 3. The minimum energy 
which this type of eigenfunction yields is — 60.39 
kcal. /mole for R=2.0 ap. 

Weinbaum® calculated the energy of the 
diatomic hydrogen molecule using the same type 
of eigenfunction, i.e., homopolar plus polar 
states, and he obtained —74.02 kcal./mole. 
Comparing the energy of the triatomic molecule 
with this, one finds that the calculated activation 
energy for the reaction H+H»—-H:+H is 
(except for zero-point energy) 13.63 kcal./mole. 
This activation energy is in fair agreement with 
the experimental value of 5.5 to 7.3 kcal./mole. 
Unfortunately the eigenfunctions are still too 


6S. Weinbaum, J. Chem. Phys. 1, 593 (1933). 
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approximate for this agreement to have real 
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significance, since the absolute value of the 
energy of the diatomic molecule with which we 
are comparing the energy of the triatomic mole- 
cule is 34 kcal./mole too large. 


3. The Wang type of eigenfunction 


In the diatomic hydrogen molecule, Wang’ 
obtained a considerable improvement over the 
Heitler-London energy value by varying the 
“effective nuclear charge,’ 2, in the atomic 
orbitals. This same type of variation may be 
carried out for the triatemic molecule without 
introducing any new integrals. Our method is 
similar to that used by Wang and by Coulson.® 

The overlap integral Jo(R) goes into Ipo(p) 
where p=2R. The nuclear attraction and elec- 
tronic repulsion integrals become the same 
integrals as functions of p instead of R except that 
they are multiplied by z. 

The kinetic energy gives rise to three separate 
integrals: 


S a(1)(— (h?/8x2m)V 1?)a(1)dr;—27e?/ 2a, 
Sb(1)(— (h?/82?m)V *)a(1 )dri—[ — 3272 o0(p) +27Joo(p) Je?/ao, 
S c(1)(— (h?/8x?m)V 1?)a(1)d71—[ — 3272 00(2p) +2°S 00(2p) Je?/ao. 


By combining the energy components it is easy 
to show that the energy in units of e?/ao for the 
pure homopolar state 1 has the form 


(EZ) .=2[(E) emit (z— 1)A 11/Si1], 


where Ay is a function of p. 
Now keeping zR=p constant, we can vary z to 
give the lowest energy: 


(0E/dz),=0=(E) 221+ (22-1) A11/Si1, 
Zmin > —}(E) 2=1911/A ut, 
Enin= — Si[(E) emi —Ai1/S11 ]?/(4A 11). 


Enin is the lowest energy that may be obtained 
for a particular value of zR but not necessarily for 
a particular value of R. However, if p is now 
varied to give the lowest value of Emin, the latter 
will agree with the lowest possible value of E 
which can be obtained by varying both R and z. 

This process corresponds to the Wang ap- 
proximation for the diatomic hydrogen molecule. 
In the diatomic molecule the variation of z led to 





an energy of —86.94 kcal./mole which is an 
improvement of 14.8 kcal./mole over that for 
z=1. In the triatomic molecule this variation 
leads to an energy of —56.16 kcal./mole (at the 
point R=1.889 a) and with z=1.059), an im- 
provement of only 3.05 kcal. This is 4.23 kcal./ 
mole worse than for the H.L. plus polar states 
value. See Fig. 3, curve I’’ for the values of Emin. 
The corresponding activation energy for the 
H+H2—He+H reaction which we calculate as 
the difference between diatomic and triatomic 
energy is then 30.78 kcal./mole. 


4. Wang plus polar states approximation 


The variation of the “effective nuclear 
charge,” z, was carried out for the eigenfunction 
in which the polar states are included. The varia- 
tion of z and of the composition of the eigenfunc- 
tion was carried out by successive approxima- 


7S. C. Wang, Phys. Rev. 31, 579 (1928). 


8C. A. Coulson, Proc. Camb. Phil. Soc. (2) 31, 244 
(1935). 
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tions. First, the composition of the eigenfunction 
was taken to be the same as that previously ob- 
tained for the H.L. plus polar states approxima- 
tion; a value of the best z for that composition 
was then obtained. The composition was then 
varied for that value of z. This last step improved 
the energy by only 0.35 kcal./mole at the point 
p=2.0; it was therefore deemed unprofitable to 
go back and readjust z. The energy values are 
shown as curve I’’’ in Fig. 3. The value for the 
Wang plus polar states approximation for the 
diatomic hydrogen molecule is — 92.24 kcal. /mole 
as obtained by Weinbaum. Since the minimum 
value which we find for the energy of the triatomic 
hydrogen molecule to this approximation is 





EYRING AND ROSEN 


— 67.09 kcal./mole for R= 1.841 a and z= 1.087, 
the calculated activation energy then becomes 
25.15 kcal./mole. This activation energy is 17 
kcal./mole larger than the experimental value. 
This approximation is the best that we can ob- 
tain without using more elaborate eigenfunctions 
and calculating new three-center integrals. It is 
not surprising that the variation of the screening 
constant on the three atoms simultaneously gives 
little improvement in the energy since one would 
expect that the central and outer atoms really 
ought to have different screening constants. 
However, such a treatment would involve the 
evaluation of additional integrals and has not 
been attempted here. 


APPENDIX 


Integrals 


The three hydrogen nuclei, a, }, ¢ lie sym- 
metrically on a straight line with a separation of 
R times the first Bohr radius for atomic hydrogen, 
do; and 1 and 2 are electrons. 

<— R~— R— 

a b C 

2 1 

If the molecular eigenfunctions are compounded 
from 1s hydrogen-like atomic orbitals each having 
the same screening constant, all of the trouble- 
some three-center integrals may be carried 
through directly in terms of ellipsoidal coordin- 
ates without resorting to ponderous or slowly 
converging infinite series. The convenient units 
e?/ay= 27.058 ev =623.96 kcal./mole for the en- 
ergy and a= 0.528 X10-* cm for the length, are 
used here. The ellipsoidal coordinates are: 





An= (tint jn)/R’, 

kn = (Tin— rin) /R’, 

¢o, = azimuthal angle of n, 
where i and j are the foci, R’ is the distance be- 
tween them, 7;, and 7;, are the distances of the 


arbitrary point ” toz and j, respectively. 
The basic integral in this work is: 


K(n, ij=r| f from exp (—7ii—171,)d71. 


all space 


(17) 
The integration of K(n, 7j) is simple in ellipsoidal 


coordinates where 7 and j are foci. It may be ob- 
tained as a by-product of Sugiura’s calculation of 


rf fre exp (—1a1— 161 — Yaz — "2)d7 1d 72. 


Schuchowsky® and Gordadse?® have also obtained it. 


K(n, tj) = LKo(R’, dn) +on?K2(R’, dn) J, 


where : 


(18) 


Ko(R’, kn) = (4R’)— [3 (An? +1) D —6An00(R’) +(6Xn —2R’) exp (—R’An)], 
K2(R’, dn) =3(4R’)—[ — (3dn?— 1) D+ On 00(R’) — (6\n+2R’) exp (—R’An)]} 
D=Ioo(R’)[ (log (An+1)/(An—1)) +.Ei(—R’ (Xn —1))]— F(R’) Ei(—R’(n+1)), 


Too(R’) = (4R” +R’ +1) exp (—R’), 
F(R’) =(2R?—R’ +1) exp (+R’). 


® A. Schuchowsky, Acta Physicochimica U.S.S.R. 1, 901 (1934). 


10 G. Gordadse, Zeits. f. Physik 96, 542 (1935). 
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Here Ei(— x) is the integral logarithm, defined by 


fee) 


Ei(-—x)= -{ x! exp (—x)dx. 


z 


Values of Ei(—«x) for x greater than 15 were calculated from the asymptotic expansion of Ei(—+x). 
Two special examples of K(n, ij) are the two three-center integrals which arise from the attraction 
of the electrons to the nuclei. K(a, bc) is K(n, ij) in which R’=R, p,=1, \,=3. K(b, ac) is K(n, ij) 
in which R’=2R, u,=0, A, =1. 

There are four three-center multiple exchange integrals which arise from the mutual electronic 
repulsions. We shall consider them in order. 


L(aa, bc) =f fran exp (—2%a1— 02 —1c2)d71d 72. 
If we integrate over the coordinates of electron 1, L(aa, bc) = K(a, bc) —T(a, bc), 
where T(a, bc) =f +ra2') exp (—2ra2—1%2—Pe2)dr. 


By putting this in ellipsoidal coordinates with a and ¢ as foci and integrating over the azimuthal 
angle @o, 


co 1 
T(a, bc) =2R*f f [(A2— pw?) + R-"'(A—) J exp (—3RA— Ru— R(A2+y?2—1)!)dudy. 
1” -1 
By making the substitution: 
f=u+(d?+y?—1)}, 
u is expressible rationally in terms of ¢: 


p= 3t—3t-1(A?—-1). 
T(a, bc) is now free from radicals: 


foe) A+1 
T(a, bc) = f f [RA -+- + +k'(A)t*] exp (—3Rd—Re)dtdd, 
1 4-1 


where k(A), ---R’(A) are polynomials. This is readily integrated to give: 


817 (a, be) = (e-*®/16)[131/3R+1670+4+2256R+624R?]— (e-** /16)[131/3R+34+160R] 
—[116/9R+116/3+44R+16R2][Ei(—2R)+2 log 2 Je-3” 
+e®®Ei(—8R)[116/9R—116/3+44R—16R?}. 


2. L(bb, ac) = mf fre exp (—2%51—fa2—1c2)d71d 72. 
The integration over the first electron coordinates gives: L(bb, ac) = K(b, ac) —T(b, ac), 
where T(b, ac) = =f +12!) exp (—2fr2—1a2—1e2)dT2. 


Again by using the ellipsoidal coordinates with a and ¢ as foci and integrating over ¢o: 
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ie) 1 
T(b, ac) = aR f (x2 — 2) [1+ R-“(A2+ w2—1)-!] exp (—2RA—2R(A2+y2—1)!) «dud. 
1 1 


The substitution s=\+ (\?+y?—1)} enables \ to be expressed rationally in terms of s, 
A=35—25-"(u?—1), 
and thereby eliminates radicals from the integral, which may be evaluated directly. 


1057(b, ac) =16R*(7+2R*)[Ei( —2R) —2Ei(—4R) ]+e-2*[16R®—8R!+64R3—40R2+52R 


+108—51/R]+e—*"[ —16R°+4R‘*— 58R?+15.5R?—8.5R+33.375+51/R ]. 


3. Liab, a= f fr exp (—%a1—1e1— a2 — Yv2) 7 1d 72. 
The integration over the first electron coordinates gives: 
L(ab, ac) =a f K(2, ac) exp (—tes—tis)drs. 
By expressing this in ellipsoidal coordinates with foci at a and c and integrating over ¢o, 


L(ab, ac) =2R* [ [CoKo(2R, \) + C2K2(2R, d) ]dd, 
~< 
1 
where C= f (\2—u?) exp (— Ru— R(\?+y?—1))dy, 
1 


1 > 
C< f (x?—p?)u? exp (—Ru—R(d?+u2—1)))dy. 


—1 











(20) 


(21) 


The same substitution used in L(aa, bc) makes the integration of Cy and C2 possible. The last inte- 
gration over \ would also be carried through analytically, but the integrand involves so many 
separate terms that it is much easier to carry through the last integration numerically. Cy and C2 as 
well as Ko(2R, \) and K2(2R, \) were evaluated for a large number of values of \ for each value of 
R considered. The final integrand was obtained by accurate curve fitting with functions of the form: 


f(A) exp (—2RA)+¢(d) exp (—4RA), 


where f(A) and g(A) are polynomials. Care was taken that these functions had the proper behavior 
for large values of \ and that they had no extraneous oscillations. The final integration was carried 
out with these functions. The values obtained for the integral are probably correct to five figures 


after the decimal point, in atomic units, i.e., 0.005 kcal. 
4. L(ab, bc) = rf fro exp (- Ta1—~T1d1 —ty2—1e2)d7\d 72. 
After integrating over the coordinates of electron 1, 


L(ab, bc) = ™{Ke, ab) exp (—fr2—1e2)dT2. 


By integrating over the azimuthal angle and putting this in the ellipsoidal coordinates with a and b 


as foci. 
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R3 co 
L(ab, bc) “7. f [MoKo(R, \)+M2K2(R, d)] exp (—2RA)dd, (22) 
1 


1 
where Mo= f (A?—y*) exp (2Ru— 2R(A?+y?—6ur+8)!)dy, 


—1 


1 
Ma= f u?(d2— uw?) exp (3 Ru—FR(A2+y?—6urA+8)!)du. 


| 


Here the substitution 
2=3\—p+(d?+y?—6urA+8)! 


makes pe expressible in terms of 2 in a rational 
manner: 


be= —32+3A+42-1(1—d?2), 


and the integrations over z may be carried out 
analytically. The final integration over A, just 
as in the case of L(ab, ac), might be carried 
through analytically but the large number of 
terms in the integrand at this point makes such a 
task extremely laborious. The My and Mz as 
well as the Ko(R, \) and K2(R, \) were computed 
for a large number of values of \ for each value of 
R considered and the integrand was accurately 
“curve-fitted.”” The curve fitting in this case 


of 


wnization. 


Energy—> Unit twice 





once nearest H atams— 


Fic. 6, Some two- and three-center energy integrals. 





was much more difficult than in the previous case 
and it was necessary to fit a number of sections 
and integrate these separately. The final numeri- 
cal values for these integrals are probably also 
correct to five places after the decimal point 
in atomic units (0.005 kcal./mole). In Fig. 5 the 
four integrals just discussed are plotted against 
internuclear separation. 

Fig. 8 shows the four three-center integrals as 
functions of the internuclear separation R. They 
each approach 5/8 as R approaches zero. The 
value of these integrals for ordinary interatomic 
distances in molecules are of the same order of 
magnitude as those for two centers although they 
become comparatively negligible for large sepa- 
rations. 

Previous to this work, the three-center multiple 
exchange integrals have always been treated by 
expanding 1/712 and at least two of the exponen- 
tial factors in the integrand in terms of a doubly 
infinite series of spherical harmonics about one of 
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Fic. 7. Some energy integrals. 
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the centers." The resulting integrations may then 


11 A. S. Coolidge, Phys. Rev. 42, 189 (1932), used this 
method for HzO. Some work along this line has been done 
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be carried through by making use of the orthog- 
onality properties of the spherical harmonics. 
However, the series which must be evaluated, 
converges slowly except for small internuclear 
distances and the terms are ponderous. Coolidge 
and James have used this method for the inter- 
nuclear separation of 1.7 a) and obtained agree- 
ment with our values to within half a kcal./ 
mole for L(aa, bc), L(bb, ac) and L(ab, bc) but 
their value of L(ab, ac) is some two kcal./mole 
smaller than ours. 

The one- and two-center integrals which occur 
in this treatment have been previously used for the 
diatomic hydrogen molecule.” The integrals for 
nuclear attraction are shown in Fig. 6, those for 
electronic repulsion in Fig. 7. The overlap inte- 
grals and all the combinations of them required 
for H; are given in Fig. 8. 

We wish to express our appreciation to the 
American Philosophical Society for financial 
assistance during the course of this work. 
by Hershtenkorn (Diss., M.I.T., June 1932) and by A. A. 


Frost at Princeton. 
12 See, e.g., N. Rosen, Phys. Rev. 38, 2099 (1931). 
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The potential energies for linear symmetrical configura- 
tions of H;+ have been calculated by the same methods 
used in the previous paper (I). The energies at the mini- 
mum in kcal./mole and the corresponding distances 
between neighboring atoms in Bohr radii for the various 
approximations are: 


Heitler-London — 91.92 at 2.00; 
H.L. plus polar —108.64 at 2.00; 


Wang Approx. — 130.89 at 1.55; 
Wang plus polar —155.35 at 1.52. 


It follows from the use of the variational method that 
the energy of binding of a proton to Hp is greater than 
46.8 kcal./mole. If we use the same Wang plus polar 
approximation for H2 and H;* we obtain for this binding 
energy 63.11 kcal./mole. The reaction H.++H2—H;*++H 
may be exothermal and cannot be more than slightly 
endothermal. 





HERE is experimental evidence that H;+ 
ions are formed whenever Hz is ionized at 


any but the lowest pressures.! This probably 


1H. D. Smyth, Rev. Mod. Phys. 3, 347 (1931); Phys. 





occurs by the following secondary processes : 


Rev. 4, 452 (1925); W. Bleakney, Phys. Rev. 40, 496 
(1932) ; 35, 1180 (1930) ; Hogness and Lunn, Phys. Rev. 26, 
44 (1925); Brasefield, Phys. Rev. 31, 52 (1928); Dorsch 
and Kallmann, Zeits. f. Physik 53, 80 (1929). 
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H.++H.2—-H;t+H, (1) 
H*++H.—H;", (2) 
H +H,.t-—H;", (3) 


of which the first is ordinarily the most im- 
portant. In this paper only the linear symmetric 
configurations are considered. The calculation 
of the energy of the triatomic hydrogen ion is 
then very similar to that of the triatomic 
hydrogen molecule treated in the paper immedi- 
ately preceding (1). The same integrals are used 
and the same types of approximations con- 
sidered. The lowest energy states for the tri- 
atomic hydrogen ion are found to be symmetric 
with respect to an interchange of atoms a and c 
in the electronic eigenfunctions as contrasted 
with the H; case. 


1. ENERGY oF Lowest BonpbD STATE 


The simplest type of eigenfunction which we 
can form from the atomic orbitals considered in 
paper I corresponds to a bond between the atoms 
a and 6 alternating with a bond between 0 and c. 


¥i = (2)~*La(1)b(2) +a(2)(1) +(1)c(2) 
+0(2)c(1) JLa(1)6(2) —a(2)(1) J. 
This state corresponds to the energy : 
a*b?+ b’ac+abab+ abbc 
1+Joo(2R)+2I002(R) 





Values for the energy components are given in 
Fig. 1. The energy is shown in Fig. 2 as curve I. 
The lowest value which we obtained for the 
energy was —91.92 kcal./mole. 


2. HEITLER-LONDON PLUs POLAR STATES TYPE 
OF EIGENFUNCTION? 


The energy for the triatomic hydrogen ion 
can be considerably improved if we include in 
the eigenfunction linear combinations of two 
polar states and two homopolar states so as to 
get the minimum energy for the combination. 


. The eigenfunctions which can be formed from 


the atomic orbitals considered, and which are 


_? Heitler and London, Zeits. f. Physik 44, 455 (1927); 
S. Weinbaum, J. Chem. Phys. 1, 593 (1933). 
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Fic. 1. The energy matrix components for H;*. 


symmetrical with respect to an interchange of 
atoms a and c are: (1) bond eigenfunction 
between a and 6 plus a bond eigenfunction 
between b and c; (2) bond between a and c; 
(3) both electrons on the central atom 6; (4) both 
electrons alternately on a and then on c. 


¥2= 3La(1)c(2) +e(1)a(2) JLa(1)8(2) —a(2)8(1) J, 
a= 2-40(1)b(2)[a(1)8(2) —a(2)8(1) J, 
¥s= 3La(1)a(2) +c(1)e(2) JLa(1)8(2) — a(2)B(1) J. 


The calculation of the energy involves the 
same formal processes as in the equivalent 
approximation for the triatomic molecule. The 
energy components are shown in Fig. 1. The 
energy for the pure states 1, 2, 3 and 4 are shown 
as curves I, II, III, IV, respectively, in Fig. 2. 
The energy of the ion is then obtained by the 
solution of a fourth-order secular equation. The 
values are shown graphically in Fig. 2 as curve 
I’. At large internuclear separations there is a 
strong interaction between the two homopolar 
states (1) and (2). It is interesting that for 
distances so large that states ¥; and ye have 
practically attained their asymptotic energies, 
their interaction energy Hy is still appreciable. 
At small separations there is a strong interaction 
between (1) and the polar state (3). Fig. 3 
shows the coefficients c; of the normalized ¥ 
giving the total normalized eigenfunction of 
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Fic. 2. Potential curves for various states of linear 
symmetrical H;*. 


lowest energy. These c; are similar to those 
defined for the triatomic molecule in paper I. 
The minimum energy that is obtained in this 
approximation is — 108.64 kcal./mole at a sepa- 
ration of 2.0a) between neighboring atoms. The 
improvement, it will be noticed, is 16.72 kcal./ 
mole over the corresponding value for the lowest 
homopolar state. 


3. WANG TYPE OF EIGENFUNCTION?® 


It is not difficult to vary the effective nuclear 
charge in exactly the same manner as for the 
triatomic hydrogen molecule in paper I. It should 
be emphasized again that because of the type 
of variation the value for E is the lowest for a 
particular value of p=zR but not necessarily 
for any particular value of R except at that 
point at which it is a minimum for both z and R. 
Thus it will be noticed in the curve I” of Fig. 2 
which shows this energy as a function of R, that 


8S. C. Wang, Phys. Rev. 31, 579 (1928). 


Values of coeticients of y/s. 


/ 
Ostance tn Bohr radi betveen neighbors in inear symmetric KH: 


Fic. 3. Coefficients of V=ciW+copetosystews WV 
and the y;’s are normalized to unity. The y;’s are not 
orthogonal. 


this method of varying z led to considerable 
improvement of the energy for separations near 
R=2.0 which is the most favorable separation 
when z= 1. However, curve I” merges into I’ for 
larger or smaller separations and for these sepa- 
rations this particular way of varying z does not 
improve the energy. The lowest energy which is 
obtained in this approximation is —130.89 
kcal./mole for the separation R=1.546a) and 
z=1.294. This is an improvement of 39 kcal./ 
mole over the lowest bond state approximation. 
It is remarkable that the variation of the 
effective charge produces such a striking change 
in this energy when it made so little difference 
in the energy of the triatomic molecule. 


4. WanG Ptus Ionic STATES EIGENFUNCTION 


By varying the effective nuclear charge in the 
H.L. plus polar eigenfunction the energy of the 
triatomic hydrogen ion is improved considerably. 
Knowing the composition of this eigenfunction 
for z=1, we may find the energy of the molecule 
for this same composition as a function of z. The 
c; are given in Fig. 3. z is then varied until E 
has a minimum value for the particular value of 
p=2R considered. This value of E can be 
improved still further by varying the composition 
keeping z constant. This variation of the compo- 
sition improved the result at the energy minimum 
by 2.87 kcal./mole. The value of the energy as 
a function of R is shown as curve I’” in Fig. 2. ° 
The minimum value of the energy has dropped 
to —155.35 kcal./mole for the separation of 
R=1.525a 9 which corresponds to z=1.312 and 
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zR=2.0a. This is to be compared with the 
energy of the normal diatomic molecule calcu- 
lated in this same manner by Weinbaum,? 
—92.24 kcal:/mole. This calculation gives an 
energy of binding of a proton to the He of 63.11 
kcal./mole. This is very nearly equal to the 
energy of binding of a proton to a hydrogen 
atom, i.e., 64.04 kcal./mole. (The zero-point 
energy is not taken into account in either of 
these values.) Since the true energy of the 
diatomic hydrogen molecule is —108.6 kcal./ 
mole instead of the calculated —92.2, the very 
smallest amount of binding energy which the 
triatomic ion could possibly have relative to a 
hydrogen molecule plus a proton is 46.8 kcal./ 
mole. However, the actual binding energy of the 
proton to the hydrogen molecule will probably 
be larger instead of smaller than 63.11 kcal./ 
mole. The reason for this supposition is that our 
experience with H; was that the accuracy for 
the three-center problem was poorer than that 
obtained for the two-center case by the corre- 
sponding method. In any event, the reaction 
H.++H2—H;t+H should take place with at 
most a small amount of kinetic energy, and it is 


probable (not certain) that it can take place 
without any relative kinetic energy. 


5. REMARKS 


From our calculations we find for the linear ion 
a moment of inertia of 2.1510-*° g cm*. If the 
same distances between neighboring atoms were 
preserved in the equilateral triangle model, the 
latter would have two moments of inertia having 
values of 3 and } of the above. It would require 
an increase of about 40 percent in the distances 
to double the latter moments of inertia. Allen 
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and Sandeman‘ found three band systems in the 
secondary spectrum of a hydrogen filled Wood's 
discharge-tube, which they interpret as due to a 
triatomic hydrogen molecule or ion with a 
moment of inertia of 1.810-*° g cm?. These 
bands are enhanced in the arc spectrum and 
under conditions favorable to the formation of 
H3*. 

The equilatéral triangle model for the tri- 
atomic hydrogen ion was considered by C. A. 
Coulson.’ However, Coulson did not calculate 
the three-center multiple exchange integrals, 
and without them he was not able to place the 
energy of the molecule within any significant 
limits. He believed that the triatomic ion would 
be more stable in the equilateral triangle con- 
figuration than in the linear. His argument was 
based on the fact that the binding energy which 
two electrons exert on three centers is increased 
if the configuration is triangular. While this is 
undoubtedly true, the nuclear repulsions become 
much larger at the same time. For a separation 
between neighbors of 1.5a 9 (which must be nearly 
the best separation for both cases) the energy 
of nuclear repulsion is 1040 kcal./mole for the 
three atoms in a linear symmetrical configuration 
and 1248 kcal./mole for the equilatera! triangle. 
Thus only if the electronic binding for the triangle 
exceeds that for the line by over 200 kcal./mole 
can the former configuration be the more stable. 
This point requires further investigation. 

We wish to express our appreciation to the 
American Philosophical Society for financial 
assistance during the course of this work. 


oy and Sandeman, Proc. Roy. Soc. A114, 293 
(1927). 

5C. A. Coulson, Proc. Camb. Phil. Soc. (2) 31, 244 
(1935). 
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A very general expression is set up for the potential energy of the tetrahedral molecule X YZ;. 
The normal vibration frequencies are given as functions of the constants in this expression. 
The equation for the parallel frequencies is shown to be essentially independent of the shape 
of the tetrahedron, while for the perpendicular case the two angles determining the molecular 
configuration must be taken into account. Connections are given between the general potential 
energy constants and the constants for the special forces used by Dennison and Sutherland. 





UBSTITUTED tetrahedral molecules of the 

type XYZ; form a large and important 
class, with many representatives among deriva- 
tives of methane, silane and other regular 
tetrahedral molecules. This class is of consider- 
able interest in the study of molecular vibrations. 
A great deal of infrared and Raman spectrum 
data is available, from which it should be 
possible to derive definite conclusions concerning 
intramolecular forces. This requires, however, 
the existence of perfectly general expressions for 
the vibration frequencies. Such formulas will be 
derived in this paper. 

The pentatomic molecule XYZ;, in which 
three equivalent atoms are attached to the 
central atom Y, has geometric symmetry corre- 
sponding to the point group C;,. According to 
the well-known results of group theory,':? the 
molecule has three single vibrations of the 
parallel type, belonging to the symmetric repre- 
sentation A,; and three doubly degenerate vibra- 
tions of the perpendicular type, belonging to the 
representation E. All six frequencies are per- 
mitted by the selection rules in both infrared 
absorption and Raman effect. 

If the Hamiltonian is written in terms of 
variables of the proper symmetry, the secular 
equation for the normal frequencies of vibration 
factors into parts corresponding to the two 
representations. One cubic equation derives from 


*Sarah Berliner Research Fellow of the American 
Association of University Women. 

** National Research Fellow in Chemistry. 

*** A large part of the work described in this and the 
following paper was done at the University of Michigan. 

1E. B. Wilson, Jr., J. Chem. Phys. 2, 432 (1934). 

2 J. E. Rosenthal and G. M. Murphy, to be published in 
Reviews of Modern Physics. 


A,, and two identical cubics ensue from the 
variables transforming according to E. Symmetry 
coordinates which lead to the desired factoriza- 
tion have been described by Howard and Wilson. * 
For the purpose of formulating a _ general 
potential function, geometrical symmetry co- 
ordinates? are more convenient, and they are 
employed here. The geometrical symmetry co- 
ordinates satisfy all the conditions set down by 
Howard and Wilson, except that of conservation 
of angular and linear momenta. 

The notation used is closely related to that 
introduced in a study of the vibrations of the 
molecule YX,.4 Thus 6r; indicates the mutual 
displacement between the central atom Y and 
one of the other atoms, while 6q;; represents a 
similar quantity for two peripheral atoms. The 
numbering is shown in Fig. 1. The shape of the 
tetrahedron is taken into account by means of 
two angles. @ represents the angle formed by a 
Y—Z bond with the axis of the YZ; pyramid, 
while ¢ is the angle between a line connecting 
X with Z and the same axis. M is the mass of 
the Y atom, while m and 7 are, respectively, the 
masses of the Z and X atoms. Since the parallel 
and perpendicular frequencies represent separate 
mathematical problems, they are treated in turn 
below. 


PARALLEL FREQUENCIES 


Take as the three independent geometrical 
symmetry coordinates 


3 J. B. Howard and E. B. Wilson, Jr., J. Chem. Phys. 2, 
630 (1934). 
4 J. E. Rosenthal, Phys. Rev. 46, 730 (1934). 
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$1=6r1+6re+5r;, 
So=6Pr4, (1) 
$3 = 6912 +6913 +6q23. 
Then the most general quadratic potential 
function may be written 
2V =Gii5 1? +Go2252? + G35" 
+2(Gi251S2+Gis51S3+Gegses3). (2) 





The kinetic energy written in terms of these 
variables and the masses is 


2T =m(M+3m+n)—'{3(M-+n) sec? 08)" 

+(n/m)(M+3m)s?+[(M-+n) sec? 6/9+ 4m ]s3? 

+2[n sec 08:82—3(M-+n) tan 0(sec 0/+/3) 8183 
—n(tan 0//3)8e83 |}. (3) 


Eqs. (2) and (3) lead by the usual methods to 
the cubic equation for the parallel frequencies.® 


m3 — Rym?d?+ Romdx— R3=0, 


where 


R,=(3/M)(3m cos? 6+ M)Gii+(m/Mn)(M+n)Ge2—(6m/M) cos 6G12+6(4/3) sin 0Gi;+9G3:3, 
R:-=(3m/Mn)(3m cos? 6+ M+n)(Gi:Ge2— G2?) + (9m /Mn)(M +n) (GooG33— Gos”) 
+ (27 cos? 6/M)(3m+M)(Gi1G33— Gis?) + 6(4/3) sin 6(m/Mn)(M-+n) (Ge2Gi3— Gi2Ge3) 
+18(4/3) sin 6 cos 0(m/M)(Gi1Ge2—G2Giz) +54 cos 0(m/M)(Gi3Ge3— Gi2G33), 


and 


R; =27 cos? 6(m, ‘Mn) (3m +M-+n) (G31G22G33+ 2G 12G13Go2 — G; 1Go3” - Go2G}3" = G33Gi2"). 


In making calculations it is convenient to 
substitute \=w*, w being given in cm~', and to 
use the masses in atomic weight units. The 
constants G thus determined are of the dimen- 
sions mass (cm~')?. However, if the constants 
are desired in dynes/cm, it is necessary to set 
\=47°c’w* and to express the masses in grams 
per atom. To transform the force constants 
from the mass (cm~')? units to dynes/cm, it is 
only necessary to multiply by 0.05863. 

Formulas (4) seem to imply that the parallel 
frequencies depend explicitly on the angle @, i.e., 
on the shape of the tetrahedron. This, however, 








is not the case. For if we introduce new constants 
F such that 


Gi,=3 sec? OF ii, 

Go2= Fy, 

G33= tan? 0F\,+ F33—2 tan 6Fi3, 

Gie= (4/3) sec OF i2, 

Gis= — (4/3) tan 6 sec OF 11+(4/3) sec OF 3, 
Go3= —tan OF \2+ Fos, 


it is found upon substitution that the angle @ 
vanishes from Eqs. (4). This linear transforma- 
tion of the constants corresponds to the use of 
the following geometrical symmetry coordinates : 


$1’ =(4/3) sec 0s;—tan 053, Se’=S2, S3’=S3. (6) 


The F;; are then the constants in the general 
potential function written in these new variables. 
We have preferred to use the variables indicated 


5 As a check the problem has also been set up in sym- 
metry coordinates of the type described in reference 3, 
and the formulas for both the parallel and perpendicular 
frequencies have been verified. The two methods represent 
different ways of treating the kinetic and potential energies, 
but are of course mathematically equivalent. 





136 5s &. 


in (1), for then the coefficients in the potential 
energy are susceptible to a simple physical 
interpretation which will be discussed in the 
following paper. 


PERPENDICULAR FREQUENCIES 


While the treatment of this problem is exactly 
analogous to the one used for the parallel case, 
the resultant formulas for the frequencies are 
much more cumbersome. Their dependence on 
the angles @ and ¢ is inherent in the problem 
and cannot be eliminated by any substitutions. 
We use as the geometrical symmetry coordinates 


(7) 


The most general quadratic -potential energy 
function consistent with C;, may then be written 


$4=6re—6r3, S5=5q21—6q34, S6=45qi2—5q13. 


2V = Gaasa? +Go555? + Geese” 
+2(Gas58s+GaeSsSe+Gsesss6); (8) 


while the kinetic energy takes the form 
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2T = (m/to) ( tassa?+bs585? + loose? 
—2(tasSa8st+tagssSettsessse) ], (9) 
where 
ty=3(cot g—cot 0)?+(m/n) cosec? 0 
+(m/M) cosec? ¢+3m(3m—n—M)/Mn, 
tsg=4 (cot? o+3+3m/2n) cosec? 8, 
ts, =4(cot? 6+3+3m/2M) cosec? ¢, 
tes=2(cot g—cot 0)?/9+3(m/n) cosec? 0 
+4(m/M) cosec? 9+3(m?/nM), 
t4,=4(cot ¢ cot 6+) cosec ¢ cosec 8), 
tys=[4 cot g(cot g—cot 6)+3(m/n) } 
Xcosec 6/+/3, 
tss=[4 cot 0(cot @—cot ¢)+3(m/M) | 
Xcosec ¢/+/3. 


The frequencies are given by 


m*\* — Pm?)?+ PomxA—P3=0, with 
P,=(3m sin? 0/M+2)G44+4 cos (8@— ¢)Gas+(3m sin? 9/n+2)G55;+3Geet+2(+/3) sin 0Gae 


+2(y/3) sin Gis, 


P2=3(3m sin? 6/M+2—sin? 6) (GasGe6— Gag?) +2(4/3) (3m sin ¢ sin? 0/M-+2 sin 9 
—2 cos (@—¢) sin 0)(G4sGsg—GasGag) +6(2 cos g cos 6+sin ¢ sin 0)(G4sGee— GasGse) 
+3(3m sin? ¢/n+2—sin? 9) (GssGes—Gse?) + 6m sin? 6/M+4 sin? (0— ¢) 
+9m? sin? ¢ sin? 6/Mn+6m sin? ¢/n ](GasGss— Gas?) 
+2(4/3)[2 cos g sin (@— ¢)+3m sin? ¢ sin 0/n ](GssGas—GasGse), 
P3=3[2 sin? (@— ¢)+3 sin? o(2m/n+3m? sin? 6/Mn—m sin? 6/M—™m sin? 6/n)+6m sin? 6/M | 
X (GssG5sGeo+2GasGaeGse — GasGse? — GssGae? — GoeGas”). 


It should be pointed out that the general 
formulas given here can be applied to special 
types of potential energy functions such as 
valence forces, Urey and Bradley® forces, or 
other types based on simplified physical models. 
To obtain the corresponding equations for the 
frequencies, it is only necessary to find the 
connection between our general G;;’s and the 
special constants and substitute it into Eqs. (4) 
and (11). For example, the relations 


e +i Urey and C. A. Bradley, Jr., Phys. Rev. 38, 1969 
(1931). 





Gi=(1/9)a sec? 6, 

Ge=c, 

G33 = (1/27)(b—2d tan 6+<a tan? 6), 
Gi2=0, 

Giz3= (1/94/3) sec 0(d—a tan 8), 


Ge3=0, 


(12) 


which are derived from the four constant po- 
tential function assumed by Sutherland and 
Dennison,’ lead to the special formula they have 
given for the parallel frequencies of this molecule. 


7G. B. B. M. Sutherland and D. M. Dennison, Proc. 
Roy. Soc. A864, 250 (1935). 
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On the assumption of the invariance of the molecular 
groups YZ; and YX formulas are given connecting the 
potential energy constants of the molecules YZs, XYZz:, 
ZYX:, and ZX,. The validity of the assumption is tested 
for the series CH,4, CICH;, HCCls, and CCl; the con- 
sistency of the results definitely shows it to be a very good 
approximation. The application leads to a determination 
of the general potential functions of all four molecules. 
Frequencies of HCCI; and DCC; are calculated from those 
of CH,4, CICH3, and CCl, and are found to be in excellent 
agreement with the observed data. Predictions are made 
for the frequencies of CICD;. A certain ambiguity exists 


for the constants for the perpendicular frequencies of both 
HCCI; and CICHs, since two different chemically accept- 
able sets will account for the observed frequencies. A 
discussion of the methane constants shows that, while 
individual constants are quite sensitive to various methods 
of determination, certain linear combinations of them, to 
which a definite physical meaning can be attached, are 
practically invariant. The methane constants are used to 
redetermine the constants computed by Sutherland and 
Dennison for ethane; they yield an entirely different set, 
which, however, is equally suited to an explanation of 
the frequencies. 





CORRELATION of the vibrations of 

polyatomic molecules with intramolecular 
forces has been attempted many times on the 
basis of various simplified models, and some of 
these—notably the valence force assumption and 
its various modifications such as the Urey- 
Bradley field—have met with a certain amount of 
success. With reference to these simplified 
models, it has been possible to speak of force 
constants characteristic of certain bonds. How- 
ever, a molecule can only be completely described 
in terms of a general potential function, and it 
is only the corresponding generalization of the 
bond constants that can have any definite 
meaning. 

Spectroscopic observations on deuterium de- 
rivatives have enabled determination of the most 
general quadratic potential functions for certain 
simple hydrides; but in molecules where the 
hydrogens are of minor importance or completely 
absent, isotope shifts are too small to be of much 
help. An enticing alternative method of deter- 
mining general potential functions is the assump- 
tion proposed by Sutherland and Dennison,! and 
others, that the forces remain constant within 
certain molecular groups in passing from one 


* National Research Fellow in Chemistry. 

**Sarah Berliner Research Fellow of the American 
Association of University Women. 

1G. B. B. M. Sutherland and D. M. Dennison, Proc. 
Roy. Soc. A864, 250 (1935). 


molecule to another. The validity of this assump- 
tion is a question of considerable physical inter- 
est, and, if substantiated, should enable calcula- 
tion of the frequencies of molecules for which ex- 
perimental data are not available. 

In the present paper it is proposed to apply this 
method to the series YZs, XYZ3, ZY X3, Y X4, and 
to make a specific numerical test of the assump- 
tion for the special case CHy, CICHs, HCCI;, CCl. 
Stated more explicitly, our assumption means 
that the potential energy for the YZ; group is the 
same whether it is attached to an X or to another 
Z atom; and furthermore, that the force between 
Y and X is not changed by the nature of the 
other atoms in the molecule.? While there is no 
reason for supposing that the methyl group in 
methyl chloride is identical with that in methane, 
or that the CCl; radical remains entirely un- 
changed in passing from carbon tetrachloride to 
chloroform, it is believed that such an assumption 
is not far from the truth. Modern valence theory 
supports this idea, and experimental data appear 
to verify it. Thus the recent electron diffraction 
experiments of Sutton and Brockway,’ in contra- 
distinction to earlier observations, indicate that 


2 This does not imply that the force between Y and X is 
the same as in the diatomic molecule YX. However such 
an extension of the assumption is not too unreasonable, 
and it might be worth while to check it by comparison 
with experimental data. 

3L. E. Sutton and L. O. Brockway, J. Am. Chem. Soc. 
57, 473 (1935). 
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the Cl—C—Cl angle in HCCI; is very near the 
tetrahedral value, and Johnston and Dennison‘ 
conclude, from the moments of inertia, that the 
CH; group has the same dimensions in methane 
and in the methyl halides. Because there is no 
obvious connection between potential functions 
and dipole moments, the difference between the 
dipole moments of chloroform and methyl chlo- 
ride does not seriously indict the assumption. 


CONNECTIONS BETWEEN XYZ3 AND YZ; 
CONSTANTS 


An analytic formulation of the assumption 
leads to a functional relationship between the 
constants’ A, B, C, D and E of YZ,, the con- 
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2V=(9A /16+4B+4E—3D)s,2+(2/3)(A/4+4B—2D)s32+2(2/3)1(5D/2—3A /8—4B)s15; 
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stants® Gi, Gu, «++ Gee of XYZ, the constants 
Gu, --- Gog of ZYXs3, and the constants A, B, C, 
D and E of YX. 

We start by separating the quadratic potential 
function of YZ, into three parts V;, V2, V3 cor- 
responding to forces within the YZ; group, the YZ 
group, and the interactions between the two. 
Thus 


V = Vi(6r1, 672, 573, 5gi2, 6913, 6G23) 
+ V2(6rs2)+V3, (1) 


where V; is a function of displacements of the 
type 6qgu, and also contains terms linear in 67;. 
For the notation, see I. In terms of the variables 
$1*++S,, the potential energy of YZ, becomes’ 









+(34/8)s2+(A/16+B+C—D/2)s62+2(3/8)(D—A /4)sise, 


2V2=(A/16+4B+4E+D)s.’, 


(2) 


2V3=(A/16+B+C—D/2)s;2+2[(4B+4E—3A /16—D)s152+(2/3)(A/8—4B+D/2)s983 


By comparing Eqs. (2) with Eq. (2) of I, it is seen 
that the assumption stated above, which also 
postulates invariance of the tetrahedral angle in 
the group YZz, leads to the following relations: 


Gi=9A/16+4B+4E-—3D, 
Gi3=(—3A/8—4B+5D/2)(2/3)}, 
G33 = (2/3)(A/4+4B—2D), 
Gu=3A/8, (3) 
Gis=(D—A/4)(3/8)!, 
Gos=A/16+B+C—D/2, 
Goo=A/16+4B+4E+D. 
Similar relations for the quantities G;; may be 
obtained by interchanging A and A etc. No con- 
nection is expected between the constants Gy, 


Gos, Gas, Gss, and Gsg and the coefficients in the 
expression for V3; for YZ,, since, obviously, the 


4 3) Johnston and D. M. Dennison, Phys. Rev. 48, 868 
(1935). 
5 Cf. J. E. Rosenthal, Phys. Rev. 45, 538 (1934). 


+(3/8)}(D—A/4)sys;+(A/16+B-—C—D/2)s55¢ }. 





forces between an X atom and a YZ: group are 
different from the forces between the same group 
and a Z atom. The expression for V3 has, how- 
ever, been given explicitly to permit a compari- 
son of the magnitude of the interaction terms in 
YZ, and XYZ;. It may be stated in this connec- 
tion that the applicability of the formulas (3) 
(with the exception of the expression for Gye) is 
not limited to the X YZ; molecule but extends to 
other molecules containing the YZ; group. 

Substitution of relations (3) into the expres- 
sions for the parallel and perpendicular frequen- 
cies (Eqs. (4) and (11) of I), where 6 has been 
set equal to its value for a regular tetrahedron 
leads to Eqs. (4) and (5) given below. For the 
parallel frequencies 


m®\3 —m*)\?R,+mrR2— R3=0, 
with 


6 Cf. the preceding paper on the vibrations of the XYZ; 
molecule. This will henceforth be referred to as I. 

7It should be pointed out that the variables s are not 
geometrical symmetry coordinates for the molecule YZ:, 
since they only satisfy the symmetry requirements of the 
group C;,, and not those of the group 74. 
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Ri=m(1/M+1/n)(A/16+4B+4E+D)+(1+3m/M)(3A/16—D)+4(1+m/M)B 


+4(3+m/M)E—2(m/M)Gi, 


R2= —4m(1/M+1/n)[9Ge3?/4+2Gi2G23(3/2)!]+ (1 +3m/M)(AB/2 —D?/2+2AE 
+32EB—16ED) —4(m/M)[4G12B —Gi2D —G2;(8E—4B+4+3D/2)(3/2)*] 


—m(3/M+3/n+m/Mn)G22+(m/Mn)(A/16+4B+4E+D)[(M+n+3m)(3A /16—D) 


(4) 


+4(M+m+n)B+4(3M+3n+m)E], 


R3=m(1/M+1/n+3m/Mn)[(A/16+4B+4E+D)(AB/2—D?/2+2AE+32EB—16ED) 
— 3G23°(9A /16+4B+4E—3D) —(G12?/2)(A +16B—8D) —(3A/2+16B—10D)G12G23(3/2)!]. 


The corresponding expression for the perpendicular frequencies is 


m*\3 — m?)\2P,+md\P.—P;=0, 


where 


P,=(3m sin? ¢/n+2)G5;+(4/3)(2v2 sin g+cos ¢)Gy,+2(3)! sin oG5¢ 


+(7/16+m/M)A+3B+3C+D/2, 


P2=(5/4+3m/M)[A(B+C) —D?]+3(3m sin? ¢/n+1+ cos? ¢)[(B+C)G5;—Gse? ] 
+sin g[(v2/3+v2m/M—cot ¢/12)A+4(cot ¢+v2)(B+C) —(7v2/3+4v2m/M 


+2 cot ¢/3)D]Gi,+ (3)! sin gL (2m/M+2/3—- 


v2 cot ¢/12)A —(v2 cot ¢+2)D ]Gi¢ 


—8 sin® g[ (2 cot e—v2/2)?/9+m?/Mn+3m/4n+2m(1+cot? ¢)/3M]G,;? 
+4(3)3 sin? gL (1/3) cot p—(2/3)v2 cot? g—V2m/n ]GasGig 
+3 sin*g[(17 cot? ¢/72—5v2 cot ¢/36+17/1444+2m/3M+2m cot? ¢/3M+7m/16n 
+m?/Mn)A +34 cot? g—43v2 cot g+m/2n—1/2)D ]G;;, 
P;=sin? g[(2v2 cot g—1)?/3+9m/n+8m/M+4m(3m—M —n)/Mn+8m cot? ¢/M] 
X[}(AB+AC—D*)G;;+(D—A/4)GasGs6(6)'! — (A /8+2B+4+2C— D)Gy;? — 3A G5¢2/4 ]. 


DISCUSSION OF THE METHANE CONSTANTS 


All the theoretical expressions given above 
refer to the w’s, the frequencies with infinitesimal 
amplitudes. Since a sufficient number of combina- 
tion bands is not known for any of the molecules 
treated to enable us to derive the w’s from the 
experimentally observed frequencies v;, the only 
way to proceed is to set v;=w;, thus introducing 
into the calculations an error of the order of 
magnitude of the anharmonicity corrections. 
However, since the anharmonicity is usually of 
the same sign, we may expect at least a partial 
cancelation of the errors, with a final result cor- 
rect to about five percent. The constants dis- 
cussed below are so adjusted as to give the best 





possible fit for the v’s. With these reservations we 
pass to the discussion of the methane constants. 

The general formulas for the fundamental vi- 
brations of methane 5 involve five force constants; 
and, with only four frequencies, an additional re- 
lationship is necessary for evaluation of the con- 
stants. This was obtained from the interaction 
between vibration and rotation by Dennison 
and Johnston,* who determined the constants 
and used them in the formulas of Rosenthal 
to calculate the frequencies of the deuterium 
derivatives of methane. They obtained two alter- 
nate sets of constants and gave two sets of fre- 


8D. M. Dennison and M. Johnston, Phys. Rev. 47, 93 
(1935). 
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quencies. The infrared absorption frequencies of 
CH;D observed by Ginsburg and Barker® were 
found to agree well with one of the alternate sets 
calculated.!° The observed frequencies may be 
used to recompute the five methane constants." 
However, in the following the constants derived 
by Dennison and Johnston have been used, since 
they reproduce the CH;D frequencies as well as 
the set calculated from the best individual ob- 
served frequency (2205). These constants are 


A=7.67 X10° dynes/cm, B=0.4760, C=0.3409, 
D=1.2779, E=0.1336. 


They are based on the methane frequencies 


Vvo= 1520, 
V4 1304, 


y= 2915, 
v3 = 3014 


for which the assignments are quite certain. 

As a matter of fact the two sets of constants 
are not, for our purposes, as widely divergent as 
they at first appear to be. While the individual 
constants calculated by the two methods may 
differ by as much as 13 percent, certain linear 
combinations of them which appear prominently 
in our potential functions are found to be prac- 
tically invariant. Thus A/16+4B+4E+D is 
equal respectively to 4.88105 dynes/cm and 
4.92, a difference of less than one percent. This 
particular quantity on our assumption corre- 
sponds to the constant of the C—H bond. 

With the values of the methane constants 
known, it is possible to calculate the three con- 
stants a, b, d used by Sutherland and Dennison! 
to characterize the parallel vibrations of the CH; 
group. The analytic relationships may be found 
directly, or by way of Eqs. (12) I and Eqs. (3) of 


(1938) Ginsburg and E. F. Barker, J. Chem. Phys. 3, 668 

10 Recently MacWood and Urey (J. Chem. Phys. 3, 650 
(1935)) have given the Raman spectrum of CH;D. Their 
results agree moderately well with the infrared absorption 
spectrum except for an intense line at 2315 cm™, which 
they interpret as an overtone. The appearance of a strong 
overtone in the Raman effect is most unusual, but the 
force constants derived on the assumption that it is a 
fundamental are incompatible with the other CH;D fre- 
quencies and the parallel frequencies of ethane. Resonance 
with the 2200 fundamental may help in explaining the 
intensity of this overtone. 

There is an error in the constants calculated by 
Ginsburg and Barker. For instance the values derived 
from CH;D 2205 should be A=8.69X10° dynes/cm, 
B=0.437, D=1.33 instead of A = 7.30, B=0.430, D=1.33. 
Whenever their constants are referred to, it is the corrected 
values that are implied. 


this paper. They are 


a=9A/16+4B+4E-—3D, b=8B+32E, 
d=(v2)(8E+3D/2—4B). 


On substitution of the methane constants the 
values a = 3.67 X 105 dynes/cm, b= 13.95, d=3.59 
are obtained. These are to be compared to the 
values a=8.24, b=8.01, d=+5.07 which were 
chosen to fit ethane by Sutherland and Dennison; 
or a=8.47, b=0.83, d=+5.39 chosen to fit 
methyl chloride. In spite of the marked differ- 
ences, the two sets of constants, when used in the 
formulas of reference 1, fit the observed parallel 
frequencies of ethane and methyl chloride equally 
well, and lead to the same values of the X —CH; 
bond constant c. This is another example of the 
situation mentioned above—only certain com- 
binations of the constants (here ab—d? for in- 
stance) are important in the formulas for the 
frequencies. The sensitivity of the individual 
constants is due largely to the manner in which 
the potential function is defined. However it 
should be pointed out that the CH; constants 
a, b, d derived from the methane constants are 
more reasonable physically, since they indicate 
larger stretching and smaller bending forces than 
do the constants derived by Sutherland and 
Dennison. 


EXPERIMENTAL DATA 


The frequencies employed are listed in Table I. 
Those for CICH; are from infrared absorption of 
the gas.!?: 3. 4 731.5 cm is the weighed mean of 
the two isotopic lines for this vibration. 2924 is 
the average of the two observed lines 2880 and 
2967, taken under the plausible assumption of 
Adel and Barker" that the pair is produced by 
resonance between v; and 2¥,=2920. There are 
few data for chloroform or carbon tetrachloride 
in the gaseous state, but since the frequencies of 
methane and methyl chloride are those for the 
free molecules, it has been thought advisable to 
correct the observed Raman spectra for the 


liquids HCCI, and CCl, in order to secure fre- 


12 i‘ Adel and E. F. Barker, J. Chem. Phys. 2, 627 
1934). 
(1934) H. Nielsen and E. F. Barker, Phys. Rev. 46, 970 
934). 
4 E. F. Barker and E. K. Plyler, J. Chem. Phys. 3, 367 
(1935). 





yf Fr AS wh mec ad AS 


a 


VIBRATIONS 


TABLE I. Vibration frequencies in cm, assignments, and 
corrections which have been added to secure 
gaseous frequencies. 








CICHs 
2924 
1355 . 
731.5 
3047, 
14605. 


1020 


HCCls 


43033 15 

683 (16) 
‘376 (10) 
= (10) 
» 776 (16) 


(10) 


DCCl; CCla 
2267 (11) », 467 (8) A, 
666 (16) vw. 222 (5) E 
376 (10) 

916 ( 8) 

753 (16) 


271 ( 9) 








3 318 (5) 
782 7 


\\ 
~271 








quencies for the gases throughout. The alteration 
‘in passing from liquid to gas is usually a shift 
toward greater frequency,’ and that has been 
the uniform assumption in applying corrections. 
On the basis of the known displacements for cer- 
tain CH, and CICH; lines, as well as for 782 of 
CCl, and 3033 of HCCI;,'° small corrections to 
the Raman spectra of the liquids!’ have been 
made. In doing this, the approximate modes of 
vibration for these molecules described by Kohl- 
rausch!* have been employed to characterize lines, 
and the dipole moments of molecular groups have 
been accepted as the chief cause of the displace- 
ments. The corrections, which have been added 
to the observed spectra for the liquids to give the 
values listed, are also shown in Table I. Those not 
in parentheses are based upon direct experimental 
knowledge. Actually these corrections do not 
affect the calculations or conclusions to any ex- 
tent, and their inclusion is a matter of minor 
importance. The frequency 782 for CCl, is the 
mean of the two resonating components!’ 768 
and 797, observed in infrared absorption by 
Schaefer and Kern.”° The. assignments of fre- 
quencies are rendered quite certain by the 
observed character of the bands in the case of 
infrared absorption,'*: '* and by polarization meas- 


% See however W. Buchheim, Physik. Zeits. 36, 694 
(1935), who shows that a shift in the other direction is 
possible. 

© Dr. Ginsburg kindly measured the infrared absorption 
“ chloroform vapor, and obtained the value 3033 cm= 
or v4. 

17 Data for HCCI; from Kohlrausch, Der Smekal-Raman 
Effekt (1931), p. 116. CCl, ibid., p. 305. DCCI; R. W. 
Wood and D. H. Rank, Phys. Rev. 48, 63 (1935); also 
Redlich and Pordes, Naturwiss. 22, 808 (1934). 

18K. W. F. Kohlrausch, Zeits. f. physik. Chemie B28, 
340 (1935). 

9 Cf. Horiuchi, Zeits. f. Physik 84, 380 (1933). 

°° Schaefer and Kern, Zeits. f. Physik 78, 609 (1932). 
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urements" for the Raman lines. For CCl, the 
assignments are indicated in the table by the 
representation of T4 to which the vibration cor- 
responds. 

The perpendicular frequencies depend upon 
the molecular geometry, and it is necessary to 
know the relative bond distances before compu- 
tations can be carried out. The values C—Cl 
= 1.77,22 C—H=1.08 78 have been used. 


THE PARALLEL FREQUENCIES 


We are now ready to apply the formulas de- 
rived above. In doing this, it is sometimes found 
that several alternate solutions for the constants 
are possible. The criterion for choosing between 
these, in addition to that of general self-con- 
sistency, is that the constants shall be physically 
reasonable. Thus large negative values, and those 
with large imaginary parts are considered suspect. 
Small imaginary parts are not sufficient cause for 
the rejection of constants, for they may be due 
to the uncertainties arising from anharmonicity, 
etc. Only small values of Gy and G23 are to be 
expected, for physical reasons. 

The first application is to methyl chloride. 
When the methane constants are substituted in 
formulas (4), it is possible to solve for the three 
remaining constants Geo, Giz and G3, by using the 
values of the three parallel frequencies of methyl 
chloride as given in Table I. Throughout unbarred 
constants refer to methane and methyl chloride, 
while barred constants refer to carbon tetra- 
chloride and chloroform. The masses used are 
C=12.00, H=1.008, Cl= 35.46.24 The formulas 
for the frequencies of methyl chloride lead to a 
quartic equation for the constants. Of the four 
sets of constants which result, only one is satis- 
factory: Ge=3.61X10° dynes/cm. G»2=0.560, 
G2; = — 0.030. The other sets have large values of 
Gy, and for them G2». does not even approxi- 
mately equal the bond constant c of Sutherland 


*1Cabannes and Rousset, Ann. de physique 19, 229 
(1933). L. Simons, Soc. Sci. Fennica, Comm. Phys. math. 
6, 1-58 (1932). Cf. also reference 18. 

2. Pauling and L. O. Brockway, J. Chem. Phys. 2, 
867 (1934); Cosslett and de Laszlo, Nature 134, 73 (1934); 
reference 3. 

*3 Mark, Zeits. f. Electrochemie 40, 413 (1934); L. 
Pauling, Proc. Nat. Acad. Sci. 18, 293 (1932); cf. also 
references 8 and 9. 

24 Use of this value for the mass of Cl is an approxima- 
tion, but since the isotope shift for chlorine is very small, 
it is quite adequate for the present purposes. 
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TABLE II. Potential energy constants. The values given 
under the headings CCl, and CH, are the constants in the 
potential energy functions of these molecules when they 
are treated as though they were of the type XYZ;. The 
numbers are given to the number of decimal places used 
in the calculations even though no significance can be 
attached to the last figure. All quantities are given in 
dynes X 10°/cm. 








CHa CICHs HCCls CCl 


Gur 3.639 3.639 2.850 2.850 
Coss 4.916 3.610 4.916 3.610 


G33 0.844 0.844 0.819 0.819 


0.442 0.560 “_ 1.421 


1.295 —1.295 \ —1.110 — 1.110 


0.250 —0.030 —" —0.661 


2.876 2.876 1.357 1.357 


0.628 c’ 0.570 
0.657 0.892 0.555 0.563 


0.657 0.657 » ~ 0.563 0.563 


— 0.360 —0.155 
—0.392 ~0.735 0.072 —0.161 





— 0.392 —0.392 —0.161 —0.161 


c 0.013 0 
0.657 d 0.367 0 0.563 








and Dennison! (c = 3.44 for CICH;) Furthermore, 
with these other sets it is not possible under the 
present assumption to correlate the frequencies 
of carbon tetrachloride and chloroform. The com- 
plete set of constants for CICH; is given in 
Table II. 

With the value Goo= 3.61, the relationship be- 
tween G2. and A, B, D, E expressed in Eqs. (3), 
and the four known frequencies, it is possible to 
obtain the five constants in the general potential 
function of carbon tetrachloride. Thus the values 
A=3.619, B=0.402, C=0.256, D=0.643, 
E=0.283 are found. Another possible set is 
A=2.53, B=0.840, C=0.256, D=—1.04, E 
=0.283, but this has been rejected because it is 
incompatible with the data for chloroform. The 
potential function for CCl, determined in this 
manner shows certain interesting similarities to 
that for CHy, as may be seen by comparing the 
values of A, B, C, D, E. Urey and Bradley”> have 
given a special potential function for carbon 
tetrachloride. Their assumption requires’ D 
=10B—2C—8E. This relation is not satisfied 
however, for 10B—2C—8E=1.243, while 
D=0.643. The two potential functions seem, 
therefore, to give rather different descriptions of 
the molecule. 


931) C. Urey and C. A. Bradley, Phys. Rev. 38, 1969 
(1931). 


The parallel frequencies of HCCl; may now be 
investigated. As three frequencies and only two 
unknown constants are involved, a test of the 
present assumption is possible. Substitution of 
the frequencies in Eqs. (4) (in this case barred and 
unbarred constants are interchanged) leads to 
equations with imaginary roots. However the real 
values Gis = —0.482, Gos =0.407 give a very close 
reproduction of the frequencies, as is evident on 
examination of the values under a in Table III. 
If, following Sutherland and Dennison’s treat- 
ment of the parallel frequencies of methyl chlo- 
ride, the further assumption that G12. and G23 are 
zero is made, the results shown under 0 in the 
table are obtained. These also agree very well 
with the experimental values, and there is in 
fact little to choose between the two alternatives. 
The results with the second are particularly grati- 
fying, since they represent frequencies for chloro- 
form which have been calculated entirely on the 
basis of data from methane, methyl] chloride and 
carbon tetrachloride. The insensitivity of the fre- 
quencies to the values of Giz and Gy; is due to a 
partial concelation of terms involving these con- 
stants in the equation, and to the fact that the 
constants themselves usually occur with smaller 
coefficients than do constants of type Gu. 


TABLE III. Experimental and calculated frequencies 
for HCCls;. 








CALC. a CaLc. b Oss. (corr.) 
Giz = —0.482, Gos =0.407 Gie, Gos =0 


3034 3010 3033 
675 681 683 
380 388 376 











THE PERPENDICULAR FREQUENCIES 


Eqs. (5) for the three perpendicular frequencies 
of both HCCI; and CICHs; involve three unde- 
termined constants and no further test of the 
major assumption is possible. However, the con- 
stants may be calculated and may then be used to 
compute the frequencies of the corresponding 
deuterium compounds. In each case algebraic 
elimination leads to an equation of the fourth 
degree, and in each case two of the four sets of 
constants can be rejected for physical reasons. 
However, there remain two alternates for Gs, 
Gis, Gsg and two for Gs5, Gis, Gse. These are all 
listed in Table II, as there seems to be no way at 
present of rejecting any of them. It is to be noted 
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in this connection that, while all constants for the 
perpendicular frequencies should be moderately 
small (less than 1 X 105 dynes/cm), those for cross 
terms, such as G;;, need not be small compared 
to those for squared terms in the potential, such 
as G;;. This is simply a consequence of the way 
the potential function has been defined (cf. 


Eq. (8) I). 


FREQUENCIES OF THE DEUTERIUM DERIVATIVES 


It is a simple matter to calculate the frequen- 
cies of deuterochloroform and deutyl chloride 
with the constants which have been determined. 
Substitution of the constants, with the mass of 
deuterium set equal to twice that of hydrogen, in 
Eqs. (4) or (5) leads to a cubic equation; the 
three roots are the desired parallel or perpendicu- 
lar frequencies. The values thus computed are 
listed in Table IV. c and d and c’ and d’ correspond 
to the two alternate sets of constants in Table IT. 
The frequencies computed for deuterochloroform 
agree very well with those observed by Wood 
and Rank,!’ and this agreement is a further 
substantiation of the major assumption. For the 
perpendicular frequencies the agreement is so 
much better with the set of constants d’ that c’ 
could be eliminated on this basis. 

The frequencies of CICD; have not yet been 
observed. The predicted values are given in 
Table IV. Here the two alternate solutions for the 
perpendicular constants also yield frequencies 
which differ considerably, and experimental ob- 
servations should permit elimination of one of 
the sets. 
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TABLE IV. Computed and observed frequencies of deuterium 
derivatives in cm=, 








CICD; 
Calc. 


DCCl; 
Calc. 


Obs. (corr. 
to gas) 





b 
2219 
659 
386 


d’ 

919 916 
745 753 
269 271 


2081 
1110 
645 


2267 
666 
376 











CONCLUDING REMARKS 


The assumption that potential functions within 
molecular groups are preserved has been vali- 
dated insofar as it has been possible to test it. 
It has led to consistent results, and has enabled 
determination of the general quadratic potential 
functions for methyl chloride, chloroform and 
carbon tetrachloride. However, a word of caution 
is necessary with regard to the potential con- 
stants given. Because of anharmonicity, incom- 
plete knowledge of all the gaseous frequencies and 
the limitations of the primary assumption, these 
constants cannot be exact. In particular those for 
cross terms in the potential function are to be 
considered uncertain, since they are quite sensi- 
tive to small changes in the frequencies. But in 
spite of these reservations it is felt that when the 
sets of constants are used to calculate frequencies 
for molecules containing CH; or CCl; groups 
they will yield good results. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


The Relative Values of the Four Butane-Butene-Hydrogen 
Equilibrium Constants! 


Eidinoff and Aston* have given a very convenient method 
for calculating the rotational entropy of molecules with 
internal rotating groups. This method, after some further 
development, has been used to calculate the following 
values for the translational plus rotational entropies at 
25°C. 

n-butane 75.366 2-butene, cis 


isobutane 71.249 2-butene, trans 
1-butene 72.98 isobutene 


The value given for 1-butene is an estimate which should 
be reliable, but which will subsequently be replaced by an 
exact value. An accurate calculation of vibrational en- 
tropies for these molecules is not yet possible. It is probably 
a good approximation to assume that the vibrational 
contributions for m- and isobutane are the same, and also 
that the contributions for the four butenes differ only in 
the part coming from torsion about the C=C bond. This 
part may be estimated by taking the fairly reliable figure 
of 790 cm“ for the torsion frequency in ethylene to de- 
termine the force constant, and by using the ratio 21.5 : 1 
for the moments of CH; and H about the C=C axis. 
A simple calculation gives 572 cm~ for the torsion fre- 
quency in 1-butene and isobutene, 244 cm™ in both 2- 
butenes. The entropies calculated in this way may be 
combined with heats of hydrogenation determined by 
Kistiakowsky, Ruhoff, Smith and Vaughan.* The most 
direct calculation is that of the cis-trans equilibrium, since 
here the torsion cancels out. The equilibrium constant 
is [cts ]/[trans]=exp [—949/RT+0.425/R] or 0.61 at 
400°C. The experimental value of Frey and Huppke‘ is 
0.64. The 1-butene, transbutene calculation is more in- 
volved; the rotational entropy of 1-butene is 5.179 e.u. 
higher than that of transbutene; the torsion correction of 
1.67 e.u. lowers this to 3.509 e.u. The heat content of 1- 
butene is 2720 cal. higher at 82°C; the torsion correction 
gives 2617 cal. at 400°C. Hence [1-butene]/[transbutene ] 
=0.49, against the experimental value of 0.56. The ratio 
of the dissociation constant for isobutane to isobutene 
against that of n-butane to transbutene depends upon the 
rotational entropy differences between n- and isobutane, 
and between trans- and isobutene, the same calculated 
torsional entropy and heat content differences as in the 
preceding example, and the hydrogenation heats. The iso- 
butane dissociation involves a 2.862 e.u. greater entropy 


not later than the 15th of the month preceding that of the 
tssue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


increase and a 695 cal. greater heat input at 400°C. Hence 
Kiso/K trans = 2.51, in somewhat fortuitous agreement with 
2.56, the experimental value of Frey and Huppke. The 
experimental determination of [1-butene ]/[transbutene ] 
varies less with temperature, and that of Kiso/K trans Varies 
more, than the theoretical values. Frey and Huppke have 
stated, however, that the accuracy of their data would 
not justify a calculation of AH from the temperature 
coefficient. The agreement of calculated and observed 
values is thus as good as could be expected. 

It is of interest to note the extremely simple basis of 
these calculations. The bond distances and angles are 
semispectroscopic in origin, but might have been de- 
termined entirely from x-ray diffraction. The individual 
translation-rotation entropies depend upon the quantum 
constant h, but the entropy differences actually used are 
ratios of purely classical phase-space volumes. The treat- 
ment of the torsjonal vibration is the only part of the 
calculation that Gibbs or Boltzmann would have been 
unable to make. 

A more complete account of this work is contained in two 
papers, of which one has been, and the other shortly will 
be, submitted to the Journal of Chemical Physics. 

Louis S. KasseEv*® 


Physical Chemistry Section, Pittsburgh Experiment 
Station, 
U. S. Bureau of Mines, Pittsburgh, Pennsylvania, 
December 20, 1935. 


1 Published by permission of the Director, U. S. Bureau of Mines. 
(Not subject to copyright.) 

2 Eidinoff and Aston, J. Chem. Phys. 3, 379 (1935). 

3 Kistiakowsky, Ruhoff, Smith and Vaughan, J. Am. Chem. Soc. 57, 
876 (1935). 

4 Frey and Huppke, Ind. Eng. Chem. 25, 54 (1933). 

5 Physical chemist, Pittsburgh Experiment Station, U. S. Bureau of 
Mines, Pittsburgh, Pa. 


The Mutarotation of Glucose in H,O—D.O Mixtures 


In an earlier communication! it was shown that the 
dependence of the specific velocity constant for the mu- 
tarotation of a-d-glucose could be described satisfactorily 
by the equation 


kobs = kH,0CH,0+kupOCHDO+kp,0CpD,0. (1) 


Eq. (1) furnishes evidence for the kinetic influence of the 
decrease in basicity of the solvent with increase in D con- 
tent but implicitly denies any influence of the isotopic ex- 
change which takes place in the substrate. 
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TABLE I. 








Pw  k(10)5 


0.050 976 
0.056 979 
0.112 925 
0.142 891 
0.145 869 
0.179 826 
0.190 829 
0.265 762 
0.369 677 
0.384 653 
0.404 654 
0.416 644 
0.507 560 
0.575 526 
0.624 496 
0.652 467 
0.663 452 
0.831 356 
0.890 324 
0.895 328 
0.960 293 
0.960 293 


(H20) 


50.03 
49.91 
43.82 
40.96 
40.69 
37.61 
36.64 
30.38 
22.70 
21.69 
20.35 
19.57 
14.17 
10.68 

8.50 

7.31 

6.88 

1.81 

0.79 

0.71 


(HDO) 


5.20 
5.78 
10.76 
13.14 
13.36 
15.74 
16.46 
20.68 
24.52 
24.88 
25.34 0. 
25.56 0.765 
26.23 0.540 
25.74 0.414 
24.66 0.345 
23.94 0.303 
23.58 0.292 
15.06 0.120 
10.62 0.073 
10.14 0.070 

0.025 

0.025 
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SS ae 
T= Py Ea. b + Po (Ba. 2) =(1040 —kops) /(1040 -273); 
0.83/Rw 


1+0.83/Rw- 


Rg (Col. 6) = 


Pg (Col. 8) = 


If one adopts as a working hypothesis the alternate view 
that the decrease in velocity is due only to replacement of 
H by D in that position in the glucose molecule which is 
involved in the mutarotation and does not depend kineti- 
cally upon the relative abundance of H,O0, HDO, and D.O, 
‘we may write 


hobs = ku,0 — (kH,O—kp,0) Po: (2) 


P, denotes the fraction of a-d-glucose in which the OH 
on the terminal carbon atom has become OD. This ex- 
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change may be assumed as more rapid than mutarotation. 
This hypothesis may be tested readily by calculating from 
Eq. (2), values of P, over the range HXO—D.O, using the 
experimental values of the k’s. P, is readily converted into 
R,, the ratio of the concentrations of heavy glucose and 
light glucose, which may be combined with R,=(HOH)/ 
(HOD) to give the exchange constant 


_ (D—G)(HOH) _ 
~ (H—G)(HOD) — 


R,Rw (3) 

The earlier data! together with more recent results have 
been combined in Table I. P,.=As/0.1079 is the excess 
specific gravity ratio. The K’s are constant within experi- 
mental error. The average K=0.83 agrees exactly with 
that for tetramethyl glucose found directly from exchange 
experiments? but disagrees with the value for glucose (0.69). 
This is to be expected since the exchanging hydrogen atom 
in tetramethyl glucose is the only one which is involved in 
mutarotation catalyzed by bases while the constant re- 
ported for glucose is an averaged constant (K,K2K;K4K¢)!", 
where subscripts refer to the positions of the five exchange- 
able hydrogens in the glucose molecule. 

When ops is plotted against P, as calculated from 
K=0.83, Eq. (2) is obeyed within the experimental error. 
The single experimental value for the isotopic exchange 
with tetramethyl glucose must be confirmed before the 
hypothesis advanced may be accepted as conclusive. 
The correct explanation may lie between the two extreme 
views. 

W. H. HAMILL 
V. K. LA MER 
Fordham University, 
New York City. 


Columbia University, 
New York City. 


1 Hamill and La Mer, J. Chem. Phys. 2, 891 (1934). 
2 Hamill and Freudenberg, J. Am. Chem. Soc. 57, 1427 (1935). 





